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1. INTRODUCTION 


The boundary value and expansion problems for the equation of the nth 
order with boundary conditions at two points have been studied by Birkhoff.t 
Béchert has suggested the generalization of these results to the equation 
with auxiliary conditions at more than two points. Such generalization of 
the essential properties of the differential system has been carried out by the 
author, and in this paper is given the proof of the convergence of the expansion, 
which may be studied quite independently of the other results. The formal 
development of the boundary problem and a more detailed discussion of the 
form of the series will be presented in other papers. 

The differential equation is taken in the form 


n—2 


+ * + Pr goa + °° + Pn U + ru = 0, 


d” u 
dx” 


in which po, p3, -**, Pn, are functions of the real variable x, continuous 
together with their derivatives of all orders in the closed interval (a, b). 
The auxiliary conditions apply ‘to k points of the interval (a, b), of which 
the first and the last are the ends of the interval. These points, arranged in 
order of ascending algebraic magnitude, will be denoted by 


G1 = @, d2, M3, +++, Ap-1, a = D. 


* Presented to the Society, April 29, 1916. 

t Birkhoff, Boundary value and expansion problems of ordinary linear differential equations, 
these Transactions, vol. 9 (1908), pp. 373-395. See also Tamarkine, Sur quelques 
points de la théorie des équations différentielles linéaires ordinaires et sur la généralisation de la 
série de Fourier, Rendiconti del Circolo Matematico di Palermo, 
vol. 34 (1912), pp. 345-382; and Birkhoff, Note on the expansion problems of ordinary linear 
differential equations, Rendiconti del Circolo Matematico di Palermo, 
vol. 36 (1913), pp. 115-126. 

t Bécher, Boundary problems in one dimension, Proceedings of the Fifth International Con- 
gress of Mathematicians, Cambridge, 1913,vol. 1, pp. 163-195; see p. 166. 
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The auxiliary conditions may then be written 
Wi (u) = Wi (u) + Wei (ue) +e) + Wis (u) = 0 (6 =1,2,-+-,n), 
in which 
W ji (%) 
and aj, ai, @;, +++, af”, are constants, real or complex. 

It is possible to define explicitly in terms of the auxiliary conditions and a 
fundamental system of solutions of the differential equation, a function 
G(x, 8;), having the properties of a Green’s function of the system. These 
properties will be discussed elsewhere with the exception of a single one, 
which forms the subject of this paper, namely that 


aj u(ajs) + aj; wu’ (aj) + aj, ul’ (aj) + eee + of? u™ (4;), 


Il 


a b 

(1) f(«) = tims G(2x,8;rA)f(s)dsdr,, 

paged, we ae A 
where I is a contour in the \-plane enclosing the first r poles of G(x, 8;X), 
and f(x) satisfies certain conditions of continuity, in some cases of the sort 
usually imposed, in others more restrictive, according to the nature of the 
auxiliary conditions. The auxiliary conditions themselves are subject to 
restrictions beyond those found necessary in the two point case. This integral 
may be regarded as the sum of the first terms of an infinite series of residues, 
corresponding to the series of orthogonal functions in the familiar cases, and 
these residues can in many cases be expressed linearly in terms of the char- 
acteristic functions of the differential system. Since the integral will be con- 
sidered here on its own merits this interpretation will not be considered 
further. 

2. PRELIMINARY FORMULAS 


A fundamental system of solutions of the differential equation may be 
chosen analytic in A. Let yi(a,A), yo(@,A), -*+, Yn(X,A), be such a 
system. Then let D(A) be the determinant which has W;(y;) for the ele- 
ment in the ith row and jth column. Its roots are the characteristic numbers 
for the system, but will concern us only as points where an expression having 
D(X) in its denominator ceases to be analytic. Let Di(s,) be the deter- 
minant which has the (n — 7)th derivative of y;(s) for the element in the 
ith row and jth column, and let D.(2,s;X) be this determinant with the 
elements in the first row replaced by the y; (a2) with corresponding subscript. 
Then let 
(2) G(x,s;}) = +34D/D,, +ife>s, —ifa<s. 


If now N(x, 8;X) is the determinant formed by bounding D(\) on the 
top by the y; (2) with corresponding subscript, and on the right by the W; (G) 
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with corresponding subscript, and filling in the upper right-hand corner with 
G(x,8;X), then the function G(x, s;X) mentioned above is given by the 
formula 


(3) G(x,8;4) =(-—1)"N/D. 


The proof of the convergence of the integral is based on the use of the 
following asymptotic solutions of the differential equation, due to Birkhoff.* 
Make the substitution \ = p". Then on any region S of the p-plane, defined 
by the relation 


lr (l+1)7 
—=—sites , 


nl n 


in which / is a positive integer less than n > 1, there exist » linearly inde- 
pendent solutions, analytic in p, and such that 


Ey 
“ 1 ppi(z—a) _—* 
yi(x, p) =ui(rt,p) +e r “om 
dy; ‘ai du; 4. epwilz—a) Ey 
(4) dx dx "lie (¢=1,2,---,n), 
n—1 4,. n—l i 
d™ y; au; 4. eo(e—a) _£, 
dy" _ dvr € p™ ati ’ 
where 
Wii (x) Uim (x) 
: — i im 
ui(2, p) = eens of 14 p oes + p™ " 


Here m is any positive integer, the w’s are the roots of the equation 
w"+1=0, 


the functions u;; (x) are finite and continuous for all values of x in the interval 
(a,b), and E is a generic notation for functions of p (and other variables) 
uniformly bounded for large values of p. In part of the work it will be neces- 
sary to use a more explicit form of the functions w;; (a), due to Milne, viz., 
ui;(x) = ;(x)/(w;)’. The functions ¢ are then the same for all the n 
solutions y;(2,). 

By the substitution \ = p” the d-plane goes over into two adjacent S 
regions, S; and S.. The integral {1) becomes 


1 b 
vif | ro G(x, 8; p")f(s)dsdp 
——— 


(5 a . 
ton f f np” "G(x, 8; p")f(s)dsdp, 
ice y2°7a 


* Birkhoff, loc. cit., p. 381. 
t Milne, Note on asymptotic expressions in the theory of linear differential equations, Bulle - 
tin of the American Mathematical Society, vol. 23 (1917), pp. 166-169. 
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in which y; and y2 may be taken as arcs of the same circle with center at the 
origin. 

Before substituting these solutions in the auxiliary conditions, the condi- 
tions are normalized as follows: 

First, the number of conditions of order n — 1 at the points a and b is re- 
duced, by means of linear combinations, to two at most; then, the remaining 
conditions of order n — 2 at these points are reduced to two at most; and 
so on. 

It is convenient to introduce Birkhoff’s notation, by which [w] represents 
an expression of the form 

a b_e¢ E 
ee $345 F ++: +=, 
eS a p 
in which w, a, b, etc., are independent of p but need not be constants, and E 
is a function of p and other variables bounded for large values of p. 

The result of substituting the asymptotic solutions into the auxiliary con- 

ditions is 


Wii (ys) = (pw) Lau], 


Wei (ys) = (pw;)!™ [agi] ere, 


Wri (ys) = (pw;)™ Lani] errer-. 

Because of the normalization of the auxiliary conditions hy; = hi; S hj, 
d ] 

j =2,3,+--,k —1. So let us set hy; = hy; = hj, and we have hj; = h; 


= his1, and no three successive h;’s are equal. 
In the odd order case let n = 24 — 1. Then the indices 1 to n can be 
determined in such a way that for any p in S 


R(pw,) S R( pwr) S R(pw;3) S --- = R( pwr), 


where R(w) indicates the real part of w. Hence for S 





R(pw;) < 0 (¢=1,2,-++,n2—1), 
> 0 (t=n+1,0+2,---,n), 
= 0) (i =nand arg» = "+? y), 
and it follows that 
Wi (ys) = (pw;)"lar] (j=1,2,--+,u—1), 
= (pw;)™[ ax; ] er (j=u+1,u+2,-++,n), 
k-1 


( pw; yh [ a; | + ( pw; yh [ Qh ] epm(b—a) abe > ( pw; yh [ Qn; ] ePes(u-a) 
(j =u). 


t=2 
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In the even order case let n = 2u. Then the indices 1 to n can be deter- 
mined in such a way that for any p in S 


R(pw;) < 0 (¢=1,2,3,---,n—1), 


> 0 (i=n+2,---,n), 
: . lox +1 
= () =y,u+1, and either arg p = (OF are p = -w), 
t 


and since — w, = W,41, we have 


Wi(ys) = (pw;)* [ax] (fj =1,2,+++,u-1), 
= (pw;)"[ az; ] er (j=u+2,---,n), 
k-1 
a (pw, yr [ a1; ] + (pw, ys [ Qi ] eh (b—a) + 2. (pw, yhei [ Oy; ] ¢P'u (ar—a) 
t=2 
(j =4), 


( — pw, y [ a1; | + ( — pw, ym [ aK: ] e7 Pu (b—a) 


k-1 
+ 2 (= pw, YM Lons]ermeor? = (j =n +1). 
l=2 


3. THE FUNCTION D ()) 

Preparatory to taking up the main convergence problem it is necessary to 
study the distribution of those points where the function D(X) vanishes or 
becomes nearly equal to zero. This study is simplified by considering separ- 
ately the cases where n is odd and where n is even. 

Odd order case. When n is odd, the substitution of the asymptotic ex- 
pressions for the W’s in the determinant D (\) gives 

n n 
D(r) = JL o™ TT A, 
: i=1 3 3=j=yn+l 
in which A is that determinant which has for the element in the ith row and 
jth column, [ a1; ] wj forj = 1,2, --+,u — 1; [ax] wi forj =np+1,44+2, 
-, n, and finally 


k-1 
B; = wih (Leas) + Lars]? + D0 (pw, [cs] erator), 
t=2 
for 7 =p. Let v be the highest power to which pw, enters in any B;; then 
A can be expanded in the form 


k—1 
(6) A = [6:1] + [Ox] eo 9 + (pw, )”) 2 [01] ema . 
l=3 


We shall now prove the 

THEorEM. If in (6) neither 0; nor 0, is zero, and if for v= 1 neither 0, 
nor 0,1 1s zero, then D(X) has an infinite number of roots, and there exist two 
absolute constants w and Q such that the number of roots of D(X) in any circle 
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about the origin of radius r is less than Qr and greater than wr and, finally, for 
values of \ uniformly away from the roots, D(X) remains uniformly away from 
zero. 
This last statement is made more explicit by inequalities (16) and (21) below. 
We assume, then, that neither 6, nor 6, is zero. The number v may be 
negative, zero, or positive. If it is negative A can be written more simply 


[ 0; ] + [ Oy ] ey (b—a) . 
Birkhoff* has shown that the roots of this are asymptotically equal to 
l ( 91 ) 
% Mar V— 1 ~~ * O } 
pt w, (b — a) w,(b—a)’ 


and that A remains uniformly away from zero if p is uniformly away from a 
root. 
When + is zero, A becomes 


k-1 


(7) [ 6; ] -- [ Oy ] ey (b—a) +. ue [ 0, ] ew, (ar—a) , 


and now at least one other @ in addition to @, and 6 is assumed to be different 
from zero. In treating this case the distribution of roots of the function 
obtained by suppressing the brackets is: determined, and then the roots of A 
are proved asymptotically equal to these. 

It is convenient to make a change of variable. Let pw, (a2. —a) =z, 
and divide through by 6,. The equation to be considered is, then, 


(8) f(z) = 1+ Qoe + Qice™ + Qre* + --- + Qr2e* =0. 


‘The Q’s are complex constants, and in particular Q,-2 and at least one other 
are different from zero. The c’s are real constants such that 1 < ¢c; < ¢e 
< +++ <epe. Let z =2+ yi, and let R(@) denote the real part of ¢, 
and I (¢) the coefficient of the pure imaginary part. Let ¢ be an arbitrary 
small positive quantity, and choose X so large that forx = — X, f(z) differs 
from 1 by less than e, and at the same time for x = X, f(z)/e** differs 
from Q,~2 by less than e. The region between x = X and « = — X will 
be called the region R throughout this paper. We shall suppose that « is 
chosen small enough so that f(z) does not vanish outside of R. 

Now f(z) is analytic in any finite region. Hence in any interval, no 
matter how small, along the y-axis, it is possible to find a line y = constant 
along which f(z) does not vanish. Let y = Y; and y = Y2 be any two such 
lines. The number of roots of f(z) inthe rectangler = X,y = Y1,2 = —X, 
y = Y2 is determined by following the change in arc f(z) as z describes this 


* Birkhoff, loc. cit., pp. 383-5. 
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rectangle. Now 


I(f) 


are f(z) = sin™ fl” 


and for y constant J (f) has the form 
T(f) =dhe? +dje™ + dee" +--- + dp2e*", 
in which the d’s are real constants. The finite roots of J (f) are the same 
as those of 
To(f) = do + dy e&™* + dy O™ 4+ +++ + dyed, 


But between any two roots of this is a root of the derivative 


I) (f) = di(ex — 1) e&"* + dz (e. — 1) e&"™ 
_ one = dy—2 (cr—2 ae 1) efo-2— — 


while its roots are the same as those of 
Ti (f) =dy(¢1 — 1) +2 (e2 — 1) &-™ + +++ + dys (Cre — 1) 9”. 


The argument is repeated till finally J;,2 is reached, which has only one term 
and does not vanish for z finite. Thus J (f) has at most (k — 2) finite roots. 
Should any of the J’s have multiple roots the argument must be slightly 
changed, but the result is, of course, the same. 

Now as z traces the line x = — X, f(z) remains very nearly equal to 1, 
so that the resulting change in arc f (z) is less than some small positive quantity 
mi which approaches zero with e€. As z traces the line x = X the change in 
are f(z) differs but slightly from cy,-2(Y2 — Y1), as is at once apparent 
from the formula ) 

+ are ac 


Cy2% 


arc f(z) = arce 


In other words, the change in are along this line can be written 
cr-2(Y2 — ¥1) +m, 


where 72 is a small positive quantity which approaches zero with e. Finally 
along the lines y = Y; and y = Y, the function sin {are f(z)} vanishes at 
most k — 2 times, whence arc f(z) increases by less than (k — 1)7 + 73 
and (k — 1) + m respectively, where y3 and 74 are again small positive 
quantities which approach zero with e. To sum up, as z traces the rectangle 
the resulting change in are f(z) is not greater than 


cr-2( V2 —¥;)+2(k-—1)r+7% 
nor less than 
Che(Y2—- V1) —2(k-1)r— (n=m4+m4tn4+m), 


if Y; — Y; is large enough to make this expression positive. Since 7 can be 
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made arbitrarily small by taking ¢ small enough, it can be dropped from these 
inequalities. So finally in any interval of R of length / the number of roots of 
f(z) is not greater than cy. 1/27 + k — 1, nor less than cy_21/2r —k +1, 
if this number is positive. In particular, by making / very small, one sees 
that a multiple root of f(z) cannot be of order higher than k — 1. 

It is now necessary to prove that if z remains uniformly away from a root 
of f(z) then f(z) remains uniformly away from zero. This is proved by 
means of the following auxiliary theorem: 

THEOREM. (Given f (2,21, %2,°**,2%n), continuous in all of its arguments 
and analytic in z for a; = x; = b; and z in a closed finite region S of the complex 
plane, and such that for no set of values of the x’s is the number of roots of f(z) 
greater than a given constant N; if for any set of values of the x’s, z 1s a point at a 
distance greater than 6 from a zero of f(z) and from the boundary of S, then 


f(z, a1, 22, o++, Qn) (ZX, 


where d is a real positive constant independent of the x’s. 

It is no restriction to assume that 6 is small enough so that for each set of 
values of the x’s the corresponding region of variation for z, defined by the 
hypotheses, actually contains points. This follows from the assumption 
that for any set of values of the x’s the number of roots of f(z) is less than NV. 
Now for any particular set of values of the x’s there is a \, such that for these 
x’s and for z restricted to its corresponding region of variation, 


f(z, %1, 22, +++, an) |ZA, 


while there is a value of z for which the equality holds. We have thus defined 
a function of the 2’s, \ (a1, 22, -++,2n), which, by means of the easily proved 
fact that the zeros of f(z, 21,22, ---,2n) are continuous functions of the 
x’s, can be shown to be continuous in the closed region a; = 2; = b;. Since 
(a1, %2, ***,%n) is continuous in a closed region, it has a lower limit which 
it attains. Call this limit \. Then for some value of z and the 2’s, |f| =. 
Hence \ is not equal to zero, for f does not vanish in any of the regions for 
which the function of which \ is the lower limit was defined. Hence the 
theorem is proved. 
First apply this theorem to the function 


f(z, b1, 2, +++, bee) = 1 + Qoe? + Qi ete +--+ + Qe Crt, 


The region S can be that portion of R for which 0 = y S 27, augmented by a 
strip of width 6 surrounding it. Obviously the breadth of R, that is, the 
value of X , depends only on the Q’s and not on the ¢’s. It is also clear that 
the function satisfies the conditions of the theorem. Let X indicate the 
number determined by the theorem corresponding to 6. 
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Consider, then, the function (8) in R. Any interval of R, 
Qe Sy S2(14+1)z, 


can, by a transformation 7 = y + 2lz be brought to the interval 0 = y = 27, 
and f becomes 


f(z 4 2I77 ) = 1 + Qo e? + 01 efit—Uaeit + oa em + Q,-2 eckst— tego 


But, without changing f, each expression 2/rc; 7 can be replaced by a jt 
where 0=@=27. So f(z+ 2lri) has the form above, and for z at a 
greater distance from a root than 6, |f(z+ 2lr)|=X. By giving / all 
integral values, all of R is covered by these intervals, and so finally, for z in R 
and at a distance from a root of f(z) greater than 6, |f(z)|=A>0. 

We can now prove that the roots of (7) are asymptotically equal to those 
of (8). First it is clear that by choosing ¢ small enough we can be sure that 
all the roots of (7) lie within the region R. Let us denote (7) by ¥(p) and 
(8) by f(p). Thenin R 


¥(e) = OF (e) +=, ¢ 


from which it is obvious that the roots of ¥(p) are very nearly equal to those 
of f (p) for large values of p. Now 


E 
) Vio) =Hf(0)(1+55705)5)> 


whence 


E 
arc (p) = arc 0, + arcf(p) + are (1 os Foy) ; 

Let 6 be an arbitrary small positive quantity. Then according to the 
previously determined distribution of the roots of f(p), it is possible to 
enclose these roots in groups of less than k in simply connected regions of 
major diameter less than k6, and such that no point of the boundary is at a 
distance from a root of f(p) less than 6. On the boundaries of these regions 
we have then |f(p)|2=X, where is determined by the theorem above. 
Hence for large values of p, the expression E/(6:f(p) p) is very small on these 
boundaries, and so the resulting change in are (1 + E/(6:f(p),)) asz traces 
the boundary of a region is zero, and the change in arc ¥(p) is the same as 
the change in are f(p). So the number of roots of ¥(p) in any such region 
is for large values of p the same as the number of roots of f (p) in it; that is, 
the roots of ¥ (p) are asymptotically equal to those of f(p). 

From this together with (9) can be deduced the fact that for large values 
of p in R the function ¥(p) remains uniformly away from zero when p is 
uniformly away from a root of ¥(p). Outside of R the function y(p) of 
course remains uniformly away from zero. 
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It is convenient to consider next the case in which A reduces to 
[61] + [0] er + (pw, )” [02] err , 


which, in particular, it always does if k = 3. It is assumed that none of the 
6’s are zero, and thaty 21. Then 
A = [62] eu“ Ao, 
where 
(10) Ao = (pw, )” + [Jeu + [Bq] rar , 
and in this 


[65] = 1) 
“" [62]? 
It is necessary to consider separately the part of S in which R (pw, ) is positive 
and the part in which it is negative. 

R(pw,) =0. It will now be shown that in that part of S where R (pw, ) 
= 0, the roots of Ap are asymptotically equal to those of 


f(p) = (pwoy)” + 0, eM, 


To find the roots of f(p) set pw, =x + yi = re*, and — 6, = re. 
When f() is equated to zero the second term may be transposed, and then, 
when absolute values and angles are equated, the two following equations 
result: 


r= r,e7r , vp = da — Y(d2—a)+2mer (m=0,+1,+2,4+3,---). 


The zeros of f(p) are the points of intersection of these two curves. Since, 

in that part of the p plane under consideration, x = 0, y > 0, the first curve 
p 

is equivalent to 


° 


a Ziv p—2x(aq—a) /v __ 2 
y — Vr? € x ’ 


and this has but one branch which is in character much like the exponential 
curve, and which for large values of p approaches parallelism to the line 
R(pw,) = 0,i.e., 

lim (dx/dy) = 0. 

p=@ 


It is also easily shown that 
lim (2/y)? =0, 
p=a 


which is to say that for points on this curve 


lim @ = 7/2 + 2lnr (l=0,+1,+2,+3,---). 
p=e 


When this is substituted in the second curve it is found that the roots of 
f(p) are asymptotically equal to points whose y coérdinates are given by 


_ ba + 2mn — v (30 + Ar) 


U 
y ae ‘ 
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and since v is an integer this is equivalent to 
ga + 2ma — dvr 


1 
y ay — a 


So the roots of f(p) lie along the curve 

(11) | (pw, )”| = |B0 Pra | 

at intervals whose lengths approach 27/(a2 — a) as p becomes infinite. For 
large values of p these roots are all simple, for it is obvious that f (p) and f’(p) 
cannot both vanish for the same value of p when p is large. 

In order to prove that in that part of S where R(pw,) = 0 the roots of Ag 
are asymptotically equal to those of f (p), it is necessary to show that |f (p)| 
remains uniformly away from zero when p is uniformly away from a root. 
For this purpose let 7 be the region included between the curves defined by 
the relation 


|( or, )”| = [_ e-raioraree |, 


Let U be the region between 7 and the line R(pw,) = 0, and let V be the 
remainder of S for which R(pw,) <0. Choose o as any convenient positive 
real number. The following inequalities are immediate: 


for U, | (pw, )” | = e7 |0, ce evultz— | , 


4 | —g | — ag— | 
for V, | (pw, )” | S e~? |0, eMule | 5 


from which it follows that 
for U, lf(p)|=(er —1)|0.e°%"= |, or (1 —e%)|(pw,)’|, 
for V, If(p)|= C1 — € 7) | 0, Pu | | 
For p in T it is necessary to proceed differently. Let po be any root of 
f(p) = 0, and call those points of 7 which are at least as near to po as to 
any other root of f(p) = 0, the domain of pp. Now the strip T is asymp- 
totically of constant width, and it is easily shown that this width is 20/(a,.—a). 
On the other hand the distance between two roots of f(p) = 0 was shown to 
be asymptotically equal to 27/(a, — a), so that finally if po is far enough 
out the distance of any point of the domain of po from pp will be less than 
2(¢0 + 7)/(a2—a). In this domain f(p) can be written 


f(p) =f(e) —f (0) = (pw,)” — (po Wy)” + 04 CP" — Og EP ule) 
(13) = w',(p — po) (p> + p”? po +--+ + po?) 


+ 6. oP y (42-2) (1 <se e6P—Po) (42-4) ) . 


(12) 


In 7 the absolute value of the first term is less than 


2(e+n) | e+e) 
| | ’ 


dg —a ado — a 
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and for sufficiently large values of p this is less than a constant times | po|". 
Moreover, since in 7’, | pw,|” is between two constant multiples of |¢~°"#“— |, 
it follows that the first term is less than a constant times 


|e Kv 1iv low, (e2—4) | : 


In the parenthesis in the last term let 
(p — po) wv, = r(cosd?+isin®g), where —-rtS¢S7. 


The parenthesis becomes 
( 1 _— gi’ s-*) (eos +i sip ¢) ) . 


which can vanish only at the points given by 
2lr 


r(ad2 — a) 


cos @ = 0, sin @ = 


and the point pp. In the domain of pp there is only one of these points, namely 
po itself, and for that portion of the domain whose points are at a distance 
not less than 6 > 0 from pp the absolute value of the parenthesis as a function 
of r and ¢ is continuous and positive in a closed region. Hence it has a 
positive minimum which may be denoted by 7. The last term of (13) is then 
greater than a constant times 


| — pie, (42—4a) 


é Lal 


for that portion of the domain of po at a distance from po greater than 6. 
This, combined with the inequality on the first term, gives the result that 
for large values of p in T and at a distance from a root of f(p) = 0 greater 
than or equal to 6 > 0, 

(14) If (o)| > lew" | Ko, 


where Ko is a positive constant. Similar inequalities have been proved for 
any p in U or V, so that this, with a new constant K,, may be assumed to 
hold for all that part of S for which R( pw, ) S 0 for large values of p at a 
distance from a root greater than or equal to 6. 

From (10) is obtained the relation 


‘ E 
Ay — f (p) — ep (t2—-4) [ 6, | ePu(o—42) . 
which can be factored into 


— pty (t2—a) | 4, ] oP u(o—a2) ) 


* on Ee 
(19) ao = f(0)(1 +: fip) * fie) 


Consider the portions of 7’, U, and V , bounded by two circular ares |p| = p’ 
and |p|= p”, (p’ <p”), such that no point of either arc is at a distance 
less than 6 from a root of f(p) = 0. Call sucharegion S. It will be shown, 
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by considering the change in arc Ao as p traces the boundary of S, that the 
number of roots of Ap and of f(p) in S is the same, provided the innermost 
boundary of S is taken far enough out. Consider first the last term of the 
parenthesis in (15). Throughout the region S the numerator is bounded. 
In U the denominator, f(p), can be replaced by (1 — e~’) (pw, )” from 
inequality (12), and so for large values of p in U this last term is small. From 
the second of the inequalities (12) it can similarly be shown that this last 
term is also small for large values of p in V , while from inequality (14) the 
same result is obtained for the boundary lines of S where they cross 7’, pro- 
vided the innermost boundary of S is far enough out. On the other hand, 
from inequalities (12) and (14) the expression ¢~°”«“*~” /f(p) is bounded on 
the boundary of S, and when it is multiplied by E/p the result is small when 
the innermost boundary of S is far enough out. Under these conditions, 
then, the parenthesis in (15) can be written as 1 + 7, where 7 is small enough 
so that the are of 1 + 7 will return to its original value when p traces the 
boundary of S. It follows that the change in arc Ao is the same as the change 
in are f(p), so that Ap and f (p) have the same number of roots in Ss. 

Now exactly the same result is obtained if p traces the circumferences of 
circles of arbitrarily small radius about the roots of f (p), and this, combined 
with the previous result, shows that in that part of the region S for which 
R (pw, ) = 0, the roots of Ap are asymptotically equal to those of f(p). 

An important inequality for A can be deduced from formula (15). For p 


large and at a distance from a root greater than 6 inequality (14) holds. For 
T and V the exponential can be replaced by a constant times (pw, )”, while 
for U the resulting inequality for f(p) has already been obtained in (12). 
The parenthesis in (15) can be made nearly equal to 1, and in particular, 
greater than 1/2, by taking p large enough; so that the inequality is obtained: 


|Ao| >| pw, |” Ke, 


in which Ke is a constant depending only on 6. The result of applying this 
to A is (see (10)) 
(16 JA] > [Cory oem? | Ke, 


in which K; is a constant depending only on 6. This holds for p in that 
part of S where R( pw, ) is less than or equal to zero and greater than some 
fixed p, and at a distance from a root of f(p) at least as great as 6. Equally 
important is the result, easily deducible from the same formulas, that under 
the same conditions A is greater than a constant depending only on 6. 
R(pw,) 20. By entirely similar reasoning it is found that for that part 
of S where R( pw, ) = 0 the roots of A are asymptotically equal to those of 


(pw, )” + 0 ero, 
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| | w, (b—ae) | 
| ( pw )” | = | eP"! 2) | 


at intervals whose lengths approach 27/(b — a2), and which are all simple 
for large values of p. Inequalities on f (p), Ao, and A, similar to those above 
are found also in this case. 

Finally the odd order case with k points is in general, i. e., when certain 
@onstants do not vanish, merely a combination of the cases already con- 
sidered. The characteristic equation A = 0 is as given in (6). The coef- 
ficient of (pw, )” can be brought to the form (8) by dividing by e?“«“?-® 
and setting pw, (a3; — a2) =z. It is now assumed that none of the con- 
stants 6;, 02, 0,1, 0%, is zero. It follows then, that the roots of this coef- 
ficient are all in a certain region R lying along the line R( pw, ) = 0, and that 
for p at a distance from a root of this coefficient not less than 6 the coefficient 
is in absolute value greater than |e?"«“*-|K,, in which Ky, is a constant 
dependent only on 6. For convenience let this coefficient be denoted by A. 
Then the inequality is 
(17) 
and A is 


[Ai | > | eh nlee—*) \K4, 


A= [ A; ] + [ 0, | eh" u(o—a + (pw, y¥ A). 


This can be written 


6 6, oP, (b—a) 
(18) A= (pw) (1+! il+[ jon). 


(pw, )” Ai 


Now in R the exponentials are bounded. Hence, if in R any region be drawn, 
such that no part of the boundary is at a distance from a root of A, less than 
5 > 0, then on the boundary A; remains uniformly away from zero, and if 
the region be far enough out the parenthesis in (18) has, on the boundary, 
the form 1 + 7, where |n|< 1. From this can be deduced the fact that the 
roots of A in R are asymptotically equal to those of A;. 

Outside of R for that part of S in which R(pw,) < 0 the roots of A are 
asymptotically equal to those of 


6 
f(p) = (pw,)” +a err) = 0, 


and for that part of S in which R(pw,) > 0 they are asymptotically equal to 
those of 


0; 
(pw, )” +a eter) = Q, 
O.-1 


To prove the first of these statements, write 


[A:] [ 0, ] eP™uo-@) 
— gh y(a2—a) paces 
Ai 6, 


A 
— =fle) + A, 


Ai 
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k—1 
[0] — A »» [ 0; ] ep" u(a-2) 


j=3 [ 0, ] erm noo) 


(19) =f(p) +> Alb ~ Ps 


k-1 
1] py (4-22) 
[0] 0: 2.185] ¢ » len) 


“101 +5a7@) ~ aT Aif(e) 

The roots of f(p) have already been found, for it is the f(p) of (11). To 
prove the roots of A asymptotically equal to those of f(p), it is obviously 
sufficient to prove that the parenthesis above has the form 1 + 7, for p on 
the boundary of any region outside of R and in that part of S where R ( pw, ) 
< 0, and such that the distance of any point of the boundary from a root of 
f(p) is not less than 6 > 0, provided that this region is taken far enough out 
from the origin. For such values of p the inequalities (12), (14), and (17) 
hold, and they will be used to show that the parenthesis has that form. From 
(14) and (17) is obtained |f(p)A:|> Ki Ky, and since the symbol [0] is 
the same thing as E/p it is evident that the term [0]/A; 62 f() is small for p 
large. In the next term the same inequality on f (p) A; is used for 7 and V, 
since in these regions, for p large, R( pw, ) is large and negative, and the 
numerator, and so the whole term, is small. For U the inequality |f(p)| 
= (1 —e*)|(pw,)”| from (12) is used for f(p). The remainder of the 
term may be written in the form 


Eom? ul e2-2)— (4-42) } = Ee? ula2—2) 1 (aaa) / (42-2) ] P 
Now in U the inequality 
| (pw, )”| > e” Oa EP yp (42—2) | 
holds, from which it follows that 
’ 


|e PM nCoe—*) [1 Cosa) /(o2—-#))] | = | (pw, y7—(as—a2) (42—4)] \K ; 


where K is an absolute constant. The exponent of this power of p is cer- 
tainly less than v, so that because of the factor |(pw,)”| from f(p) in the 
denominator, the term approaches zero as p becomes infinite. By reason of 
inequality (17) the last term of the parenthesis can obviously be written 
E/f(p). InT and V inequality (16) shows that this is small for large values 
of p, while the inequality |f(p)|= (1 —e*)|(pw,)”| from (12) serves 
the same purpose for U. Thus the parenthesis has been shown to be of the 
required form for the values of p in question, and so the roots of A are 
asymptotically equal to those of f (p) for that part of S outside of R for which 
R(pw,) < 0. 
The proof for the remaining part of S on the other side of R is similar. 
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It is to be noted that the above reasoning applies to the whole of that part 
of S in which R(pw,) S 0, even including the strip R, provided that p is 
at a distance greater than 6 from a root of f(p) or A,, and is large enough. 
Hence from (19) can be obtained the following inequality which holds for all 
p’s in absolute value greater than a certain p and at a distance from a root of A 
greater than 6: 


(20) |A| >|Ai|-|f(p)|Ks, 


in which K; is a constant depending only on 6. 

In particular this shows, by reason of inequalities (14) and (17), that if in 
this region p remains uniformly away from a root, then A remains uniformly 
away from zero. 

Another inequality for A can easily be deduced from the above. First 
replace |A;| by |e#“*-®|K, from inequality (17). Then for the region U 
replace |f(p)| by (1 — e~”)|(pw,)|” from (12), and for the regions 7’ and V 
replace it by |e~?"«“?-®|K, from (14). But in the regions T and V the 
exponential can be replaced in this inequality by a constant times |(pw,) |”, 
so that finally the following is obtained: 


(21) |A| >| (pw, )” e"#= | Ke, 


in which Kg is a constant depending only on 6, and the inequality holds for 
all values of p greater in absolute value than a certain p and at a distance 
from a root of A greater than 6. 

Even order case. When n is even the substitution of the asymptotic ex- 
pressions for the W’s in the determinant D () gives 


n n 
D (r) = II pl Il eps(b—a) A, 
ix1 j=ut2 
in which A is that determinant which has for the element in the 7th row and 
the jth column, [ a; ] wj‘forj = 1,2,---,u—1,and[a,;] w}forj =p +2, 
uw+3,---,n, and 


k-1 
Di=w' (Lari) + Lars Je" + (prey) [ous ]ema®) (j=), 
(21a) me, 
H;=(—w,)™ (Lars) +[ ans eo" +2 ( — pw, )—™[ oxy; Jer ) 
7 (j =#+1) 
where n = 2y. 
It is evident from the definition of the bracket symbol that in general this 
determinant expression for A can be expanded in the form 


(22) A = P;(pw,)” + P2(pw,)” + Ps, 


where the highest and lowest order terms of the P’s are 
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P;: [O11 ] oP u(te-1—42) and [ O15 | ey (tk-1—72) : 


P2: [0.,] e""°-™ and [O.2]e °°"), if b—as>ay1—a, 
(23) ' ' 
or [63,]e?"#" and [63.)e7°"#%) if = b—ax<ay_1—a; 
P3: [Oy Je" and [642] ec P%#- , 


The number + is the exponent of the highest power of p appearing in a D; 
or H; e 

We shall now prove the 

THEorEM. If none of the @’s in (23) are zero, then D(X) has an infinite 
number of roots, and there exist two absolute constants w and Q such that the 
number of roots of D(X) in any circle about the origin of radius r is less than 
Or and greater than wr , and finally for values of \ uniformly away from its roots, 
D(X) remains uniformly away from zero. 

This last statement is given more explicitly by inequality (25). 

Should the number vy be negative, A reduces formally to Birkhoff’s even 
order case,* and the problem of distribution of roots will not be considered 
farther here. If v is zero, A becomes merely a sum of exponentials which 
when multiplied through by the non-vanishing factor e?”#@>® all have real 
parts of the same sign, and so are reduced to the form (8) already treated. 
In either of these cases the formal expansion (22) is not needed. But when v 
is positive we make use of this expansion, which tacitly assumes that at least 
two of the D;’s and H;’s actually involve this power of p. This is, of course, 
equivalent to assuming that the highest order derivative in the auxiliary 
conditions at points other than the end points appears in at least two con- 
ditions. We also assume that none of the @’s in (23) are zero, which is to 
say that certain of the determinants from the matrix of coefficients of the 
auxiliary conditions must not vanish. The method here used can be applied 
equally well with very slight modifications to many of the cases thus excluded, 
but a detailed study of the more specialized cases is beyond the scope of this 
paper. 

The more general even order case, then, is treated as follows: The function 
P, can be reduced to the type (8) by multiplying through by e?”#“*>-) and 
making the proper substitution. Its roots, then, all lie in a region R, and the 
proof that the roots of A in R are asymptotically equal to those of P; is im- 
mediate. To find the roots outside of R write 

> 
p= (een) + (pw) GE +B. 


The line R( pw, ) = 0 is now not the bisecting ray of each of the S regions 


* Birkhoff, loc. cit., p. 387. 
Trans. Am. Math. Soc. 28 











432 CHARLES E. WILDER | October 


as in the odd order case, but the line of demarcation between two of them. 
Hence the region R overlaps the t wo S regions into which the \-plane is carried 
by the transformation \ = p”. Consider first the region S; in which R (pw, ) 
<0. The above expression can be written in the form 


A = . [ As2 ] ——w om} , \v [B42 —pw avo | 
PT { (ee) T F612] sei (ew) ‘ee 


+ (pw, )” Qi + Qe. 


Suppose that a; — a is not equal to b — a,_; and, for convenience, let az — a 
> b —a,y-1. Then the term of highest order in Q; is 


[0] Pep (O41) , 
and in Q2, 
[ 6] eM pha) a 


The expressions Q; and Q» are merely those terms which must be added to 
the product of the two factors in order to make the result equal to A/P;. 

The roots of the two parentheses are found by the methods already de- 
veloped, and it can be shown by means of inequalities of the type (12) and 
(14) that the roots of A in S; and outside of R are asymptotically equal to 
those of the two parentheses. 

Similarly the roots of A in S, and outside of R are asymptotically equal to 
those of two parentheses of the same type, but with positive instead of nega- 
tive exponents. 

The following inequality is easily found for large values of p which are in 
S, and at a distance from a root of A greater than 6: 


| | [ Os2 —pw,,(ae—a) | [ O42 
|A| >| (pw, )” +e peple—*) |. | ( pw, )” +194 


in which K; is a constant depending only on 6. If 


e7Py(>—a-1) | ‘ | Pi | K, ’ 


A= Aer™n(o—2) | 


the inequality for A is 
|A| > | Pi Pl t-1—42) | " | (pw, )” Pu (%2—4) 


(24) [ As2 ] [ Ase ] 
+ T2] [52] 


It can also be proved that for large values of p in S; which are at a distance 
greater than 6 from a root of A, the absolute value of A remains greater than 
a positive constant depending only on 6. 

When az — a = b — a,_; the two systems of roots lie along two curves of 
the type (11) in which both have the same exponential, although in general 


| ‘ | (pw, y Pu (o—24-1) + |_Kz. 
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the coefficients are different. If these coefficients should have the same 
absolute value the roots all lie along one curve, and if they are identical the 
roots are all double. An inequality like (24) holds in this case also, the only 
difference being in the 6’s and the value of Kz. From (24) it is easy to deduce 
the inequalities 7 
B| > | Cpu) ert ota | Ky, 
(25) > | (pw,)” ea | Ky, 
>| (pw, )” eraO—-) | Kio, 


which hold under the same conditions as (24), and in which Kg, Ky, and Kyo 
depend only on 6. 


4. THE EXPANSION PROBLEM 


The expansion theorem is as follows: 

THEOREM. Given f(x), any real function which together with its first v 
derivatives (or if v < 1 with its first derivative) is continuous in the closed interval 
(a,b), then its expansion in terms of the residues of the Green’s function (3) 
converges to f(x) at every interior point of (a,b), if the auxiliary conditions 
have the following properties: 

In the odd order case 

Neither 6; nor 6; in (6) is zero. 
If v = 1 then in the coefficient of (pw, )” neither 62 nor 0;,—; is zero. 
In the even order case 
None of the 6’s in (23) are zero. 
The intervals b — az and a,-, — a are longer than any other of the subinter- 
vals of (a, b) included between any two of the points a, a2, a3, +++ ,a,~-1,b. 

The essential purpose of this last condition is merely that the second point 
from either end shall be farther away from that end than the first point from 
the other end is from thatend. The conditions on the @’s require merely that 
certain determinants from the matrix of the constants in the auxiliary con- 
ditions do not vanish. The condition on the 6’s in the even order case is 
stated in the above form in order to have it as simple as possible. If v < 1 
it is not necessary that all these 6’s be different from zero. 

To prove the theorem it is necessary to consider separately the odd and 
even order cases. But before that is done it is convenient to introduce a 
new form of the expression for G(x, 8;X) given in (2). If G be expanded 
in terms of the first row of the numerator determinant it can be written 


= 1 bad 
G(2r,s;) = +5 2 yi(#) G(s), 
in which 
(1 HOD... Yn—1(8) z). 


oF 7. ry, ees =< : nal 
(26) (2) = np” pw; (pw; )? bd id (pw;)""" " p” 


Jorv Qu p 4/7, 
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This expression for the g’s in so far as it differs from the one given by Birk- 
hoff* is due to Milne.¢ In the work to follow, if »y < 1, Birkhoff’s form would 
be sufficient, and the convergence proof is no more complicated than in the 
two point case. The complications necessary in case v = 1, are due to the 
necessity of using the more explicit formulas of Milne. 

The above form of G gives for W;(G) 


n n e 
Wi(G) = — 5 Le Waslys)Gs(s) +5 2 Was (ys) G(s) + - 
“aa . . : 

+ 5 2 Weas(ys) di (8) +5 2 Wes (ys) G5(8), 
in which the plus sign is taken for values of s to the left of the point at which 
the term in question applies, and the negative sign for values of s to the right. 

Odd order case. R(pw,) =0. The convergence of the integral (5) is 
now to be studied. It has already been separated into the sum of two inte- 
grals, one in S,; and onein Sz. In the odd order case it is necessary to separate 
each still farther into the sum of two integrals, one in that part of S where 
R(pw,) <0, and the other in that part where R(pw,) >0. The suc- 
cessive paths of integration are arcs of circles |p| = constant extending across 
both the regions S; and S2, and so chosen that no point of such an arc is at a 
distance less than a preassigned 6 > 0 from a root of D. That such arcs 
can be drawn in such a way that the sum of the orders of the roots between 
two successive ares is not greater than 2k when 6 is chosen small enough, 
follows from the distribution of the roots of A as determined earlier in this 
paper. The parts of these arcs in S; and S, will be denoted by y: and 2 as 
in (5). Let S represent either of these regions. Then let B be the point 
where the line R (pw, ) = 0 in S cuts one of these ares, and let A be the other 
end of the are that crosses the region where R(pw,) <0, and C the other 
end of the arc that crosses the region where R(pw,) > 0. The convergence 
of the integral over the arc AB will be studied first; that is, the integral 


(27) 254). f np" "G(x, 8; p")f(s)dsdp. 
Before making use of the expression (3) for G (2, 8; p") it is convenient to 
alter the last column by adding to it each of the columns 1 to » multiplied by 
the term 39; (s) with corresponding index, and each of the columns yu + 1 to 
n multiplied by the term — 39; (s) with corresponding index. The resulting 

last column is: 
1st element, yy (2) 9: (8) fory <2,or—- > yi(2)9i(s) for y>2; 


oe it=yu+l1 


* Birkhoff, lo loc. cit., p. 391. 
t Milne, loc. cit., p. 167. 
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(1 + 7)th element, 
— LD Ways) 9i(s) 


" — 2 Wai (ys) ai (8) + hi 5 ec 
goes P 
(28) p> Wei (ys) 95 (8) 

j= 


— 2 Wea, s(ys) G(s) ; 

j=ut+ = 

+ i" +2 Was (ys) 5i(8), 
Ds Wes, i (ys) 9 (8) ~ 
j=l 


in which the upper term of each brace is taken for values of s to the right of 
the point to which that term applies, and the lower for values of s to the left. 
The asymptotic solutions are now substituted in the expression G , and then 
the factor p” of D is divided into the (1 + 7)th row of the numerator, and 
the factor e?”°-® into the jth column. The factor — 1 is divided into the 
last column. Thus is obtained _ 
N 
G(z,8;p") = np" BD? 
in which N is the determinant A bounded on the top by e@-® [1] with 
corresponding index for the columns 1 to yw, and by e?”@—» [1] with corre- 
sponding index for the columns » + 1 to nm. On the right the ith row is 
bounded by F; and the upper right-hand corner is filled in with F, where F; is 


n 
? hy+ (a— 
n kl a ) os [ ani wt je" oc 
j=euti 
Ay+l o.(a— - 
ani } [ ay; wt jer"a 8) +> pe™ : 
OTe j=urtl §=3 7 h 1 7 2 
(29) » > any, whit? | ermdar—*) 
j=1 


yu 
+ Do [ores wit] rrr), 
j=! 
and F is 
n 


B 
+e [w;], fors <2, and — >) e&”*-)[w;], for s > 2; 


j=l j=utl 
and for convenience we give B;, the ith element in the wth column of A: B; is 


1 


k— 
wy! ( [ a; | + D> ( pw, yuh [ an] eh ul) + [ On: ] ens-0)) ; 
i= 


This expression for G (2, 8; p”) can now be substituted into (27), and since 
s occurs only in the last column of N, the function f(s) can be multiplied 
into this column, and the s integration performed. It is then to be noted 
that all terms of the determinant save possibly those in the wth and last 
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columns are bounded. Moreover the cofactor of F is — A, and so the deter- 
minant can be expanded to give 


I . 
I4n = sil — Ff(s)ds (a) 


taf (LEE Enereos ® a fF f(s)ds ) dp (b) 


2m, =I r=1 I=1 


toa) LE. ) ip ee — wii (c) 


=~a+l1 r=1 l=1 


is > 


It will be remembered that £ is a generic notation for functions of p (and 
other variables) uniformly bounded for large values of p. In particular, 
in (b) and (ec), E jr: is identically zero for r = 1. 

The integral (a) will be treated first. It is necessary to separate it into the 
sum of two integrals, one from a to x and the other from x to b. The result 
of integrating these by parts with respect to s gives (u/4n)f (a) + € and 
((u — 1)/4n)f (x) + € respectively, for, as is easily seen, 


f p* eP"k2-) dp, for J <b; 
AB 


f p* es) do, for I>; 
AB 


approach zero as the radius of the are AB becomes infinite, for a any real 
constant, and x distinct from a and b. One also easily shows that 


eh pe—®) d 
(32) { p 
AB p 


approaches zero as this radius becomes infinite. 

To evaluate (b) and (c) it is only necessary to notice that 1/A is bounded, 
while both D, and the s integral are certainly less than a constant times p”. 
Hence (b) and (c) reduce merely to bounded functions multiplied into integrals 
of the type (31), and so approach zero as the radius of AB becomes infinite. 

The evaluation of (d) is more complicated. We now make use of the 
more explicit asymptotic solutions of Milne, and his expressions for the g’s, (26). 
In the formula for F; in (29), the square bracket [ a; w;"t'], for example, 
has the more explicit form 


Air (8 
[ ani w yt) = -” wit ( au + > wr a’. p" =) ’ 
i (pw; ) 
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f F; f (s)ds. 


It can be broken up into the sum of several integrals, since F; isa sum. Of 
these integrals consider first the one with the exponential involving a; (/ + 1, 
k). When the more explicit form of the bracket symbol given above is 


Consider now 


used this becomes 


a, 
pin ( [ 7 Qi wit! epwdar—) F (8 ) de 
v/a 


j=! 


b n 
+ - 2 acy; Wit CPA" F( 9) ds (a) 


a, j=uti 


+ fod: + ir (s) epuar—8) wy" f(s) $) ds 


jai rl 2 (ps) 


+f - + er (s) errsar-s) wi)" f(s) ds (b) 


j= de r=1 2 Cow,)" 


+o b 2 senses wt? f(s) ds 
a 


j=1P 


oh 
- . Gs _— d epufer wi" f(: $ ) ds ) : (c) 


a j=utiP 


Now integrate (a) by parts and add the results, obtaining 


2V (a1) 0 (= wh) + Dw fay enero 4 usaf (bene® 


j=! j=1 j=util 


ay 
+ > wysernsar—) f’ (s ) ds + Me, wi ehti(ar-*) - ( 8 ) ds). 


a@ j=l a j=purl 


Since n > hy; > 0, the coefficient of f(a,) is zero. The remaining integrals 
can again be integrated by parts,* and again the terms involving a, will drop 
out. This is repeated h;; — h; times, which will give the factor a;;/p’~” 
outside, and this is just enough to cancel the factor in (33). Then absolute 
values are taken, f“*-*—-” (s) is replaced by the upper limit to its absolute 
value, and the exponentials are again integrated, giving E/p. The total 
result may be written 

(34) bi Ej; i ePrs(ar—2) + > Ej; put ePy(4r—>) + Ei 


j=zil j=util p 


Next consider the rth term of (b). The factor 1/p’ can be removed from 


*It is to be noted that the number of derivatives here required has been provided for in 
the hypothesis of the expansion theorem, and also that if v < 1 all this work with Milne’s 
formulas is unnecessary. 
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under the sign of integration, and a:;,(s)f(s) can be considered as one 
function. The integral can then be treated as was (a), but it requires only 
hy, — h; —r integrations by parts and one further integration. The result 
can then be written in the same form as (34), though in this case, to be sure, 
the E’s are not only bounded but actually approach zero. The terms (c) are 
obviously of the form E/p, so that the whole of (33) can be written in the 
form (34). Finally a single integration by parts reduces the terms involving 
the points a and 6 to the form E/p, so that the integral 


f F;f(s)ds 
k—1 


nm n E; 
35) E(B remene 4 SF Byypert emer) 4B, 
j=! 


i=e \j= j=utl p 


can be written as 


When this is substituted in (30d) the last term here gives a result of the form 
(32). The other terms for 7 + u give results of the type (31) with the value 
a; for x , and since for these terms a; is distinct from a and b, the result of the p 
integration approaches zero as the radius of AB becomes infinite. The only 
terms left from (35) then are 


k—1 
> E, ph ereuu-) 
l=2 


which are multiplied by (E/A)e?"«*-® for the p integration. If it can be 
proved that 


hy—hj—-1 ePmy(41—42) 


t=! Fp 


is of the form E/p, then this last term will give a result of the form (32), 
and (30d) will have been shown to approach zero as the radius of AB becomes 
infinite. Now inequalities (16) and (21) prove this for »y > 0, while for 
v = 0 the power of p is negative while the exponential is bounded and A re- 
mains uniformly away from zero, and the same result is obtained. 

To sum up the results, for the radius of AB large 


Inn = if (2) +. 
The work for the case R(pw,) > 0 runs along parallel to this in every 
detail, and the result is obtained that for the radius of BC large 
Inc = if (2) +€. 


This result holds for both S; and S;, so that (1) converges to f(x) in the 
odd order case. 
Even order case. In the even order case the paths of integration are the 
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same as in the odd order case, but now there is no need of dividing farther 
the ares 7: and yz, since throughout each region S the sign of the real part 
of pw, is unaltered. The region S, in which R( pw, ) = 0 will be considered 
first, that is the integral 


1 b 
I, ~ sa f np" 1G (2,8; p")f(s)dsdp. 
Y, ° a 


In making use of the expression (3) for G(2, s; p") we now alter the last 
column by adding to it each of the columns 1 to w multiplied by 39; (s) 
with corresponding index, and each of the columns » + 1 to n multiplied 
by — 39:(s) with corresponding index. The asymptotic solutions are then 
substituted, and the factor p™ of D is divided into the (1 + 7)th row, and 
the factor e’”-® into the jth column. The factor e?’«+1©°-® jis also divided 
out of the (uw + 1)th column and out of D, thus leaving in the denominator 


A = ePrulb—a) A, - 


since W,41 = —w,. Thus G(2, 8; p") is found to be equal to 1/(np"" A) 
multiplied into a determinant of order n + 1 in which the element in the 
first row and jth column is e?°@—® [1] for 7 =m and e”@- [1] for w+ 1 
=jz=n. The last element is F, given in (29). In the (1 +7)th row and 
jth column the element is [a1; w;‘] for 7 <u, [ax: wj'] for pw +1<jen, 
while for 7 =, » +1, and n+ 1 the elements are respectively D;, | 
and F;, where D; is as given in (21a), while H; is the H; given there, multi- - 
plied by e°"#-® , and F and F; are as given in (29). 

The variable s appears only in the last column of this determinant, whence 
f(s) can be multiplied into this last column and the s integration performed. 
Since all terms, save possibly the D’s, H’s, and F’s, are bounded, the expansion 
of the determinant takes the form 


1 -— 
Land J — Ff(s)dsdp (a) 
D; Hy f Fif(s)as' 
ual 
D, Hy SF,f(s)de (Lorre 
| t=1 
D, H, S P-i(s)ds| 


1 n 
2m 


ao 


% i,j,r=1 





+D E eon) dp (b) 


f=p+1 


, oP w(2—a) [ 1 ] EP um +1e—>) [ 1 ] 0 | 
H, J Fif(s)ds| Bs gy, (c) 


H, S Ff (s)ds 
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By comparing (a) with the similar integral in the odd order case it is obvious 
that it approaches the limit f(a2)/2 as the radius of y:; becomes infinite. 
In (b) the determinants can be written in the form Ep*, where a < 3n. 
Since in the parenthesis the exponents have negative real parts on y:, (b) falls 
under the type (31), and so approaches zero as the radius of +; becomes 
infinite. Finally (c) is expanded in terms of the top row of the determinant, 
and the result is treated in a manner similar to that used in the odd order 
case. The expanded form is 
«Sf Fi fds | 1|Di Jf Fifds 
(36) FE, e"«#?-® pi! : + Kh, emasie-o FI 
jd H, SF; i fds | ad \D, JF; fds | 
The formula for SF: f (s) ds is given in (35), but a more explicit form for the E 
terms in that expression is now needed. From (28) can be obtained directly 


a formula for S Fif (s)ds. It is 


p" oof F; f(s)ds = 5 > W (ys) f 9;(s)f(s)ds 


j=urtl 


bil 


37), [- D www rf a(ors(sdas 


+ fons +2 Wu (ys) f 9;(s)f(s)ds. 
— Lw ui ( (ys) f O5(s)$(4)de - ™ 


By integrating by parts vy times, 


T p—put—a) cor 
f e~P"l— F (9) ds =“ . uy U(r) +" = fe) +2 Fare 
where 
(s 7) rw) (7) 
f(r) f(r) +E + fo, 


and [f(o)] has the same form. And now this can be applied to each term 
of 9; (s) in (26) with the result 


np" o e the oe Pwkt—a) e7 Puto—a) E 
mf" as(s)f(s)ds = h(n) + ho) + 
where the symbol [h(7)] has the more explicit form 
hi(r) , he(r) hy_, (7) 
st ali dal pw; (pw;)?*  —— * (pw;)-"” 


and [h(o)] is the same with the value o substituted in place of tr. The 
function h is the function f, but a different notation is introduced because 
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the functions h, that enter are not precisely the derivatives of f. The essential 
thing in the expression is the fact that the functions h, are independent of j. 
When this result is used in (37) and the asymptotic solutions are substituted, 
the result must be of the form (35). Then by comparing the two forms the 
more explicit form for the £ terms in (35) is found: 


k-1 “ 
[ Ff (oyds = > (Le [anu] | [h(a) | p™ hi] pp; (a,—a) 
a l=: “es 


(38) " E 

+ 3 whan] [h(b)] prt corer) 4, 

j=urti p 

in which [ a; ] is identical with the [ a; ] appearing in D; (i. e., when the bracket 
symbol is written out in full for any value of m, all corresponding terms as 
well as the E terms are identical). In the same way [a] is identical with 
the [a);] appearing in H;. This evaluation of the integrals is now substituted 
in (36), and each determinant there can be written as the sum of several, 
according to the terms of (38). Since from (25) the expressions D,/A and 
(H;/A) are bounded on 71, the determinants which involve the last column 
E/p are themselves of the form E/p, and that part of (36) is of the form (32), 
and so converges to zero as the radius of 7; becomes infinite. On the other 
hand the terms involving exponentials other than w, or w,,, reduce to the 
form (31), with one of the constants a; (1 +1,k) in place of x. Hence 
these also approach zero. All that is left is then 


iT. 1 H; N; ( o> | \D; N; 
er w (x—a paces + E, ee uti t— aie , 
ci A\H,; N; cIAID; N; 
where 

k-1 
Ni = Do wh Lous] (h(a) ] pt eres 


l 


Il 
to 


$5 ws [ae] [A () | pA erate 
l=3 
By inequality (25) the expression (N;/A) is of the form E/p. So when the 
values of D; and H; are substituted from (21a) the first two terms of each 
give expressions that reduce to the form (32) and so approach zero. What 
still remains is 
E, e?*™ue-® Ayi Bai + Boi} By ee utie@-)| Aci Bai + Bai 


(3) ————_ ¥ | ++ ——-—_ :, 
A i,j Ajj Baj + Bo;) A Ag; Ba; + Bo; 


where 


Agi = ma [ ay; | e?"#—™ ( pw, )"—™ wi 


2 
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[us] er" ( — pug) *(— wy), 


Bai ; [ax] [h (a) ] tiie oP pu-®) wr, 


Bs; = [ a; | [ h (b) ] p-a-1 eP™ w(o—av) ( — W, )**, 


These can be expanded into the sum of a large number of determinants. 
Such determinants as have the same exponential factor in both columns 
are zero, for example 


[an | eP™u(o—a) (- pw, -™ ( — wu, yh wr [an] [h(b)] pi! Pu (o—a2) 
[aij] "9 ( — prog) ( — wy) wh [airs] Lh (b)] ples ertate—ee 
for when the common factors 


P 1 
eu) [h(b) ](— w,)™ (w,) pre aie” 


(it is to be noted that [h(b)] is independent of the subscript) are removed, 
we have 

lau] [ar] 

[aij] [ar] 
which is obviously zero, since the two terms in the same row of the deter- 
minant are identical. 

Since a2 — a and b — a,_, are the shortest of the intervals ending at a or b, 

the terms that approach zero most slowly are e?"#“?— and e?"#-*) | and 
each of the other determinants in (39) can be written 


Ep? er ulee—® eP—e-1) | 


This divided by A can, by inequality (25), be written E/p, whence (39) is 
of the form (32), and so approaches zero as the radius of yy; becomes infinite. 

To sum up, 

I, - 2f (x) + € 
for the radius of 7; sufficiently large. 

The treatment of the integral in S., that is over 72 is, step for step, parallel 
to the above, and the result is obtained that for the radius of 2 sufficiently 
large 

Ty, - af (x) + €. 
By combining the two results, the convergence of (1) is established in the even 
order case. 





A MEMOIR ON THE DOCTRINE OF ASSOCIATED FORMS* 
BY 
0. E. GLENN 


The theory of associated forms, founded by Hermite,t and developed by 
Clebsch, and Sylvester, and by many others, has as its primary object the 
study of sets of concomitants of a system of forms F;, F2, --- , such that every 
concomitant of the system is quasi-reducible, on multiplication by a power 
of a certain concomitant of a set, in terms of the members of this set. By 
expansion of each form of a system of binary forms, in terms of a pair of linear 
covariants as arguments, there is obtained the so-called typical representation 
of the system in which all coefficients are invariants; and every concomitant 
of the system is rational in the associated set consisting of the totality of these 
invariants, the two argument forms, and a multiplier which is the resultant 
of the two arguments. Extension to typical representation of binary forms 
of even order in terms of three quadratic covariants as arguments has been 
made by Clebsch.t 

In the first section of the present paper is given a generalization of the above 
theory, in which the system of ground forms is represented typically in terms 
of two quadratic covariants, and, under restrictions, in terms of two nics. 
This generalization is both natural and comprehensive, and gives remarkable 
syzygetic relations among the forms of a complete system. 

The second section is devoted to associated and to complete systems of 
concomitants of a binary form f under a general linear transformation f¢, 
which are functions of the coefficients of t as well as of the coefficients and 
variables of f. These concomitants have various unique properties; for 
instance, the factor r of a concomitant C in the equation expressing the in- 
variant property, C’ = rC, is not, in all cases, a power of the determinant of ¢. 

In the third section one of a system of ground forms is assumed to be a 
universal covariant of a special group. If the usual domain of rationality of 
concomitants is then altered by adjunction of the multiplier invariant of a 
certain associated set (an absolute constant), this set furnishes a fundamental 
system under the special group. Certain obscurities in the theory of in- 


* Presented to the Society, under different titles, October 28, 1916, and December 27, 1916. 
j Hermite, Journal fiir Mathematik, vol. 52 (1856), p. 1. 
t Theorie der bindren-algebraischen Formen (1872), p. 410. 
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variants and covariants of special transformations, which have persisted 
since the time of Boole, are made clear by thus identifying their theory with 
the general doctrine of associated forms. 


I. GENERALIZATIONS IN THE THEORY OF ASSOCIATED FORMS 

1. Formal groundwork of the theory. I showed a few years since that the 
binary form f of order n(m +1) — 1, whose coefficients are arbitrary vari- 
ables, has a unique expansion* in terms of any two binary nics fin, fon, 

fin = gait nayattast:++, fon = Boat + nBiay'aet+-::, 
provided only that the resultant R of fin, fon does not vanish, in the form 
m 

: m — ee 

(1) 6=2 (7) dee Saws 
° i=0 
The coefficients ¢;,-1 (4 = 0, ---, m) were determined as binary quantics 
of order n — 1 in 21, x2, with coefficients linear in the coefficients ag, ---, 
Gn(m+1)—1 Of f and rational in the quantities ao, «-+, @, Bo, «++, Bn; integral 
also in the latter, save for a common denominator Rk”. These conclusions 
conformed to the theory of a symbolical basis of a finite expansion, according 
to which f is represented symbolically as the mth power 
— — _ nt 1) —sm 
f = (Zi fin + =2 Jon =~, = =“, » 

and the quantics @in-1 = Zi‘ =; (i = 0, +++, m) are all of the same type, 
in the sense of covariant types; satisfying the equations 


a wv 
ée (a 55) din = (m — 1) di+in-1, 


0 , ) ) — 
- (s 5) bins = Ij—-1n-1; ( B 5.) = Bo rd +-:- +46, yy 
The infinitude of symbolical expressions 
F = Dik ( BE’)? (EE")¢ +++ By By ee, 
constructed by analogy with the symbolical concomitants of a binary mic, 


represent rational and integral isobaric functions of the forms ¢jn-1, fin, fon, 
and all F’s are annihilated by the three isomorphic operators 


0 0 
Fo — fin 5» 


i] 


(a) 


(6) mbna s— + (m — 1) ee Oe ee 
) MOin— ;_ on— incline mn— .. n ’ 
is Opon—1 P 2 ’ Obin-1 ; Obm—1n—1 , Of on 


*These Transactions, vol. 15 (1914), p. 72. 
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( 0 4 0 4 ‘ 0 ) 

Cc) —_ a=, as: we'd An > ’ 

. Bo dB; dB, 

and therefore by [8(0/da)]also.* This is a sufficient as well as a necessary 
condition for the invariance of the F’s under the formal alterations of the 
quantics fin, fen, given by the equalities 


fin - Ai fin + M1 fon ’ Sen nas Ae fin + Me fon ’ 


where (Au) + 0; fin = ao at + nay ay tae + +: 
In this paper much use is made of the explicit expressions for the forms 
¢in—1.in terms of the a’s and f’s. If we assume 


— n—9? ~— ‘ 
Pin-1 = Pio Xi . + Pil L, ~ Xe S i we Pin-1 Y2 l (1 =0,°-++,m), 


and regard equation (1) as an identity, it gives just a sufficient number of 
linear equations in the quantities p;; for their unique determination by 
Cramer’s rule. The absolute terms of this linear system are respectively the 
coefficients a9, 41, *** , Gn~m+1y—1 Of f. The matrix of the system is the matrix 
of the coefficients of the n(m-+ 1) forms of order n(m +1) — 1 below, 
save for binomial multipliers common to elements of columns, which we omit 
(ef. (1)): 
x} ah fin ’ x} x —_ fon, dais’ zi wt fe 
(j=0,---,n—ljt=n—-j-—1). 

The determinant of this matrix is a constant times #m™"t) + 0. 

Example. If n = 2 and f is any form of odd order 2m + 1, the matrix is 


Poo Po. << Pa "+> Tet 
ap 0 -+> att Bi 0 -+- O ao 
| m m—i i 
bay m =a bai Bo m—t Qi Rave m a 
1 ay | 1 ay * Bo 0 1 


(2) . . . . . . . . . . . . . . . ° . . . 
3" am ge! an a at* B; gn" ax Bi gue! Br 
fo a tae a: 
mn mn — 1 mn ‘mn — | mn — 1) "2" 
s™ an 7" “a B s™* Be 
0 — sss O (oe ae Aom+1 
mn mn mn 


where the column on the right is the column of absolute terms, and 


a 3 F a 
= 2 3 2 2 Gan 
= 2a aT + ot 5g, + Prag, + P56, 





° White, Am erican Journal of Mathematics, vol. 17 (1895), p. 235. 
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If we operate upon the determinant D of this matrix by either one of the oper- 


ators 
(12). (+8 
B5a F a5;), 


the result is a sum of 2m + 2 determinants each having two columns alike 
and therefore zero. For, evidently, 


0 i pi ° m—i— i+ 
(s5.) 12 don = (m — 2) fi, "Sid's 


d pPm—i fi spm—i ri— . 
(«3,) 12 foo = Ute at hy (i =(Q,--+,m). 


Hence D is a combinant of fis, foo, and, as the fundamental system of com- 
binants of two quadratics contains only one invariant, their resultant R, 
we have D = kRi"™*») , i: being numerical. 

For similar reasons the determinants of the augmented matrix corresponding 
to the solution by Cramer's rule for the respective quantities 

P00, Pio, ***s Pros ***» Pmo;s 

are all derivable from that corresponding to poo by successive operation by 
— [a(0/08)] upon the latter. Also, successive operation by — [a(0/08) ] 


upon the determinant of the augmented matrix ccrresponding to the solution 
for po. gives the determinants of the augmented matrix corresponding to 


Pols Pilly *** 5» Pity *** sy Pmi- 

A similar type of expansion of a binary form f,, of order m, with arbitrary 
coefficients, of which we shall make use, is one in terms of argument forms of 
unlike orders, fin, , fon,. We prove the following 

Lemma. If fin,, for, are any two given binary quantics of respective orders 
n,, No, and if 


m—ng m—ng—l m—ns ° 
Pim—n, = Pio xy ” + Par ™ v2 + dls + Pim—n; x2 ” (4 


then a necessary and sufficient condition in order that the linear system in the 
quantities p;; obtained by making 


(3) Je = a + em—n, Sons (nm + ne — m+ 1 ‘ 


identically, should be consistent, is that the resultant R of fin, and fon, should 
be different from zero. 
To prove this we assume 


a -1 oll 
fin, = Go 2 + a 2] te + °° fang = Bo xi? + Bi zi 2+ -°, 
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and, as R +0 precludes ap = 0, Bo = 0 simultaneously, we let ap + 0. 
The matrix of the linear system is now 


Pio Pu Pim—n, P20 P21 


P2m—nz 


ao 0 see 0 +» 0 ao 


Q; ao -e- 60 | ses 0 ay 


0 coe Qn, cee || Om 


If we delete the first ¢ = m + 1 — mn; — ne columns of 6’s from M and apply 
to the resulting square matrix M’ properly chosen elementary transformations 
we get a square matrix in which the first ¢ elements in the principal diagonal 
are ao and all other elements in the first o rows and first o columns are zero. 
Then the complementary minor of the minor of order ¢ in the upper left-hand 
corner of M’ is the dialytic form of the resultant R, and, if D is the deter- 
minant of M’, D = ag R. Hence if R + 0 the rank of M is m + 1, and the 
lemma is proved.* 

If ny + nm. = m+ 1, theno = 0, and the expansion (3) is unique. 

Again if the order m — n; of dim—n, (¢ = 1, 2) is sufficiently large, in fact 
if m—nj +12,+ m2, the form dim», may itself be developed in an 
expansion (3) in terms of the argument forms fin, , fon, 


, , 
PDim—n; - Pie~ag~m, Jim, + Pan—ay—n, Sans 


and, so far as the orders permit repetitions, v repetitions of this process for an 
arbitrary order m gives an expansion more explicit than (3), viz., 


(4) fin = DY fing + PSI fone toes FWP Sin, Soa! + WO FE es 


wherein the coefficient forms are all of order less than nj + m2 —1. This 
expansion is not unique, in general, for, as seen above, arbitrary parameters 
remain in the coefficients of the forms ¢, yY™, except in particular cases 
analogous to the case n; + mz = m + 1 in (3). 

2. Covariant expansions. Connection of the foregoing theory with invari- 
ant principles is made by regarding the argument forms in (1), (3), or (4) to 
be covariants of a definite set of linear transformations, either universal 
covariants of this set or covariants of f under the set. The following existence 
theorem proved for expansions of type (3) holds as well for type (1) and the 
unique cases of (4). 

_ Bocher, Introduction to Higher Algebra (1915), p. 46. 
Trans. Am. Math. Soc. 29 
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THEeoreM. If fin,, fen, are covariants under a definite transformation or 
group G, and if ny + nz = m + 1 80 that the expansion is unique, then the forms 
dim—n, » P2m—n, are covariants of f under G. 

For proof let (3) be written in the form 


(5) Ju = o1( Ao, Ai, +**5 21, X2) fin, + o2(Ao, Ai, -°°5 71, X2) fons, 

where A; is the rational function of ap, a1, --+ obtained by solving (3) for 

pij (j = 0, «++, m — m) and A; is the similar function equal to po; (7 = 0, 
++, m— M2). Let fm be transformed by the general transformation of G, 


viz. 
° , , 
M1 = AT) + 2X2; 


Zo = Bor + Biz (a1 Bi — a2 Bo +0), 
and suppose f., is expanded in terms of fin, (21, 22), fone (a1, %2) by formula 
(3). Then 
(6) ; _ go: (Ao, Aj i eS ay ’ 2) fin, + 2 (Ao, Ai, Ls 5 2 , 2) fons» 


where A;, Aj are the same functions of a), a}, «++ that A;, A, respectively 
are of aj, a1, :-:. Now if the inverse transformations on the 2’s and the a’s 
are applied to all forms in equation (6) there results 


(7) fm = ¥i( Bo, Bi, +++5 t1, 2) ei fin, + W2 (Bo, Bi, ++ +3.t1, 22) C2 fon, 


in which ¢;, ¢2 are functions of the coefficients of the transformations ¢t only, 
such that 

fi (2x; ’ 22) = efi (x1, Xe) 
Moreover ¥; = ¢1( Ao, ---), Ws = G2 (Ao, «++ ) on account of the equalities 
of t. This expansion (7) must be identical with (3) since (3) is unique; and 
hence 


(8) 


’ , ’ , = 

oi (Ay, Ao, o9*5 21, 22) C| o1( Ao, Ai, “+3 a1, %e), 
, , , , oe 

2 (Ag, Ai, ***3 1, 22) cz g2(Ao, Ai, **°3 %1, X2), 


which proves the theorem. 
For expansion (1) the equations analogous to (8) are 


(9) Pin-1 (Ajo, Ain, rag x ’ a2 ) = cy” cz ' Pin-1 (Aino, Aa, eee5 Sis x2) 


(i =0,-+++,m). 


In the case ny + ne < m+ 1 in (3) and for the non-unique cases of (4) 
the above proof does not hold. In numerous particular instances, however 
(ef. § III (26)), I have found that the arbitrary parameters which remain in 
the coefficient forms ¢ in these expansions can be determined so that relations 
analogous to 


vi = o1(Ao,°::); Yo = g2(Ao, *-), 
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hold. The functions ¢1 = ¢im—n,, $2 = G2m—n, are then covariants whose 
coefficients involve arbitrary parameters. 

3. Generalized typical representation. Let S be any system of binary forms 
subject to general algebraic transformations t, and let K be a complete system 
of concomitants of S, which contains two covariants fin, fon of order n, whose 
resultant does not vanish.* Then every covariant C of S whose order is of 
the form n(m + 1) — 1 can be developed in a covariant expansion (1) and 
is thereby represented typically in terms of fin, fon as arguments, and of other 
concomitants of K of orders <n. The multiplier invariant in this repre- 
sentation is the resultant R of fin, fon. 

If n = 2, C represents any covariant of S of odd order. Thus the con- 
comitants in K of orders = 2 form an associated system for the representation 
of all covariants of odd order, in which the multiplier invariant is the re- 
sultant of any pair of quadratic covariants whose resultant does not vanish. 

If we multiply any covariant C of even order of S by a definite linear co- 
variant / in K and expand IC in the arguments f\2, fos by formula (1) we get 
an expression of C in terms of the covariants fi2 , fos and concomitants of orders 
< 2 in which the denominator concomitant is of the form /R*. 

Again, covariants of S can be represented typically in terms of two co- 
variants of respective unlike orders m;, m2, when the orders of these covariants 
satisfy conditions which make expansion (3) or, more generally, expansion (4), 
unique. The multiplier invariant is then the resultant of the argument 
forms fin,» fons- 

The simplest system S affording numerous examples of this theory is the 
well-known system of a quadratic f and a cubic g. The fundamental set K 
here contains the following fifteen forms: Invariants D = (f, f)?, (f,4)?, R, 
(f*, g?)®, (f%, gQ)®, R being the discriminant of g; linear covariants (f, g)’, 
(f?,9)®, (f,Q)*, (f7, Q); quadratic covariants f, A = (g, g)*, (f,A)'; 
cubic covariants g, (f,g)', 9 = (g,A)!. 

The problem of determining the associated system A for orders = 2 in this 
case requires only the selection of two of the three quadratic covariants as 
argument forms and expressing their resultant, the multiplier for the asso- 
ciated system, in terms of the invariants in K. All covariants in K of orders 
greater than 2 are then expressed in terms of A by expansion formula (1) (or 
(3) with ny = ne = 2). 

The resultant p of any two quadratic forms is readily determined as the 
discriminant of their jacobian. Thus, writing J = (A, f)!, 


3p(A,f) = ((4,f)', J)? = ((J, A), f)?. 


* For a proof that two quadratic covariants with non-vanishing resultant exist for every 
form of even order > 4 cf. Clebsch, loc. cit., p. 410. 
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Applying the Gordan series 


we have 
(J, A)! = 2hf - 34 (f, A)?. 


3p(A,f) = 2RD = 21 (f,4)?f. 


2p ((f,A)',f) - ZRD* = iD[(f, A), 
2e((f,4)',A) = {RD — {R[(f, A). 


Hence, 


Similarly, 


We append the resulting representation of the covariant g when the argu- 
ments are chosen to be A, f: 


2L(f*,Q) +2(f,9)(F, APIA — aR, gf 


— {RD —[(f, A)?P} 


(10) g = 


In the case where S consists of a single binary quintic f the fundamental 
system K consists of 23 concomitants, of which three are quadratics, viz., 


t=(f,f), (*#,H7), (#, 7), 
where // is the hessian of f. Selecting the first two of these as argument 
forms, all covariants in K of odd order > 2, of which there are eight, are 
represented in terms of an associated set of orders = 2. The hermitian skew 
invariant does not occur in the multiplier invariant p as its degree is 18, whereas 
the degree of p is 16. 
II. SysTEMS BELONGING TO PRESCRIBED DOMAINS 

1. Domains of rationality. The general transformation ¢ leaves fixed a 

point (ka, kb) provided that 
8 = (a, 8 + a2)/(Bos + Bi) (s =a/b). 

Thus ¢ leaves invariant the two linear forms 

fi = 2Boa1 + (81 — a1 +A) 22, fur = 2Boa1 + (81 — a1 — A)ae, 
where A = ¥(8; — a1)? + 480 a2, the roots of fi, f-1 being the poles of t. 
Also t has the covariant f; f_1, or 

J = Bort + (Bi — a1) ate — a2 23. 


Ordinary rational integral concomitants of a form under ¢ belong to a 
domain, which we designate by R;(1, 0,0), determined by the coefficients 
and variables of that form. The covariant.J belongs to a domain obtained 
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from R; by adjunction of the coefficients of the transformation t, which will 
be indicated by R2(1,¢,0). Moreover, f; and f_; belong to a domain 
R, (1, ¢, 4) obtained by adjunction of A to Re. 

The forms f;, f_1, J, and the concomitants considered below in this section, 
being functions of the coefficients of ¢, are primarily invariants of a single 
transformation ¢, but a; and 6; appear in them only in the combination 
a, — 1; hence, if € is a fixed constant, our invariants belong to a one-param- 
eter continuous group leaving two points fixed, and characterized by the 
equation a; — 6; = 2e8). An important sub-problem under the general 
case is given by assuming 


a, — fi = 2€Bo, 3 = — Bo. 


Then J, deprived of a constant factor, is a universal covariant of the set 
given by 

% = (2€Bo + Bi) 21 — Bo a2, 
* "ae = Boi + Bias, 


of which the orthogonal transformations are a particular case. All of our 
systems, particularized for ¢t;, belong essentially to the domain R;(1, 0,0). 

Both orthogonal and boolean systems are special cases of the general 
systems given in this section. We shall revert to this phase of the question 
in Section III. 

2. Complete systems in R, (1, ¢, A) of the transformation f. If a binary 
form f = ao 2{' + ma; 2 22 + --- be represented typically in terms of f,, 
f-1 as argument forms, the multiplier invariant, a power of which occurs in 
the denominator of the expansion of f, is 4894. This is an absolute constant 
and is, essentially, the quantity adjoined to define R,(1,%,4). Writing the 
expansion in the form 


(11) f = dm r + Mom—2 - = + inti + mo—m—2) fi fr" + o-mfri; 

we have in ¢m—oi (1 = 0, ---, m) a set of m+ 1 linearly independent in- 
variants of ¢ belonging to R,, linear in the coefficients ao, --+, Gm. The 
factor r of @m—2; in the equation expressing the invariant property 


, 
Pm—2i = Thm—2i » 


is not a power of the determinant of ¢, (a1 61 — a2 Bo), but is a product of 
powers of the factors pi, p_: of this determinant in R,, 


(12) p= 3(a+6i1+A), p-1=3(a+8,-—A). 


t 


In fact, writing c: = p;', ¢2 = p—}, it results from I (9) that 
(13) bn—2i = Ci” C2' m2: 


= pi ** (a1 Bi — a2 Bo)* bm—2i (¢ =0,--+,m). 
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Moreover there foliows the 

THEeorREM. The m+ 3 forms fi, f-1, dm-2i (0 = 0, -++°' , m) of the 
associated set constitute a fundamental system of concomitants of f under t in 
Ri(1,t, A). 

The explicit systems in R, for m = 1, 2, 3 are given below, in which, to 
secure isobarism, we write 8; — a, = Y1.- 


TABLE 
m=1 
2? Bo Adi = (— ¥1 +A) ao + 28a, 
— 2? Bo Agd-1 = (— y1 — 4) Go + 2B a1. 
m=2 
(¥? + 2a2 Bo — y1 4) ao + 280( — v1 +A) ai + 268 a2, 
2Bo A? do = a2 Qo + ¥1 41 — Boa, 
23 B2 A? d_2 = (7? + 2a2 Bo + ¥14) ao + 280( — v1 — A) ai + 285 ae. 
m=3 
— 2¢ 83 A® bs = (vi + 38071 a2 — (77 + Boar) A) ao 
+ 3Bo( — vy? — 28002 + 71 4) a; + 383 (v1 — 4) a2 — 285 a3, 
— 2! B3 A® di: = Boar (71 — A) ao + Bo (vi — 2Boa2 — 714) a1 
+ B ( — 37: + A) az + 265 as, 
24 B3 A® g_1 = Boar (v1 +A) ao + Bo( yi — 2Boae + 714) a 
+ B32 ( — 371 — A) a2 + 265 as, 


23 82 A? o2 


ll 


2! BS A? o_s = (7? + 38071 a2 + (7? + Boar) A) ao 
+ 380 ( — v? — 2Boa2 — ¥1 4) a; + 363 (v1 + A) a2 — 285 as. 


3. Systems belonging to R.(1, ¢, 0). It will be shown hereafter that 
complete systems of f in R2 may be constructed from associated systems by 
taking the set S of ground forms to be the quantics f, J. Since, however, 
the forms ¢m~2; do not satisfy some of the hypotheses (cf. (13)) fulfilled by 
ordinary invariants of f under the group of general linear algebraic trans- 
formations, it is desirable to develop a special method for systems of f under ¢ 
in R,, and later to identify the results obtained by the two methods to be 
employed. . 

Concomitants in R; or Rz may evidently be constructed by forming linear 
expressions y in products of powers of fi, f-1, dm-2i (¢ = 0,---,m).* A 
necessary, though not sufficient, condition in order that such an expression 
may be an invariant function in R, is that the exponent v of p; (cf. (13)) in 
each and every term of the equation expressing invariancy should be 0. The 
exponent v for every term of any concomitant is equal to the sum of the 


*Cf. Elliott, The syzygetic theory of orthogonal binariants, Proceedings of the 
London Mathematical Society, vol. 33 (1901), p. 226. 
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subscripts of the ¢’s diminished by the exponent of f; and increased by that 
of f_; in any term. For pure invariants therefore, 


y= >) (m — 2). 


If d is the degree of such an invariant y and if we revert to the notation for 
the ¢’s used in I (1), dm = $00, m—2 = G10, ***, the sum w of the (first) 
subscripts in each term of y is 2i. Thus 


vy = dm — 2w, w=%3(dm-—vr). 


We have thus established a one-to-one correspondence between invariants 
in the @’s, of degree d and factor py, and isobaric invariant functions of de- 
gree d and weight w. The number of invariant expressions y with factor p} 
is therefore the number of partitions of w denoted by 


(4 (dm — v);d,m) = (3(dm+y);d,m). 


Now in the domain R2(1,¢,0) we are concerned with the case v = 0. 
The actual invariant products of the ¢’s for which v = 0, each affected by 
that power of D which is the multiplier in the invariant relation for the product, 
may be found from the caylean real generating function 


jj 


(14) G(f) =4 IL (1 = m2: pt Di) | = (D = a1 6: a2 fo). 
i=0 


The expansion of this function in ascending powers of p; has for the terms 
free from p; in it the sum of all of the products of the ¢’s which are invariants 
with vy = 0, each product affected by the power of D which is its multiplier. 
The invariants of this set which are asyzygetic are those possessing the same 
power of P as multiplier. 

Not all of these invariants belong to R.. But the products may be arranged 
in conjugate pairs (¥,, y_,), ete., and by replacing ¥,, Y_» respectively by 
Vn + Vv-», Vp — Y—p we get the complete asyzygetic sets which, neglecting 
factors involving A, belong to R,. In fact y_, is derived from y, by changing 
all factors dm—2; into @-~m—2i); and that the power of D figuring as multiplier 
in the invariant relation for y_, is the same as the power for y, is easily proved. 
For, we have 


Vp _ Pn—2%, Pin—2to ike Y—p _ Pr\m—2i;) Gr m—2iz) res 


The multiplier for y, is D?**?#*+~, and that for y_p is D™@—)tm™— ete | 
But if v = >; (m — 27) = 0, then 


pPi(m — %) + po(m —%2) + oe) = Mit tpet: 


The generating function for simultaneous invariants is I1G (f) , the product 
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including the G(f) function for each form involved, and covariants of a 
single form are simultaneous invariants of the form and a linear quantic. 
Hence the generating function for covariants of f in Rp is 


m -1 
(15) T(f) = {Ic — dm—2i pt D*) (1 — fipr') (1 — fips D*) | ; 
Thus, we have the 

THEOREM. The products of functions dm—2; occurring in the portion of the 
expansion of I'(f) free from pi, which are affected by the factor D* , furnish a 
complete asyzygetic set of concomitants in Rz of multiplier D*. These products 
have only to be arranged in conjugate pairs (Wp, ~—p) and the terms of each pair 
replaced by Wp» + ¥-p, Vp — -p respectively. 

We shall now illustrate this theorem by exhibiting the explicit systems 
given by these methods for the forms of orders 1, 2, 3. 

m= 1. The generating function for invariants of a linear form f is 


G(f) = {((1 — d1pi) (1 — dip? D)}". 


When this is expanded it is found that the terms free from p; are 
1+@¢10D+¢5¢,0P?+---=(1-—¢:¢1D)". 


Hence a linear form has only one invariant in R2, viz., @; ¢-1 of index 1. 
Moreover the terms free from p; in G(f), when separated from the rest, are 
all given by the formula* 
1 
1—¢@¢1D° 


This formula will be called the real generating function for invariants in the 

domain R,. The form of the denominator of such a function shows what 

invariants should be selected to form a complete system in this domain. 
The generating function for covariants and invariants both is 


r(f) = {(1 — dips) (1 — Gi pr' D)(1 —fipr') (1 —firpi D")}-. 


Separation of the terms free from p; from I'(f) gives the terms furnished 
by the expansion of the following real generating function for concomitants 
of the linear form in R2, 
L—G@igififin —_— 
(1-—¢:¢1D)(1—@fi)(1 — dn fa) Cl -fifirD*)’ 


Accordingly the following four forms constitute a fundamental system of f 
under ¢ in R.(1, ¢, 0) (ef. table for m = 1 in the preceding paragraph): 





* Cf, Elliott, Proceedings of the London Mathematical Society, 
vol. 33 (1901), p. 236. 
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(16:)) A=Gidu, Me*fiteafau, Ar=Gfi-—gufu, 
w=fifa, 
I = 48) A2X = a2 a2 + ¥1 G0 a1 — Bo aj, 
(16) Cy = Xu = (— 91 G0 + 280 a1) 1 + (2a2 ao + ¥1 41) 22, 
Cy = =ar1+ a2 =f, 
J = (480) w = Boat + 7121 te — a2 2. 
These forms are connected by the syzygy 


S = (48) a2 + yi) f? — C2, —41J = 0. 


where 


m = 2. The generating function for concomitants of a quadratic f is 
T'(f) = {(1 — dept) (1 — do D) (1 — 62 pr” D*) 
(1 —fipr')(1 —f.ipD")}". 
The real generating function for concomitants in R2 is found to be 


1 — $2 o-2fif-1 = 
(1 — do D)(1 — 22D?) (1 — off) (1 — 2 f21) (1 —fifuD")’ 


Hence the fundamental system of f in R: consists of the following five forms 
from which extraneous factors involving A are to be deleted: 


(171) W=¢, A=Gebo2, M=Gfitgeerfii, 
Ai=¢efi-Gefa, w=fifa. 
This system is given in explicit form below:* 
Io = 289 A? Xo = G2 do + ¥1 G1 — Boas, 
I = 1683 A*d = af a? — 4ae Bo al + Bi az + 2a2 1 a0 a1 








+ (71 + 2a2 Bo) ao a2 — 2Bo ¥1 41 a2, 
Cy = A? = (( 77 + 2ex Bo) do — 2Bo ¥1 a1 + 285 a2) ai 
(17) + 2( — a2 710 + 4e2 Bo a1 + Bo ¥1 G2) 21 Xe 
+ (2a do + 2oe yi a1 + (Yi + 2a2 Bo) a2) 22, 
Cy = Adi = (— 7140 + 2B a1) ai + 2 (a2 ao + Bo a2) 21 2 


Q + (2a2a1 + ¥1 42) a2, 
J = (480) pw = Bort + ¥121 22 — a 23. 


*In place of Cy we may use f since A? f = C, + 2JoJ. 
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A single syzygy connects these forms, viz., 
S = 41 J? — C? + (4B ae + vi)C?, = (). 


m = 3. Ina similar manner, employing the real generating function for 
concomitants of the cubic in R,, we derive as a fundamental system the fol- 
lowing thirteen quantics (cf. table in II 2): 


39-3, d19-1, bPitgdsdi, 3¢51-O3¢1, fifi, 
(18) @ifit dif, bsfitosfii, 
dorfitgsoifri1, sO ifitosdif-r. 


III. INVARIANT SYSTEMS FOR SPECIAL GROUPS 


1. Complete systems derived from associated systems. The systems 
treated in section II may be derived from certain associated systems, and by a 
method which is of importance in the invariant theory of various special sets 
of transformations. 

Construct, in terms of the forms f,, f-1 as arguments, the typical repre- 
sentation of the system S of ground forms consisting of an mic f and the 
covariant J. The representation of f is given by (11), that of J is J itself 
since f, f.1 = 48) J. The multiplier of the associated system is the absolute 
constant 485 A which we adjoin to the domain. The fundamental system of f 
under ¢ in R, then consists of an associated system of S = {f, J}, with A, 
as a constant, omitted. We can construct this associated set explicitly by 
making the following substitution in the associated system of an mic f and a 
quadratic J; = bo aj + 2b; 21 22 + be 23, in which the multiplier invariant is 
the square root of the discriminant of //;: 


(19 a; by b; bo 1 Xe ial . 
. dm—-2i 0 1 0 fi fr 7=0, »>m). 


Thus, in case m = 2, the covariant 
(ao by — ay bo) ai + (ao bz — a2 bo) a1 te + (a,b. — agb 1) 23 


becomes 


defi — d2f21, 


and this is the fourth form of (17,). Similarly the sets (16;), (17,), (18) 
may all be derived from the corresponding associated systems of f and J,. 
2. Orthogonal systems. An obscurity which has long persisted in the 
theory of orthogonal invariants may now be cleared up. To quote from the 
paper by Elliott on orthogonal binariants referred to above: “ Absolute 
orthogonal invariants of a p-ic are invariants of the p-ic and xj + x} and 
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vice versa. It has therefore been hastily supposed elsewhere that the complete 
irreducible invariant system of a p-ic and a quadratic produces exactly the 
complete (absolute) orthogonal system for the p-ic, when 2j + 2} is taken 
for the quadratic. A first case of the redundancy of the former system for 
the latter purpose is exhibited above [the reducibility of the covariant Q of 
the cubic in terms of the orthogonal system corresponding to the set (18)]. 
The search for complete absolute orthogonal systems is not identical with 
the search for invariant systems of forms one of which is a quadratic.” 

The inference here is that orthogonal systems require a special method of 
treatment not comprehended in the ordinary theory for integrally complete 
systems under the general binary group G. We can now show, however, 
that the orthogonal invariant theory is essentially identical with the theory 
of associated systems under G. For, by particularizing the transformation 
t; of II 1 into the orthogonal transformations 


, . , . , 
x1 = 2; cos 8 — a sin 8, %2 = 2; sin 6 + x, cos 8 


(making e = 0, Bo = sin 6, B; = cos 6) J becomes essentially 27 + x3, while 
the concomitants exemplified in the table in II 2 are freed from the coefficients 
of the transformation and become the non-absolute orthogonal concomitants 
of the respective forms f, the multiplier A of the associated system now be- 
coming a constant times V— 1. At the same time the sets in R., (16,), 
(17,), (18) are also freed from the coefficients of the transformation and 
become the complete absolute orthogonal systems. 

Thus the absolute orthogonal invariant system of a form f is given precisely 
by the associated system of the set f, xj + 23, where the multiplier invariant 
is the square root of the discriminant of x] + 23. 

The lowest instance where this associated set is smaller in number than the 
ordinary integrally complete simultaneous system of f and 2j + 23 is where 
m = 3. In fact, in the general case of the system (18) of the transformation ¢, 
the covariant Q of the cubic is reducible, on multiplication by a power of A, 
in terms of the set (18). 

3. Boolean systems. If we particularize ¢, of II 1 by the substitutions 





sin a 


(20) Bo = 


sin B 
== €= — cosw, 


sin w’ sin w’ 


and assume w = B — a, we get the following transformations of determinant 
unity, 


sin (w — a) , sin (w — B) , sin a@ , sin B 
21) a = 2; —— +2 : Ye = 2 H+ 2 
(21) a ’ sin w 7% sinw ” ' sin w + me sin w’ 


representing rotation of a pair of oblique coérdinate axes inclined at an angle w 
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into another pair inclined at the same angle. Then J becomes a constant 
times the universal covariant (qua transformations on 2;, x2 and not on w) 


q = 21 + 2 cos wx t2 + 23. 


Hence if we make the substitutions (20) in the systems (16,), (171), (18) 
these become the boolean systems of the forms of respective orders 1, 2, 3. 

Moreover, the invariants in the table in II 2 become, under (20), systems 
of a species of boolean invariant in which the factor in the invariant equation 
is a product of powers of p;, p_: as below (cf. (13)), 


pi = cosa+isina, p1 = cosa —isina. 
The two invariants of the linear form, for example, become 
(22 (cos w — isin w)dy — a, (cos w + 7 sin w)dy — a. 


4. Kleinian groups. The formal nature of the invariant systems treated in 
section II shows that they may be interpreted when the transformation ¢ has 
complex coefficients a;, a2, Bo, 8:. In particular, the multiply infinite set 
represented by 

ay = (2€Bo + 81) %1 — Bore, 
th: - J 
%2 = Bor + Bite, 
where ¢ is a fixed and 8», 8; arbitrary complex numbers, has the universal 
covariant 
q = xi — 2exri 22 + 23; 


and the systems tabulated in II, particularized for t;, become complete 
systems with complex coefficients, free from the arbitrary coefficients in ¢;. 
Note that the latter transformation includes all types of substitutions with 
complex coefficients excepting the parabolic type, i. e., includes the elliptic, 
hyperbolic, and the loxodromic types.* 

5. The general cyclic group on » variables. The following transformations 
generate the so-called cyclic group C on n variables: 


Ai, Ao, °°, An 
C: (x)= Ress *** 5 ieee (x). 
ho, As, °°, A 


The determinant D of C is a linearly factorable n-ary form of order n. For, 
let us multiply its columns by 1, g‘, ---, g“»*, respectively, where g is an 
arbitrary primitive nth root of unity, and ¢ = n — 1, and add all to the first 
= Klein-Fricke, Theorie der elliptischen Modulfunctionen (1890), p. 163. 
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column. The result is immediately of the form 


Re he, coe, Xn 


Ta Pe Te mee 


, 


gest, As, wee, De 
where 
Fi(A) =A. + g' Xe + g’' Xs fee tH gg? * ha; 
and hence 


D=T[ F(a). 


The following theorem is easily proved: 


THEOREM. The forms F,(x) (t = 0,---,m—1) are universal covariants 
of C,, satisfying the invariant relations 
(23) F,(2’) = Foye (A) Fi(2). 


The general quantic f of order m in n variables, whose coefficients are arbi- 
trary variables, can be given a typical representation by expanding it in terms 
of the forms F,;(z) (t = 0,--+,n—1) as arguments. The multiplier of 
the resulting associated system, i. e., the invariant, a power of which occurs 
in the denominator in the expansion of f, is the eliminant of the n linear 
arguments.* This resultant, being a function of nth roots of unity only, is a 
constant which we may adjoin to the domain of rationality. The coefficients 
in the expansion will be linear expressions in the ("*"-') coefficients of f, 
belonging to the domain R of rational functions of the nth roots of unity. 
These linear expressions, ("*"~') in number, constitute a complete system of 
invariants of f under C in R.t 

6. A transformation leaving an involution fixed. The systems in R, treated 
in section II are complete for any transformation ¢ whose only covariant in 
R,z is a power of the quadratic represented by J. There exist no covariants 
of the general transformation ¢ other than f; f_1 and its powers. This appears 
by consideration of the generating functions employed in il, but a direct proof 
of the fact is possible. Since ¢ has but two poles, any other quadratic covariant 
of t, as 

K = a? — 2ba; 22 + 23, 


must have its roots interchanged by t. But the condition 


(24) + w=; eke Sots 
‘ Bo(b + Ve? —c) + Bi 











gives Bj = — a. 
*Compare a paper by the present writer, these Transactions, vol. 15 (1914), 
p. 83. 


t For a treatment of the case n = 3, cf. Elliott, Messenger of Mathematics, 
vol. 33 (1903-4), p. 108. 








460 0. E. GLENN [October 


Under the transformation with arbitrary coefficients, 
Tit, = a2, + 02%, % = Bor — a2, 
the following quadratic forms are relative covariants: 
fie = Boti tarrz, foo = ari taza. 


These forms are left invariant by having their roots interchanged, and all 
quadratics thus left invariant belong to the involution 


L= ky fie + ke foo. 


In fact, equation (24) gives 
Boe — 2a; b — ag = 0, 


and hence the following equality is an identity: 


The covariant J = Bo ai — 2a; 2122 — a2 az is the jacobian of fiz, foe. 
Hence the 

THEOREM. A necessary and sufficient condition in order that a binary form 
be transformed into itself by T 1s that it be equal to a power of J times a product 
of powers of forms in the involution L. The double points of this involution are 
the poles of T. 

If we expand a quantic fo»+1 in the argument forms f12, fo2 according to the 
formula (1) I, the coefficient forms ¢,;; (2 = 0, ---, m) are linear covariants 
of fomy1 under 7 in R2(1, 7,0). Then it is manifest that a fundamental 
system of concomitants of fom,; under 7 in this domain is precisely the simul- 
taneous system of the m + 1 linear covariants ¢,1 (i = 0, ---, m) taken with 
the covariants of the transformation, e. g., fiz, foz. Definitive processes 
appertaining to 7 for the construction of this joint system are not known. 

We append the coefficient covariants ¢o1, ¢1: for the case where the ex- 
panded form is a cubic. If D = — aj — Bo ae, then, 


— a2 Doo = (a3 ao — 8a a2 a) + 3aj de + a1 Bo a3) a1 — Daz ae, 
(25) — a2 Dou = (a1 a2 dy + 3a Bo a1 — 3a Bo a2 — Bo a3) 21 


+ (— af ay + 3a1 a a1 + 3a2 Bo a2 — a1 Bo a3) 22. 


Note that the so-called skew-orthogonal transformations form a special 
case of 7’, but, also, that in this case the poles are real and rational, and the 
domains R,, R, coincide. 

7. Modular systems. One of the main features of the general theory of 
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sections I, II of this paper is its character as a formal algorithm under which 
a variety of invariant systems for special groups can be studied. 

Consider the group of binary transformations 7 in which the four para- 
meters are least residues modulo p, where p is a prime number. A complete 
system of universal covariants mod p of T consists of* 


L = xi %2— 212%, Q = (af ve —aat) + LD = 2fP + +) + ze, 


Regarding all equations in I 2 to be congruences modulo p, and setting 
fin, = Q, fon, = L in (3) and in (4), it appears that a binary form of order 
m = p* is reducible in terms of Q and L and two first degree formal modular 
covariants of respective orders p and p> — p—1. This raises the question 
whether the modular expansions (3) for m > p*, containing arbitrary para- 
meters in their coefficients, may always have these parameters determined 
so the coefficient forms @im—psip, Gem—p—1 are modular covariants. Assuming 
this, (4) shows that the general quantic of order m > p* — 1 is reducible 
modulo p in terms of Q, L, and its own first degree covariants of orders from 
0 to p? — 1 inclusive. 

For the case p = 2 I have answered this question in the affirmative for the 
general order m.f 

For the case p = 3 I have actually reduced in this manner the forms of 
orders m = 37, m = 10, m = 11, thus verifying the existence of the desired 
expansion in a very typical case. Thus, while the theorem given below has 
not been proved in as general a form as is stated, I venture to give it in its 
general form for the notice of mathematicians who may be interested. 

Since the proposition relates to the reducibility of the general quantic we 
may state it as a theorem on covariants. 

THEOREM. The maximum order of an irreducible covariant modulo p of any 
binary form or set of forms is p? — 1. 


For the case m = 11, p = 3, we find f = ¢15Q + dor L (mod 3), where 
$15 = ao Xi + (a; + a3 + AS) at re + (a2 + ag + ag — AT) 23 23 


(26) : 
+ (as + az + ag — AS) aj 2 + (ag + aio +AT) a1 ze + ay, 23, 


where 
S = a, + a3 + a5 + a7 + Ag, T = az + a4 + ag + Og + Ano, 


and X is any least positive residue modulo 3. The septic covariant ¢o7 is 
readily obtained as the quotient modulo 3 of f — Qdi5 divided by L. 

A complete system of formal modular concomitants of a binary form f 
could be obtained, under this theorem, as a simultaneous system of a set of 


+ Dickson, these Transactions, vol. 12 (1911), p. 75. 
7 O. E. Glenn, these Transactions, vol. 17 (1916), p. 545. 
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first degree covariants of f of orders not exceeding p? — 1. Note, however, 
that the methods by which this simultaneous system would be constructed 
would be the definitive construction processes for simultaneous invariants 
and covariants under the modular group 7’, and these definitive construction 
processes are not fully known. Nevertheless, for fundamental systems of 
pure invariants in the modular realm, Dickson’s theory of classes furnishes 
the requisite definitive methods. 
UNIVERSITY OF PENNSYLVANIA 

















THE DEGREE AND CLASS OF MULTIPLY TRANSITIVE GROUPS* 
BY 
W. A. MANNING 


If a group contains substitutions of degree u but none that displace fewer 
letters, the identity excepted, u is said to be its class.| The same term is 
unfortunately used with a different meaning in the theory of prime power 
groups. 

In substitution group theory the two problems of the order and of the 
class of groups of given degree are in large measure interchangeable.{ For 
this reason Bochert’s inequalities, limiting the degree n of a multiply transitive 
group in terms of its class u are second in importance only to Sylow’s theorem.§ 
Bochert’s theorem may be stated as follows.|| 

The class (> 3) of a substitution group of degree n exceeds 4n — 2 Vn if it 
is doubly, 4n — 1 if triply, 4n — 2 if quadruply, transitive. 

In the proof of this theorem Bochert introduces the condition that the 
class is a given number by means of the commutator of two non-commutative 
substitutions. Now the commutator is only one combination out of many 
that may be made from the substitutions of a group (@) and is by no means 
always one of the substitutions of lowest degree in G. In the hope of pene- 
trating more deeply into this matter the author has assumed the presence of 
substitutions of order 2 among the substitutions of lowest degree in G and 
thereby gains the advantages that are afforded by the simplicity of structure 
of diedral rotation groups. The principal result arrived at may be put thus: 

1. THeorEM. Let n be the degree and u ( > 2) theclass of a group that contains 
a substitution of order 2 and degree u; then if it rs 


doubly transitive, p>in—tvn—-1, 


triply - » sae 

4-ply “ , #Be(n+1), 
5-ply a) 2 w= en, 

6-ply “ > rm, 

7-ply m= ¥ p= ihn, 


* Presented to the Society, Nov. 25, 1916. 
t Jordan, Paris Comptes Rendus, vol. 72 (1871), p. 384. 
tSee these Transactions, vol. 16 (1915), p. 139. 
§ Jordan, Journal de mathématiques, ser. 5, vol. 1 (1895), p. 35. 
|| Bochert, Mathematische Annalen, vol. 40 (1892), p. 176; vol. 49 (1897), 
p. 133. 
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more than p; + po +--+ + p, times transitive, where py, po, «++, Pr are distinct 
odd primes (r > 1), 


(pi — 1) (p2 — 1) +++ (p — 1) - 1 
(pi — 1) (p2 — 1) +++ (p — 1) 





nN. 


THE DIEDRAL GROUPS 


2. The auxiliary theorem by means of which this will be proved is the 


following: 
If the order of a group of class u generated by two substitutions s and t of order 2 
and degree yw is divisible by each of the odd primes p,, po, «++, pr, its degree n 


does not exceed wp + w/(pi — 1) (pe —1) +++ (p, — 1); of the order N of st is 
odd or is divisible by 4, n does not exceed p + 2y/N; if N (> 2) is twice an 
odd number, n is not greater than p + 2u/N + 8yu/N?. 

If we erase from {s, t}, of degree n, the & letters of all the regular con- 
stituents of {s, t} which involve letters displaced by both s and ¢ the group in 
the remaining constituents is of degree n — £, of class uy — £, and is generated 
by two similar substitutions of order 2 and degree » — —. Then if the the- 
orem can be shown to be true for diedral groups satisfying the given con- 
ditions, but in which the regular constituents, if any, are of order 2 and ex- 
hibited by transpositions of s (or ¢) in letters fixed by ¢ (or s), it is a fortiori 
true when regular constituents on letters common to s and ¢ are present. 
It may be noted that if p* is a Sylow divisor of N, i. e., the highest power of 
the prime p that divides the order of st , and if p* divides the degree of a regular 
constituent of {s, ¢} it must also divide the degree of some non-regular con- 
stituent of {s,t}. If this is not so, (st)*” involves no letter not in a regular 
constituent of {s,¢} and is not the identity. Thus the degree of (st)*” 
is at least u and cannot exceed £; but if £ = uw the truth of the theorem is not 
in dispute. Consequently the deletion of the regular constituents from the 
group {s,t} leaves a group of the same order. In what follows the group 
{s, ¢t} under discussion will be assumed free from any regular constituent on 
letters common to s and ft. 

The transitive constituents of {s,¢} may be enumerated as follows. In 
s there are m, cycles that displace letters not in ¢, and ¢ has me cycles that 
displace letters not in s; let m = m, + m2. There are y; constituents of 
degree Y; and order 2Y;, Y; an odd number, j7 = 1, 2, ---. There are 2, 
transitive constituents of degree Z, and order 2Z,, Z, even, k = 1, 2, ---, 
with the generator of degree Z;, in s and the generator of degree Z;, — 2 int. 
In like manner there are z, constituents of the same order Z, with Z, — 2 
letters in s and Z, letters in ¢; let z, = z, +2;,. The numbers Yj, Yo, --- 
are all the odd divisors of N, unity excepted, and the numbers Z,, Z2, --- 
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are all the even divisors of N , with the exception of 2. We recall that a non- 
regular transitive diedral rotation group of degree n is of class n — 1 if n is 
odd, and is of class n — 2 if n is even. In the latter case the two generators 
of order 2 are not similar substitutions. Of course some of the numbers 


us w” 


Yi, Y2, ***, 21,22, °°, 21, 22, *** may be zero. The degree of s is 
2m + Dy(¥;-1)+ DxAa+ Du (Ze —- 2), 
and the degree of t is 
2m. + Do yj (¥; —1) + Doa(Z —2) + Dow Z%. 


The summations extend the subscripts through their entire range of values. 
With this understanding we may omit the subscripts j and k in future and 
introduce these abbreviations: 


E,=>dSy(Y¥-1), Fi=>doz2(Z-1), 
T,=E£e.+Fa= dD y(¥—a) +) 2(Z-a). 


and similarly 
Then the degree of {s, t} is 2m + Io which on the elimination of 2m becomes 


2u — T2; for 
— 2M.+T = dD yY — Dd (Qy¥ — 2y) + S2eZ — YS (227 — 22) 
=-Dy(¥-2)-L2(Z-2)=-k. 


3. If st is of odd prime order it is obvious that n = wp + u/(p — 1). 

Let st be of order p*, pany prime anda > 1. Let {s, t} have x; transitive 
constituents of degree p',i = 1,2, ---,a. From (st)*” we get the equa- 
tion, 

pr, =urth,, 
from (st)?*’, 
Poza t pr H=ethe, 
from (st)”**, 
Peet plait pat. =~ Ut hee, 

from (st)?, 

Pt2tpast-->+p*x, =uthe, 
while 

put pi+---+ptr,=uts, 
and 

(p—1)a1 + (p? —1)ae+--> + (p*-—1)a =p. 


Here hz, hs, --- , ha, 5, are positive integers or zero, and up + 6 is the degree 
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of {s, t}. These a + 1 equations are consistent, their eliminant must vanish, 
and when equated to zero and reduced it gives 








5 = w/p**(p — 1) — 2 hi/p*. 
Then if st is of order p* (a > 1), the degree of {s, t} does not exceed 


w+ u/p**(p—-1). 


4. Let N = 2p*, where p* is a power of an odd prime. Let y and z be the 
number of transitive constituents of degrees p* and 2p* respectively in {s, t}. 
Then from (st)*”? and (st)*” we have 


poy t+2ptz=pth, 


2m + 2p*z =pt+k— Po, 
and . 
2m + p*y + 2p*z =w+5—-To, 
with is 
m+ (p*—1l)y+ (2p* -—1)z =n -T. 


By Fy, To, and T; is meant the summation explained above but with the 
exclusion of constituents whose degree is divisible by p*. The meaning of 
h, k, and 6 is the same as before. They are positive integers or zero. The 
eliminant of m, y, and z is 


0 p* 2p* uth 
2 «0 2p" pw tk — Po 
2 p* 2p* pt+s—Te ¥ 


1 pe—1 Q2p*—-1 pl; 
This reduces to 


(p*—1)8 =p—h(pt—1) —k-— (Ey + pre), 


whence we conclude that when the order of st is twice a power of an odd prime, 
its degree does not exceed uw + w/(p* — 1). 
5. Let N = pi po «++ p,, the product of r distinct odd prime factors. 
There is at once the equation 
Ei =p, 


and there are > {=| -C; = 2” — 1 equations given by the 2” — 1 powers 
(st )*\»%~) | which is in no case the identity but is of degree » + hi;..., where 
the quantities h,;... are positive whole numbers or zero. Let 6 = hy...., a 
constantly recurring quantity since u + 6 is the degree of st or of {s, ¢}. 
Now (st)*%, say, involves letters from and only from those transitive con- 
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stituents of {s, ¢} whose degree contains p; or p; as a factor, and so on for 
other combinations of the primes pi, p2,---. The equation given by (st )*”'”* 
when r = 3, for example, is 


Pr Yi + Po Yo + Pr Po Yi2 + Pi Ps Ys + Po Ps Yos + Pi P2 P3 Yrs = w+ hy. 


There are 2” — 1 variables yi, ---, yi, «++ in the 2” equations. Replace 
these equations by an equivalent system got by replacing the first, E; = yp, 
by Ey — E; = éand the last 2” — 2 by themselves less the second, Ey = uw + 6. 
The eliminant of these equations is 
1 11 +: 1 1 5 | 
| 
Pi P2 Ps *** P2Ps-** Pr Pipros** Pr MS | 


QO po ps +++ PopPzs+** Pr Pi po+** Pr -~h+8|7° 


| 
| 
cf 
Multiply the first row of this determinant by p; pe --- p, and then divide the 


first column by 71, the second by po, ---, the (r + 1)th by p; po, ---. Then 
we have 


P2 P3*** Dr P1P3°** Pr eee P2 Pi 1 Pi P2 *** Pr b| 


1 1 me Ee EY wuts 
0 1 + 0 10 -h+8 |_y 


1 0 ao § 8 4 — he +8 | 


For the sake of convenience of reference we may call the row in which h;, 
hi;, «++ occur the h;-row, the h,;-row, and so forth. Then the law of the 
formation of the last 2” — 2 rows is evident: a term in the h,-row (except 
5 — h;) is unity if the term above it in the first row contains p; as a factor; a 
term in the h,;-row is unity if the term in the first row above it contains p; p; 
as a factor, and so on; otherwise it is zero. For wherever a zero is found in 
the h;;.-row of the determinant there was a term 7p; p; px Yi; in the original 
equation, from which was subtracted Ey = up + 6, and this product p; p; px 
we divided out from every term of the column corresponding to y;;,, thus 
leaving p1 p2-++ pr/(pi pj px) in the first row directly above the zero in 
question. The present units correspond to zeros in the original equations. 
For the sake of clearness it may be well to write out the determinant for the 
case r = 3: 
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| P3P2 PsP1 P2Pr Ps P2 Pr 1 Pi P2 Ps 6 | 








1 1 1114141 pté 
0 1 1 0010 6-hk 
1 0 104100 8-hkh 
1 1 0 10 00 b&—hs _ 
0 0 1 00 0 0 6—hy 
0 1 0 0 0 0 0 b6—hs 
1 0 0 0 0 0 0 &— hes 


Now multiply the second row by — 1, the h,-row by — (p; — 1), the hj;-row 
by — (pi — 1)(p; — 1), ---, and add them all to the first row. The sum 
in a typical column, the y2...,-column, is 

k i=k; j=k 
Pi P2*** De -|1 +20 (pi —-1)+ DL (pi—1)(pj-1) +>: 


t=1; j=i+l 


+ (pi-—1)(p2 —1) +--+ (pe — 1)) 
= Pi po*** Pe— (Pi-—1+1)(p2—14+1)-+: (pe-—141) =0. 


This gives zero in every term of the first row except one. From the addition 
of the elements of the last column we have 


5[pi pe +++ Pp - 1 — D(m—1) — D(m—-1)(m—-1) --: 
— (pi — 1) (p2 — 1) «++ (pea — 1)) 
—e+Di(pi — 1) + D ie (pi — 1) (p2 — 1) +--+ = 0. 


Hence when the order of st is the product of r distinct odd primes pi, p2, +++ 5 Prs 
the degree of {s, t} does not exceed 


o+pu/(pi —1)(p2—1)-++ (p - 1). 


6. Let N = 2pq, where p and q are two distinct odd primes. The follow- 
ing five equations can be written down: 


from s and ¢ jointly, m+T,=4n, 

from st, 2m +To = p+, 
from (st)?%, 2m + Fo =pth, 
from (st)*%, To — qy2 — 2qz2 =u +t, 


from (st)’?, To — py: — 2p = w+), 
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where 5, h, i, and j are positive whole numbers or zero. The elimination of 
m, Yi » 21, and 23 gives 


[1 p-1 2p—1 2pq-1 w—(q—1)y2— (pq—1)ys— (2g — 1) 2 


2 p 2p 2pq w+ 6 — qy2 — pay. — 2422 

\2 0 2p 2pq w+ h — 2qz = 0, 
0 p 2p 2pq— wB+t— pays 

0 0 0 2pq w+j — qy2 — pays — 2G22 





which reduces to 
—utq(p—1)6+qh+ (p—1)¢4+(¢-—1)3+49(p¢q —2p +1) y% 
+ pq(q —1)y¥3s + 2¢(qd —1)(p—1)z =90; 


whence 6 = w/q(p—1). [If the order of st is twice the product of two distinct 
odd primes the degree of {s, t} does not exceed 


w+ p/[q(p—1)]). 


This case is needed in the proof of our theorem only when one of the two odd 
primes is 3. 

7. Now let N have at least two prime power factors the smaller of which 
is greater than 3. N = p*q® ¢, where p and q are distinct primes and ¢ 
is unity or a number relatively prime to p and q. We start from 


gh wt Be, 

The powers (st)*” and (st)*” give 
To(*,0)+To(*,*) =ut+h, 
To(0, *) +To(*, *) =u +2, 


where h and 7 are positive whole numbers or zero, and I’, ( * , 0) is to be under- 
stood as equal to Dy(Y — a) + Xz(Z—a), the summation extending 
however only to those constituents of {s, t} whose degrees are divisible by p* 
and not by q®. In I.(0, *) the summation extends to the constituents 
whose degrees are divisible by g® and not by p*, while in T,(*, *) the 
summation extends to constituents whose degrees are divisible by both p* 
and q?. Thus 


l'e = T.( x, *) +T.( x, 0) +1T.(0, *) +T,(0, 0). 
Then 


h+e: T ,0 r(0, 
pn yb HE PoC. 0) + P00, #) 


+To(x, *). 
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Multiply by (p* q® — 2)/(p* ¢* ) and add to 
pts = 2 —T2(*,0) —T2(0, *) —T2(*, *) — 12 (0,0), 
with the result (writing m for p* q° ) 
w—2h +1 r+2 
9 


7 ( 
a zie —_ T?. 


B+6 =pt 2u/ar — Mien +2) 


7 2r 


Here the new symbol I is introduced (k = 0,1,2,--+). It is to mean 
the sum [',(*,0)+T.(0, *) when k=1, ITa(*, *) when k& =2, 
T,(0,0) when k = 0. It will presently be used when more than two prime 
powers are in question. Then I’, for example, shall be the sum of all the 
quantities which, like T.(*,*,*,0,---) contain three stars in the 
parenthesis. Returning to the proof of the last result, 

wr—2 


aan = 
T'o(*,0) —T2(*, 0) ->(4r- y+2)y 


2r 


2a 
+2 a 
Soe E ~>53)? 


r+2 - 4m’ r+2 
~ 9 o(z-- 3)t=- or Tse (+2) (*, 0); 


9 
Z-%+2)s= -¥(*h*r-2)s 








| 
| 
to 
+. 
to 
| 
‘eu 
N 
Il 
| 
Sy 


and similarly 
er—-2_ 2 2 ow 
ar To(*, *) —T2(*, *) = —=T,(#, *) = —7f.- 
The summations > are to be taken in accordance with the symbol Ip ( * , 0) 
or I'o(*,*). Hence in this case n S wp + 2y/z. 
8. If g®@ = 3 and p* > 3, we may proceed as follows: 


4 . 
n =m —Tr+3,0(—m—h+To(*,0) +To(x, *)) 


= ; 
+P (—w—i+To(0, *) + To(*, *)), 


«a & 
n= —T2+*——(—p—h+To(*,0) +To(*, *)) 


4 
+—(—w—-itTo(0, *) +To(*, *)) 


wr—4 
Pi3(*,0) — - Dew w—ay (0, *) 


ll 
:— 
+ 

las 

| 

| 

| 

| 


9 
—-T,(*, * ) = T, (0,0). 
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Hence n is not greater than up + 2y/7. 

9. Next, let N have at least r (r > 2) Sylow divisors each of which is 
greater than 2: N = rq, where 7 = pf p?* --- p?, the product of powers of r 
distinct primes, and where ¢ is relatively prime to 7. From the r sub- 
stitutions (st)”™* ,k =1,2,---,r, we have r equations which when added 
together give 


— 1m — Do he + 2 RTP = 0. 


Multiply this equation by (2 — 2)/rz and add to 


n=2m—)>T?, 
k=0 
whence 
9 
n= p+ —_ 
T 
lfk =0, 


*, 2k —k) ve 2 
p ee TP eeetet — a he 


k=0 't 


2rar 
ak+(r—k)xr ~’ 
and ifk =r, 
2rr 
2k+(r—k)r ™ 
since r > 2, and m — 2 < rm/3, we have fork = 1, 


er 
2+(r—l)e~* 
When 1 <k <r, 


__2rar 7 2r ace ‘ 

Seir- bie *e-k* +1)<3-4-5--- tok factors. 
It follows at once that the greatest possible value of n in this case is up + 2/7. 
The results of the last three paragraphs may be stated thus: 

If the order of st has r (r > 1) Sylow divisors, no one of which is the first 
power of 2, and whose product we call x, then the degree of {s, t} does not exceed 
w+ 2y/7. 

10. There remains to be investigated the limit of n when JN is twice an 
odd number (7) in which there are at least two prime power factors. 

If t = p* - 3, we have already before us 


+ 


4 2xr—6 vy 
5 gatas : is Ton xn—s (0, *) 


9 
2u _ 
n= ae eae 
ae” T T T 


6 ; 
— = Tys(*, 0) -2- T,(0,0), 
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where z is the number of transitive constituents of degree 27 in {s, t}. Now 
( st)” gives 
2m + Fo =u t+), 


where j is a positive integer or zero; expanding this, 
Qrz + Fo(*,0) + Fo(0, *) + Fo(0,0) =p +z — 2m. 
Substitute this value of 2z in the expression for n above, and obtain 


2m 6 


n=p+e ~ + — — = Exp (*,0) —* Salt, *) — E,(0, 0) 


— 1 Fai. » (0, 0) 


— = Fars (0) — — Frets (0, #) — * 








— 6 4i } 
7 im i 
7 ° 


2u 
2m +ToSu+—7, 


as has been proved, while fp =u. Hence 


9 


9 
2m = = ’ —=s= 4 
T 1 
and therefore 
9 
mM , 2p 
nsapt+—+-;- 
TT 17 


11. When N = 27 is twice an odd number which is the product of two or 
more prime powers each of which is greater than 3 we may use the formula 


: 2k 
> + mn k) 


(k) 
Pore /(2k+(r—k)e] = 22, 


onp+B-2—2fa, - 


k=1 


where z is the number of transitive constituents in {s, ¢} of degree 27. From 
(st)” we have 


r—1 
‘ o 2 ¢ Xk) 
— 2rz = —p—j+2m+>F 
k=0 


(j a positive integer or zero), which when substituted for 2z in the preceding 


expression gives 


hw, 2m Bk+ar(r—k) wy 
ee a 
"> rr 


Sekar t a(t b) pw j 4-2 
_ 2re/(2k—r+an(r—k)] — — 





rr . T TT =I 
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It is necessary to show that 


a. ~<2+-3-5---tok+ 1 factors 
2k—r+a(r—k) ~~ 7 
If k =0, the difference Z — 27/(4 — 1) is positive since 6 is the least 
possible value of Z. If k = 1, we must have 2rr/(2 —r+rmr— 7) <6, 
that is 2rm — 38r — 37 +6>0, or 2(r —2)(4 — 3) +2>0, which is 
certainly true. Finally for1 <k <r, 


tT Tv r 


3k —++e(7—-k) “x —-1r-F 





1 
<(14+ 44) @40 
<2(k+1) <3* <3-5--- tok factors. 


Just as in the preceding paragraph we now draw the conclusion that 
wo, 2p 
— } }. at | —,. 
eer es” 


When the order of st is twice an odd number 1 divisible by two or more distinct 
primes the degree of {s, t} is not greater than wp + p/m + 2y/n’. 

12. This result is significant in the proof of the present theorem when 7 is 
the product of two or more distinct odd primes. In this event it gives a limit 
for n that is lower than the limit found when the order of st is TIj~; px, for 


1 2 1 
Wa, th ~*~ TW m n))’ 
Thai Pe + Tia Di TIi=i (pe — 1) 
In fact, 
oe oo = 
Thin: pe © TNies pi ~ (pi — 1) Tice pe’ 








It has now been proved that if the order of st is divisible by the odd primes 
Di, Pay «++, Pr its degree is at most w+ u/(pr — 1).(p2— 1) +++ (mw — 1), 
and that if its order is 2*, a > 1, its degree does not exceed pw + u/2*". 
Then the theorem follows as stated even when {s, ¢} has regular constituents 
in letters common to both s and ft. 

The general proposition that the degree of a group of class uw generated by two 
non-commutative substitutions of order 2 and degree u is at most 3u, is covered 
by this theorem.* 


*Bulletin of the American Mathematical Society, vol. 20 (1914), 
p. 469. 
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13. Let G be a p:p times transitive group, p an odd prime. By hypothesis 
there is in G a substitution s of order 2 and degree yp: 


8 = (Ap dz) (a3 a4) +++ (Ap-2Gp1) °°, 


and by reason of the fact that G is more than p — 1 times transitive, there 
is in G a second substitution s’ similar to s: 


8’ = (ao a3) (a2 a5) (a4 a7) «++ (Gps —) +++ (@p4). 


If the blank space that follows a,_; in s’ is filled by a letter new to s , the product 
ss’ has at least one cycle of p letters. Transform s’ by all the g,-: substitutions 
of the transitive subgroup of G that leaves fixed the letters ao, a2, «++ ap-1. 
Each of these gp: transforms s’, s’’, --- has in common with s the p — 2 
letters ao, @2, +++, @p-2, and x other letters. It is easy to calculate Doz for 
these gp: substitutions s’, 3’, ---. Any letter, not ao, a2, «++, Or Gp-4, 
is found as often in the set as any other. There are in the cycle (a,_3 — ) 
and in the last 34 —3(p—1) cycles of the substitutions of the set s’, 
3’, «++ (uw — p + 2) gp-1 such letters. Hence there are 


(u — p +2) Gp-1 
n—-p+1l 

substitutions among s’, s’’, --- that contain a given letter (not ao, a2, ---, 
ap—1). Therefore 

(u — p+ 2) 9-1 
Q) Le=(e-ptl)~ ay 
But in the set s’, s’’, --- there are exactly (nm — uw) g, substitutions which 
replace a,-3 by a letter new to s, and (u — p+ 1)4g, that replace ap; by a 
letter of s. Then 


=(u—pt+1)(u—pt+2)g. 


(p —2 
(2) XD (p—-24+2) = 2—)# in — gy, 
(n—2)9p Pp 
(3) D> (p-24+2) = w(u—ptl)g. 
(u—p+l)gp 


A word of explanation should perhaps be added. Since all the substitutions 
3’, 3’, +++ fix ap4, all are non-commutative with s. Moreover, when s’ 
(say) replaces ap; by a letter new to s, ss’ contains a cycle of degree p. It 
has been proved above that in the latter event s’ has at most u/(p — 1) 
letters new to s, whence it follows that the number of letters common to s’ 
and s cannot be less than wp — p/(p —1) = (p—2)pu/(p—1). Then on 
the addition of (2) and (3), 


—2 
~zr= (P _ re (n — 2) 9p + -—p+2 (u —p+1)gp; 
p-l 2 
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put for Dz its value from (1), 
p—2 
(uP +1)(H—P +2) op => 4 (uP +1)(n— 2) op 
+($-p+2)(u-ptio, 


or, on division by (u — p +1) 4g», which is not zero because yu > p,* 


—s 
B-pt+2= —j(n—H) +5-pt2, 


p 
Pp 
e=2 t= 7. 


pt, i#=5—1™ 


2 
whence finally, 
~ 2p — ‘.. 
& op — 5 





If p = 3, that is, if G is triply transitive, 


IIV 
role 


be 


It will be recalled that the holomorph of the elementary abelian group of 
erder 2* (a > 1) is a triply transitive group of class 2*', in which all aa 
substitutions of degree 2* are of order 2. 
14. If p = 5, this limit is 
w= Zen. 


Whenever u is not divisible by 12, this may be replaced by 





3n + 4 
ae 


IIV 


be 


For those substitutions s’, s’’, --- which replace ap_3 by a letter of s cannot 
have so many as p/2 letters in common with s unless such a product as ss’ 
is of order 6. The three equations 


m + 2y + 52 =p, 


2m + 62 =u, 
dy + 6z =u, 
give on solution z = »/12. Hence if yu is not divisible by 12, no product 
ss’, s3’’, --+ is of order 6, and in consequence 


ee (342) =($ +1) - 400 


* Jordan, Traité des substitutions, 1870, p. 64. 
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whence 
3n + 4 


0 


IV 


Mm 


15. If p = 7,no product ss’ , ss’’,, - - - is of order less than 8 when s’ (s’’, --- ) 
replaces a, by a letter of s. Then for the (u — 6)g; substitutions of the 


set 8’, s’’, --+ which replace a, by a letter of s, we have 
dy (5+) Ziu(u—6)q, 
; (u—8)97 
which added to 
D> (5+2) = hu(n—-yp)g, 
. a (m— 697 
gives on reduction 
w= ign. 


16. Let it be assumed that @ is more than p times transitive, p an odd 
prime. Then besides 


8 (do a2) (a3 a4) +++ (Ap Gp-1) +++ (ap) 


we have 


, 


8 


(do a3) (a2 a5) «++ (Ap-3 ap) -+* (Gp-4). 


‘ 


The transitive subgroup G, of degree n — Pp transforms s’ into a set of gp 
Mt 


substitutions s’, s’’, --- such that all the products ss’, ss’’, --- are of order 
p or a multiple of p. Hence 


p-2 
X(p-2+2) =e — gy. 
I Pp 
Any one of the last » — p + 1 letters of s is found in 


(u—pt 1) gp 


n— p 


= (4 —pt1) 9p 
of the substitutions s’, s’”, ---. Then 


, p-2 
(u— P +1) Gri1 = poe — (P — 2) op, 


or 
‘ p—2 
(u — P+ 1)" dos = (m — Pp) (w— P+ 1) gp, 
whence 
p—2 
= = oun - 
‘=. ‘lle OR 
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and for sextuply transitive groups, 
pint. 

17. The general case may now be very briefly stated. Let the group G be 
more than o times transitive, where o is the sum of r ( > 1) distinct odd primes, 
Pi, P2,°**,Pr- When the substitution s is given, s’ , because G, is transitive, 
may be so written that the first r cycles of ss’ are of the orders p1, po, +++, Dr- 
The substitutions s and s’ have {= (px — 2) common letters in these first 
5 Li=i (px —1) cycles, and x other letters in common. Summing for all 
the g, substitutions s’, s’’, -- 


D3E: = (u- (pe - 1) )? geri; 


and 


IT (pe -1) -1 


IIV 


Yo » 


> (Lim — 2) +2) 


II (pr — 1) 


whence, putting [| for II (gm ~ 1), 


(--£en—v)fs[Mps Lem ](o-En) 


whence 


r r 


. : (JI+1)-Lp x»n-Dp 
»-Emtret et (n-rm)+——p 1, 
w— 2 (pe-1) 
Since 


n— Dive > - Li (pe -1), 


r 


IT (x —1) 





Nin~1) 


1 


Here the condition r > 1 must be taken into account because the coefficient 
of (n — Diss pe)/( mu — Dit (pe 1) ) is zero whenr = 1. Then if a group 
is more than 8 or more than 15 times transitive, u exceeds {n or fin. In 
connection with this result it might be well to note that a closer study of our 
group G when 16 times transitive has shown that instead of the limit {4n, 
the class is actually greater than $$n. Since this appears to be an isolated 


result depending on the three primes 3, 5, and 7, I omit the proof. 
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18. If G is exactly o = pi + po +--+ + p, times transitive, pi, po, +++, Dr 
being distinct odd primes in ascending order of magnitude, r > 1, we can 
prove in like manner that 


(TI (p.-1)-1)n4r-1 


II (p: -1) +n -2 





u> 


19. Let G be doubly transitive.* The substitution s of order 2 and degree 
u( > 3) is one of w conjugates. Any transposition occurs exactly 


wyu/n(n — 1) 


times in this set. The number of possible transpositions in which one letter 
is displaced by s and the other is fixed by s is y(n — yp), so that if y is the 
number of substitutions conjugate to s which are not commutative with s, 


2wu (n — uw) 
n(n — 1) 


IV 


y 


Since yp letters may be paired in u(y — 1)/2 ways, the set has wu(u — 1)/2 
pairs, and wu(u —1)/n(n — 1) substitutions of the set have a given pair 
incommon. Hence if x is the number of letters any one of the w substitutions 
has in common with s, 


dx = wp?/n, > x(a — 1) = we? (wp — 1)?/n(n — 1). 


It is assumed that yu is less than n/2. Neither y nor w — y is zero. Now 


Ee-E(r- Ee) +5 oe, 


w 


p(2-#) -maee 


- n 


whence 








Since 
t= y/2, 
y 


2\2 2\2 
E(=-£)=0(3-£). 
; n 2 


and therefore 


2\2 ay? (n — yp)? 2\2 
> (2 -£) pea okt 27. 
oa» n n?(n — 1) 2 
* This page is a brief sketch of the reasoning given in full in the Bulletin of the 


American Mathematical Society, 2d series, vol. 20 (1914), pp. 471 ff., which 
see, as well as Bochert, |. c. 
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we - w\ PP 1 wefn 2 
o(2-#)e. x (« -f)\f=24 4 (5-«). 
so that 
Ky = 1 ¥ F( n ) 
—~w—un \2 tt 


n(n — p) _ fo~ pr. ’ 
~ Qy sar =(3-#)- 


On solution this inequality gives 


Now 


n?(n — 1) 


we (n — pw)? | v(5- 


whence finally 





p>an—ivn—-1. 
STANFORD UNIVERSITY 


August, 1916 
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SOME GEOMETRIC CHARACTERIZATIONS OF [ISOTHERMAL NETS 


ON A CURVED SURFACE* 
GABRIEL M. GREEN 


INTRODUCTION 
if a surface S, whose equations are 
x= 2z(u,e), y=y(u,r), 2=2(u,r), 


be referred to an orthogonal net of parameter curves, its first fundamental 
form is 


ds? = Edu? + Gdv’, 


7 


where FE and G are functions of u and v. This orthogonal net is called iso- 
thermal, if by a proper transformation i = @(u), do = ¥(v) of the parameters 
the first fundamental form may be reduced to 


ds? = \ (i, 0) (da? + di). 


A necessary and sufficient condition that such a reduction may be effected 
upon the form as first written is, that the equation 


(1 an ¢_) 
) Oudv ~AgI™ 


be satisfied identically. 

So far as the writer is aware, a purely geometric interpretation of this 

analytic condition does not appear in the literature. Of course the charac- 
terization of an isothermal net as a net which divides the surface into in- 
finitesimal squares is not a geometric one.t The fact that these nets have 
long played a very important part in geometric investigations may lend 
interest’ to a purely geometric characterization thereof. In what follows, 
we present two such characterizations, each of which has its own peculiar 
* Presented to the Society, December 27, 1916. 
{+ Nor can Kasner’s elegant characterization of an isothermal net, as one for which the 
two-parameter family of isogonal trajectories form a linear system, be regarded as geometric. 
Cf. his note, A characteristic property of isothermal systems of curves, Mathematische 
Annalen, vol. 59 (1904), pp. 352-4. 


480 




















ISOTHERMAL NETS 481 


advantages: one is very easily derived, while the other, although it requires 
more complicated analytic machinery, leads to a more elegant result. 

It should be noted, moreover, that the second characterization mentioned 
applies to all interesting cases, whereas the first breaks down when the iso- 
thermal net consists of the lines of curvature, that is, when the surface is 
isothermic. Darboux, however, has given a very elegant, though more 
complicated, geometric characterization of isothermic surfaces,* in connection 
with his solution of the problem of Ribaucour, in terms of the envelope of a 
certain two-parameter family of spheres. The related transformations D,, 
of isothermic surfaces, due to Darboux, permit the characterization of such 
surfaces in terms of the theory of rectilinear congruences. 


1. THE RELATION R. CONJUGATE CONGRUENCES 


Let X, Y, Z be the direction cosines of a line / in space. If X, Y, Z be 
given as functions of u, v, then through each point P (2, y, z) of the surface 
S passes a line /, and the totality of these lines form a congruence which we 
shall denote by ['. We shall assume that the line / which corresponds to a 
point P of S does not lie in the plane tangent to S at P, in other words, that 
the determinant 


Ox Ox ., 
du dv 4 
Oy Oy 
a du dv : 
Oz Oz , 
du ao ” 
does not vanish for any point of the surface 8. 
Let us write the equations 
xr dx Ox y 
au - "du + av Fea, 
(2) xr ox yr Ot 4eX 
m dude “aut” ae t © - 
ox 5, OX 5, OX sail 
oF aut a tee, 


the coefficients of which are functions of uw, v which shall be determined as 
follows. In the first of equations (2), for instance, replace 2, X successively 
by y, Y and by z, Z, thus obtaining three equations in ail, which may be 


* Sur les surfaces isothermiques, Annales scientifiques de l’école normale 
supérieure, 3¢ série, vol. 16 (1899), pp. 491-508, in particular p. 504. 
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solved for a, b, ¢ since the determinant A is nowhere zero. The coefficients 
of the second and third of equations (2) may be calculated in the same way. 
These equations are obviously fundamental in the study of the congruence T 
in its relation to the surface S. They are, in fact, very similar to the Gauss 
equations for the surface S, to which they reduce if X , Y, Z are the direction 
cosines of the normal to S. 

The equations 


(3) tector v= s45e C=2ty— 
Ou’ : Ou’ 
define the coérdinates of a point M(£,7,¢) on the tangent to the curve 
» = const. at P. If P be allowed to move along a curve u = const. on S, 
the corresponding lines PM generate a ruled surface R™ , and the point M 
traces a curve on this ruled surface. The tangent at any point of this curve 
traced by M has its direction cosines proportional to 0£/dv, dn/dv, O¢/dv, 
where for instance 
OF Ox . Op Ox x 


de ae * dv du? “aude’ 


Using the second of equations (2), we obtain 


. Ou ‘ du P Ox — 
(4) Fe = (Get aa’) Get + wb) Se + we x. 


The tangent to the curve traced by the point M as v alone varies lies in the plane 
determined by the lines | and PM if and only if 


(5) B= — FP 


We have therefore determined geometrically a unique point M on each 
line PM, which we shall call the minus first transform of the point P, and 
which is defined by the equations 


6 ‘. 1 Ox 1 Oy 1 dz 
=2-s> =y¥—-asas> =Z—-Tps-: 
(6) 7 b’ du’ 19 Wau b’ du 
In the same way, we may determine the first transform of the point P as a 
point M’ lying on the tangent to the curve u = const. at P, and defined by 
the equations 
1 Ox 1 Oy 1 dz 
7 f=x2-— > ‘= y- se i ms—-->. 
(7) . a’ Ov’ 7 Y~ a’ dv’ s a’ dv 


Our designation of the points M and M’ as the minus first and first trans- 
forms of the point P is intended to indicate the analogy to the Laplace trans- 
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forms of a conjugate net. In fact, if the parametric net on S is conjugate, 
c’ = 0, the second of equations (2) becomes one of the Laplace type, and the 
points M and M’ defined by equations (6) and (7) are actually the minus 
first and first Laplace transforms of the point P. 

If the parametric net is conjugate, the points WM and M’ just determined 
are the same no matter what the congruence [ may be. If, however, the 
parametric net is not conjugate, the points M and M’ change with the con- 
gruence I. The following considerations will show why this is so. For a 
non-conjugate parametric net, the ruled surface R™ formed by the tangents 
to the curves v = const. along a fixed curve u = const. is a skew ruled surface, 
and therefore a plane through a generator thereof is tangent to the ruled 
surface at a definite point of that generator. From the way in which the 
point .V was obtained, it is obvious that the point M 1s the point at which the 
plane determined by the lines | and PM 1s tangent to the ruled surface R. A 
similar characterization obtains for the point M’, and so different points M 
and .M’ would be determined by different congruences ['. 

The line /’ which joins the points M and M’ we shall call the line conjugate 
to the line 1, and the congruence I’ which is formed by the lines /’ the con- 
gruence conjugate to the congruence T. If the parametric net is non-conjugate, 
it is easy to see that a given congruence I’ can have but one congruence IT 
associated with it, so that the congruences [ and I’ may in this case be de- 
scribed as conjugate to each other. For, the line /’ intersects the parametric 
tangents in points M and M’, and the line / is then determined as the line of 
intersection of the two planes tangent at M and M’ to the two parametric 
ruled surfaces R™ and R™ respectively. If the parametric net is conjugate, 
however, a congruence I’ will have a unique conjugate congruence I’, but a 
congruence I’—of which there exists only one for a given conjugate net— 
will be conjugate to any congruence I’. In the case of a conjugate net, 
Wilezynski* has called the line joining the minus first and first Laplace trans- 
forms of a point P of the surface the ray of the point P, and the congruence 
formed of these lines the ray congruence. 

We have elsewheref spoken of the relation between the congruences T 
and I” as the relation 2. It is obviously a purely projective relation, and 
we have in fact found the notion of the greatest importance in the general 
theory of surfaces and rectilinear congruences. 


* KE. J. Wilezynski, The general theory of congruences, these Transactions, vol. 16 
(1915), pp. 311-327. 

t In three communications to the Society, abstracts of which have appeared in the Bulle - 
tin as follows: vol. 21 (1914-15), pp. 484-5, vol. 22 (1915-16), p. 274, and vol. 23 (1916-17), 
pp. 73-4. 
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2. THE DEVELOPABLES OF THE CONGRUENCE I” 


The congruence I’ formed of the lines MM’ contains two one-parameter 
families of developables, which we now proceed to determine. Any curve on 
the surface S may be defined by equations of the form u = u(t), » = v(t). 
We seek those functions u(t) and v(t) which define the curves on S which 
correspond to the developables of I’. We note first of all that if the line MM’ 
generates a developable, then the tangent to the curve traced by M, the 
tangent to the curve traced by M’, and the generator MM’ lie in a plane. 
The tangent to the curve traced by M has its direction cosines proportional 
to dé/dt, dn/dt, dé/dt, where £, n, ¢ have the values given by equations (6) 
when wu and v are replaced by the functions u(t) and v(t) respectively. We 
have of course 

dé od& du . dE dv 
dt du dt * dvdt’ 


Og ; a oOf1 Ox b dx c y 
Ou ~ b’ du\b’ du wb’ dv b’*’ 
Ou ¢’ 


Life 
—=-—-z,|—logb’+a jJ— - TX. 
Ov b’\ av '°® du ib’ 
Equations (8) are obtained directly from equations (6), use being made of 
the first two of equations (2). We have therefore 


ff, -% a5.) 1(@ ey aS 
dt — — PT ou dt 4” a4 dt | du 


b dudx 1 ( du i) ‘ 
—~—>>—slest+es )X, 


etc., 


where 


(3) 


~ b dt dv b’ dt dt 


the same equation subsisting, of course, if & be replaced successively by 7 
and ¢, and x at the same time by y and z, XY by Y and Z. Similarly, 


dé’ a’”’ dv Ox { ( b’’ d a 


dt ~— a’ dt du t | t a” ov jdt 
1 (@a’ pr du|d@x 1 ,,dv , du y 
~ a? \ du ves dt {av a’\° dt sil a 
We have also 
10x 1 Ox 
2 ae x wea oe te 
"> b’ du © a’ dv’ saan 


We wish now to write a necessary and sufficient condition that the line 
whose direction cosines are proportional to dé/dt, dn/dt, d¢/dt, the line 
whose direction cosines are proportional to dé’/dt, dn’/dt, dt’/dt, and the 
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-/ 


line whose direction cosines are proportional to & — £’, » — n’, ¢ — ¢’, be 

parallel to the same plane. The condition sought is sali seen to be the 

vanishing of the determinant formed from the nine coefficients of d2/du, 
-s 


0x/dv, X in the above expressions for dt/dt, di’/dt, — — &’. When this 
determinant is expanded, the condition is found to be 


du du dv (dv\? 
(9) 4(G) +23 at° (|) rm 


where 


bh’ 
[ow at Fab] 4 (Sa v), 

Ou 

, 7 i 
C= -e| aa —1 y") +2 -ary |e" (S40 b’ ), 
ob’ ; ,{ da’ db’ 
{o” |» (b tree a b] + (* - 5) 
— e| a’(a’ - vy + aro | I. 


Equation (9) is a quadratic differential equation of the first order, the 
solutions of which determine the curves on the surface S which correspond to 
the developables of the congruence I’ of lines 14M’. The asymptotic net 
of the surface S is moreover defined by the quadratic differential equation 


du yy de dv ” dv -_ 
(10) v(S ny +2D'= 7 +D (7) = 0, 


where D, D’, D” have the usual significance as the coefficients of the second 
fundamental form of S. Each of the quadratics (9) and (10) determines a 
pair of directions at any point of the surface S, and these pairs separate each 
other harmonically if and only if 


(11) AD” — 2BD' + CD = 0 


A 


bole 


B 


The left-hand member of this equation is the simultaneous invariant of the 
two binary forms (9), (10). But two tangents are conjugate if they separate 
harmonically the asymptotic tangents. Therefore the net of curves defined 
by equation (9) is a conjugate net if and only if equation (11) holds. 

Now, the coefficients c, c’, c’’ of equations (2) are easily seen to be pro- 
portional to D, D’, D’’; in fact, 


e:D=a¢c:D= D” = VEG — F*/A, 


where A is the determinant introduced at the beginning of § 1. If the values 
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, 


of c, ce’, c’’ be substituted in the expressions for A, B, C, and the resulting 


quantities then substituted in equation (11), this equation reduces essentially to 


9 mw Tr ee a eee 
(12) (DD”" — D ($ - 


da’ x) 
The first factor on the left vanishes identically only for a developable 
surface. But for such a surface the only conjugate nets are nets one com- 
ponent family of which consists of the straight line generators. We may 
therefore state the following theorem: 
TuEeoremM. The developables of the congruence Y’ of lines MM’ correspond 
to a conjugate net on the surface S if and only if 


da’ Ob’ 


du av 


(15) 


If the parametric net is conjugate, equation (13) expresses the fact that 
the net has equal Laplace-Darboux invariants. The characterization of a 
conjugate net with equal invariants as a net for which the developables of 
the ray congruence correspond to a conjugate net on the surface, is due to 
Wilezynski.* 

We have found a geometric property which is expressed by equation (13) 
when the surface S is not developable. We proceed now to obtain another 
geometric property which covers the case of a developable surface, but which 
on the other hand becomes trivial if the parametric net on the surface is con- 
jugate. 

The totality of points W(£,7, ¢) form a surface S_;, which we may call 
the minus first transform of S, and the points M’(£’, 7’, ¢’) likewise form a 
surface S,, the first transform of S. From the way in which the point M was 
obtained, it is clear that the plane of the lines / and PM is tangent to the 
surface S_; at M. Consequently the tangent at M to the curve u = const. 
on S_, intersects the line / in a point whose coérdinates we shall now find. 
The said tangent has direction cosines proportional to d&/dv, dn/dv, d¢/dv, 
these derivatives being given by the second of equations (8): 


a4 a ae, ee 
4) av——sé as? +e du —b’** 


The required point of intersection must have for its x-codrdinate an expression 
of the form £ + vdé/dv, and also of the form 2 + wX. By combining equa- 
tion (14) with the first of equations (6), it is readily seen that the point of 


* Loc. cit., p. 319. 
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intersection has for its 2-coérdinate 


1 0é c’ : 
(15) é _ a, os = £ a F X fs 
5, 08 b’ +a’ b’ (<: log b’ + a’) 


Similarly, the tangent at M’ to a curve v = const. on S; meets the corre- 
sponding line / in the point whose 2-coérdinate is 


| eg cy 
a) * £ : — ) ine 
= log a’ + b a’| =-log a’ + b’ 
Ou ou 


Examination of the right-hand members of equations (15) and (16) shows 
that the points in which the line | is cut by the corresponding tangent to a curve 
u = const. on S_; and the corresponding tangent to a curve v = const. on S; 
coincide if and only if 


If the parametric net on S is conjugate, c’ = 0, and condition (13) is not 
at all necessary. In fact, in this case the two points of intersection on / coin- 
cide with the point P, always. 


3. GEOMETRIC CHARACTERIZATION OF ISOTHERMAL NETS 


We shall now suppose that the parametric net on the surface S is an ortho- 
gonal net, and that the line / is normal to S at P. Then equations (2) become 
the Gauss equations of the surface, and in addition F = 0. Some of the 
coefficients in the Gauss equations are 


c= J), ce’ = Pp’, ce’ =f)", 


p 1a ws , 10a 
a = 95,08 E, b = 55,108 G. 


The condition that the parametric net be isothermal, i. e., equation (1), may 
| b 
therefore be written 
da’ ab’ 


Ou —s- av’ 


which is precisely the equation for which we have already found geometric 
interpretations. Using the results of the preceding paragraphs, we may 
therefore state the following geometric characterizations of isothermal nets: 

An orthogonal net on a non-developable surface is an isothermal net, if and 
only wf the developables of the congruence which is conjugate to the congruence of 
normals of the surface correspond to a conjugate net on the surface. 
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If an orthogonal net on any surface S is not conjugate—. e., does not consist 
of the lines of curvature of S—then the net is isothermal if and only if each normal 
to the surface is cut in the same point by the corresponding tangent to the curve 
u = const. on the minus first transform of S, and by the corresponding tangent 
to the curve v = const. on the first transform of S. , 

If the lines of curvature form an isothermal net, the surface is isothermic; 
therefore, in virtue of the theorem of Wilczynski, a surface is isothermic if 
and only if the developables of the ray congruence of its lines of curvature correspond 
to a conjugate net on the surface. 


Harvarp UNIVERSITY 
March, 1917 




















A NEW METHOD IN BOUNDARY PROBLEMS FOR DIFFERENTIAL 


EQUATIONS* 


BY 
R. G. D. RICHARDSON 


INTRODUCTION 


One of the most important problems of mathematical physics is the study 
of boundary problems for differential equations of the second order in one 
or more variables. These equations contain a parameter which must be so 
determined that solutions exist satisfying prescribed linear homogeneous 
relations between the values of the function and its derivative on the boundary. 
The problem in one dimension has been studied in great detail.f As typical 
of the problems in one independent variable we have the equation 


(1) Ver +([r(a) +AK(a2)] e(x7) = 0, v(0) = o(1) = O. 


By means of his theory of integral equations, Hilbertt extended the theory 
of boundary problems in one dimension to the corresponding problem in 
partial differential equations of the elliptic type, of which 


(2) Ure + Uy, + (r(a,y) +AK(2z,y))u=0, u = 0 on boundary, 


is typical. By hypothesis, for both equations (1) and (2) the function & has 
one sign (orthogonal case) or k has both signs and r = 0 (polar case). It has 
been shown that there exist an infinite set of real parameter values \; for 


* Presented to the Society, September, 1914. 

t For existing methods and literature of the subject see Bécher, Encyklopaedie der mathe- 
matischen Wissenschaften, I1A7a; Boundary problems in one dimension, Proceedings of 
the International Congress of Mathematicians (Cambridge, 1912), and 
Lecgons sur les méthodes de Sturm (Paris, 1917); Lichtenstein, Zur Analysis der unendlich- 
vielen Variabeln, Rendiconti del Circolo Matematico di Palermo, 
vol. 38, p. 113. 

t Grundziige einer allgemeinen Theorie der linearen Integralgleichungen, p. 58 (Teubner, 
1912). This book is a reprint of a series of articlesin the Gé6ttinger Nachrichten. 

Refinements concerning the nature of the field R, of the coefficients of the equation and 
of the functions f (x, y) to be expanded in terms of the set of solutions constitute the chief 
progress made in this field since Hilbert published his work. See the literature cited by 
Lichtenstein, Randwertaufgaben der Theorie der linearen partiellen Differentialgleichungen 
zweiter Ordnung vom elliptischen Typus, Journal fiir Mathematik, vol. 142, (1912) 
wy i. 
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which solutions of (1) and (2) exist satisfying prescribed boundary conditions. 
In the orthogonal case this set has a limiting point at positive infinity only 
and in the polar at both positive and negative infinity. 

The boundary problem for partial differential equations of the hyperbolic 
type so far discussed has been of an entirely different nature. The existence 
of a unique solution has been shown* in case the value of the function is given 
on one of two intersecting curves and the value of the normal derivative on 
the other, neither curve being cut twice by a characteristic. The elliptic 
and hyperbolic types of equations are essentially different and are char- 
acterized by Du Bois Reymond as “ himmelweit verschieden.” On the 
other hand, in certain special cases at least, a boundary problem for the latter, 
analogous to that for the elliptic case can be solved, and functions f(a, y) 
may be developed in terms of the solutions corresponding to the various 
characteristic numbers. A simple example is given by the solutions of the 


equations 
[3(a)] Ure + Uyy + Au = 0; 
[3(b)] Ure — Uy, +ru = 0. 


If we demand that solutions of these equations shall vanish on the boundary 
of a square whose opposite vertices are (0,0), (a, 7) we find that solutions 
sin mx sin ny (m, n positive integers) of one equation are solutions of the 
other. For the elliptic equation the characteristic numbers must be of the 
form m? + n’, for the hyperbolic, m? — n?. In the latter case the set has 
limiting points at 0 and at + « while in the former this occurs at + © only. 
If the field is taken with sides of incommensurable length there may be an 
infinite number of limiting points for the set of characteristic numbers of 
[3(b)].¢ 

Including [3(a)] and [3(b)] as special cases there is a class of partial differ- 
ential equations for which known methods for ordinary equations suffice. 
For, if in the unit square we consider the equation 


= te)  @ Ou 
m(y) 5x] v(x) $e | — s(x) 5 | aw 5. | 
+[r(a)m(y) —s(a)l(y)]ut+rAs(x)m(y)u =0 


for which solutions are to vanish on the boundary, we know that the functions 
u(x,y) =v(2)w(y) and the corresponding parameters Ann, = Mn, — Yne 


* For the literature of the subject and the most general theorem see Mason, Selected Topics 
in the Theory of Boundary Value Problems of Differential Equations, The New Haven Mathe- 
matical Colloquium (Yale University Press, 1910). 

{ For example, if the opposite vertices are (0,0), (+, z/V2), the characteristic numbers 
are of the form m? — 2n?. Among the integers which occur among this set, each an infinite 
number of times, are all those of the form H? — 2, where H is an integer. 
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will fulfil the conditions, v(x) and w(y) being solutions of the equations each 
in one variable 


d dv(x)], — 

_~ p(x) — +[r(xv) +yus(r)] r(x) = 0, 
p>0,s>0, 0(0) =r(1) =0, 

d _dw(y) 

a1 49(y) rm +[l(y) +ym(y)]w(y) = 9, 


dy 
q>0org<0,m>0, w(0) =w(1) =0 


and fn,, Yn, the parameter values corresponding to solutions which vanish 
m, — 1, m2 — 1 times respectively. This is immediately evident on multi- 
plying the first equation by m(y)w(y) and the second by — s(2)v(a2) and 
adding. If g > 0 the equation is hyperbolic; if g < 0, elliptic. 

None of the methods at present in use lends itself to the discussion of the 
new boundary problem of the hyperbolic equation. The objects of this memoir 
are (a) to develop a new method by which the well-known results in regard to 
boundary problems for the ordinary and the elliptic partial equations can be 
deduced and (b) to attempt an application of this method to an investigation of 
the facts concerning the new boundary problem for the hyperbolic equation. 

The new method consists in dividing the region into equal rectangles by a 
lattice and studying the difference equations approximating to the differential 
equation at the lattice points. The size of the lattice division is decreased 
indefinitely and it is shown (§ 3) in all cases (including the hyperbolic) that 
each set of corresponding parameter values \ has a finite limit. The set of 
numbers which are limits of the characteristic values of the difference has, 
both in the elliptic and hyperbolic cases, a limit point at infinity. In the case 
of the elliptic partial differential equation it is shown (§ 5) that the corresponding 
functions of the lattice approach uniformly a function u(a,y) which is con- 
tinuous and that this limiting function has continuous derivatives of as high 
order as those possessed by the coefficients of the equation. The totality of such 
functions furnishes the characteristic solutions of the equation. 

The results obtained in the hyperbolic case (§ 6) suggest that with the 
possible exception of special cases, the boundary problem similar to that for 
the elliptic case may be possible of solution. However, the reasoning is not 
complete and the problem is left in large measure unsolved. 

Although the problem in one dimension is not treated by the new method, 
its discussion is naturally simpler than that of the elliptic equation which is 
discussed in full. All the ordinary theorems both for the orthogonal and polar 
cases of (1) can be derived in simple fashion by the new processes. The method 
can also be extended (a) to problems with more general homogeneous boundary 
conditions and (b) to problems in three or more dimensions. 

Beside a completion of the discussion of the hyperbolic case, other new 
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and interesting problems of various kinds are suggested by the results obtained. 
No problem in quadratic forms or in linear equations in an infinite number of 
variables in which the set of characteristic numbers has limit points both at 
infinity and in the finite region has as yet been solved. We can however 
set up by a method analogous to that used by Lichtenstein* for ordinary 
differential equations, problems of this nature equivalent to the hyperbolic 
case of (4) for which we know solutions exist. The hyperbolic differential 
equation may be regarded as the Lagrange equation of an isoperimetric 
problem in the calculus of variations in which there are an infinite number of 
integral sub-conditions (see footnote, p. 496). Such a problem has never been 
studied. The expansion of linear forms in terms of the solutions of the 
approximation equations (§ 1, III-IV) suggests an expansion of functions in 
terms of the characteristic functions of our equation. The author hopes 
later to report progress along some of these lines. 


1. THE LINEAR ALGEBRAIC EQUATIONS OF APPROXIMATION 


Let us approximate by algebraic equations to the most general linear self- 
adjoint homogeneous differential equation of the second order in two inde- 
pendent variables which may be written 


a Ou(x, y) ¢ Ou(x,y) 
Ox pP (x, y) Ox = oy q (a ’ y) dy 


(9) 
+[r(x,y) +k (x, y)]u(az,y) =0, 
where \ is a parameter. The notation will be simplest if the region to be 
considered is a unit square with opposite vertices (0,0), (1,1) and there 
are mp equal rectangles formed by m — 1 lines parallel to the y-axis and wp — 1 
lines parallel to the z-axis, the values of the functions p, q, r, & at one of the 
lattice points with coérdinates (i/m, j/u) (4 =0,1,-+--,m; 7 =0,1, 
-++,) being denoted by pi;, qi;, rij, kij, and the derivatives u,, u, being 
replaced by m (ui; — Ui-1j), bij — Uij-1) respectively. In this way we 
avoid the complexity of notation that arises in the general region for lattice 
points near the boundary.{ For the algebraic equations it will be assumed 
* Loc. cit. 
+ More generally the plane may be divided into rectangles by lines parallel to the axes, 
the lattice points denoted by (7,7), the derivatives du / dx at these points replaced by 
(uiz — ui-sj) / ly where 1; is the length of the interval between i — 1, 7 andi, 7, the derivative 


0? u / dx? replaced by 
2 (s+ —Uj Uj Mint) 
’ 
liga + l; liga l; 


ete. For points near the boundary the lengths of the intervals cannot be made uniform and the 
equations set up for these points must be of irregular form. While in the following discussion 
details are given only for the case that the region is a square it should be borne in mind that 
this restriction is due to a desire to retain simplicity in the notation and that the results hold 
for more general regions. 
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that p,q, 7, and k are continuous functions* of x and y and (to insure the para- 
meter values being real) that k is positive throughout the region. For boundary 
conditions we postulatet that u;; vanish on the contour, 


(6) Uo; = Umj = Un = Uj, = O 


(¢ = 0, 1,2,.°**, mj 7 @ G6, 1,2, +**, 0). 
There are then (m —1)(u — 1) homogeneous equations in as many un- 


knowns of the type 


L (ij) + Aeij uy = mM | pigs (Migay — Ui) — piz (May — UiR-1;) | 


(A) = Me | dijna (Maja a Ui; ) — qij (Vij ~~ U;;-1) | 
+ (rij + Nkj;) ui; = 0 
(t7=1,°---,;m—1lsj=1,°--,u—1). 


If we hold fixed the lattice point in question and indefinitely increase the 
number of divisions (in this paper we shall consider in general the set obtained 
by repeatedly doubling m and yw) the limit of (4) gives the equation (5) for 
this point. 

Intimately connected with equations (A) are quadratic forms 


1 . 
Q: = aa m* d Piti ( Ui41) _ Ui; )° 
(7) — we ) Vij+i ( Uij+1 — “:;)° = >» Vij ui; ], 
ij ij 
I 
Qo = — Dd hi; ui; 
= ET kus, 


in variables u;;.. From the hypothesis on k (2, y) it follows that Qs is definite. 
In seeking the minimum of Q, subject to the restrictions (6) and Q2 = 1, one 
is led on differentiation with regard to each of the variables u;; to the system 
of equations (A). Another relation under the conditions (6) is given by the 
formula 

1 
(8) Qi = —— DL (wis) wis = Qa, 

mu G 
which is established by multiplying the equations (A) by u;;, adding, and 
dividing by mp. 

In order that there be solutions different from zero of the (m — 1) (yu — 1) 
linear equations (A), \ must be a root of the determinant formed from the 
4 Eventually the important cases are those where p and q have each one sign. When they 
have opposite signs (p > 0, g < 0) we shall call the difference equations elliptic and when 
the same sign (p > 0, q > 0) hyperbolic. 

+t More generally it might be assumed that on the boundary hj; ui; + li; An (ui; ) = O where 


A, (ui) is a first difference expression corresponding to the normal derivative, and h,;, li; 
are sets of constants. The discussion in that case would be similar. 
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coefficients of that set of equations. This symmetric determinant has 
(m — 1) (mu — 1) roots for \, and from our hypothesis it follows that these are 
all real.* For, since \k;; occurs only as part of each term in the main diagonal, 
we may multiply the row and column in which it stands by 1/vk;; without 
altering the symmetry of the determinant and thus reduce it to the form of the 
secular equation.f Corresponding to each of these (m — 1)(yu— 1) values 
of X, called characteristic values, there is a solution u;; of the linear equations 
(A). These solutions shall be said to be normalized when the u’s are multi- 
plied by the constant which makes Q. = 1. The normalized solution corre- 
sponding to the characteristic \“ and which consists of a set of constants 
sach defined for a point (7/m , j/u) of the lattice, shall be called a characteristic 
function and denoted by U}. These sets of solutions U) and the values 
A” [x = 1, +--+, (m—1)(p —1)] have various properties, some of which 
we shall now develop. 

I. Characteristic functions US}, UY) corresponding to two values of are 
orthogonal, that is 


a Tix) JT) — 
(9) >> kis US) UY = 0. 

D 
To prove this consider a set of equations (A’) similar to the set (A) except 
that pi, ij, Tis, kis, Wij, XN are replaced by pj; , qi;, 7i;, kis, uiz, respectively. 
Now multiply the equation for u;; in (A) by uj; those for u;; in (A’) by — uj; , 
add and sum over 7j. There results under condition (6) the fundamental 
relation 


1 , — ; 
mg Oe (P85 + Nis — Fis —D hej) Ug; Ui; 
m 


(10) = 7 Do (pig — pis) (wig — win-ag) (wig — a3) 
yj 


M , 
+ >» (Gis — Gis) (ig — Uiz1) (Mis — Wid). 
ij 


On setting ki; = hij, rij = Tis, Dis = Piss Vig = Viz We have for two solutions 
Us), UW} of (A) the relation 


js 


(A wea ”) Zz. kj; U®) UY) = (0) 
ij 


*In the problem corresponding to the elliptic differential equation where Q, is definite 
(p >0,q <0,r <0), the roots \ are of course real even if the k’s have both signs. But if 
neither of the two quadratic forms Q,, Q2 is definite, at least a part of the \’s and of the solu- 
tions will in general be complex. Similarly if characteristic functions exist for the hyperbolic 
case (p > 0, q > 0) when k takes on both signs, at least a part of the characteristic numbers 
will in general be complex. The minimum and calculus of variations problems corresponding 
will then have no meaning. 

Tt Scott and Matthews, Theory of Determinants (2d ed.), p. 147. 
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and when “™ +, this gives the desired orthogonal property. If X is 
multiple valued (A = ™), the solutions U's), UY) may by the usual 
methods be so chosen that they are orthogonal to one another as well as to 
the others. The set of (m — 1)(u — 1) solutions thus chosen as orthogonal 
and normalized is called the fundamental set. 

II. The fundamental set UY} [x = 1, 2, ---, (m —1)(u — 1)] comprise 
(m—1)(u—1) linearly independent sets of eck —1)(u—1) constants. 
For, if there exist constants c, such that >, ¢, US) = 0 for all ij, on multi- 
plication of each fundamental set by /;; UY? and summing over 7 and then 
over x we get by using the orthogonal sepa the equation c, = 0. This 
equation holds for all values p = 1,2, ---, (m—1)(u—1). 

III. Any set of (m —1)(u—1) oseinaaie ¢;; defined for the lattice 
points is expressible in terms of the sets of constants Uf) as follows: 

— . Tix) 

in Li ee | _ (ce = 1 24,1 : sv). 
For, denoting by R;; the remainder ¢;; — X,¢, Ui; and multiplying by 
k,; U} and summing, we get Dis ki; Rij UY = 0 for ail «. Since the 
(m — 1)( — 1) sets of constants U% are linearly independent by II, the 
(m — 1) (pu — 1) multipliers k;; R;; must all be zero. And since by hypothesis 
k;; > 0, it follows that R;; is zero. 

IV. The arbitrary set of constants ¢;; is not generally expressible linearly 
in terms of v of the U;;’s [vy < (m — 1)(u — 1)] which we denote by U®. 
However, if we define the coefficients c; to be 


cz LF gis kes UP, 


~ my 5 
iy if by best 
approximation is understood that which furnishes 1;; ki; [¢i; — L;¢; U8 P 
with a minimum. 
For, in order to obtain a minimum for this last expression we equate to 


the best approximation to ¢;; is given by the formula ), ec, U“? 


zero the derivatives with regard to c; and obtain the relations 
’ 7) 7@ — 
2d kis (bis — De, UP) UG = 0 
y n 
for all k. From the orthogonal and normal properties it follows that this 
may be written 


— Fay kis UF 


mm 7 


and on taking the second derivatives these relations are found to be sufficient 


to establish the theorem. 
It follows from III that ¢;; — L;¢, U" = XL, ¢, US where « ranges over 
Trans. Am. Math. Soc. 3% 








496 R. G. D. RICHARDSON [October 


the remaining (m — 1)(u — 1) — v of the x«’s and e, is defined in a manner 
analogous to c,. Hence 


¥ a , 
— Do kis bi5 — Lo UGP = 2G. 


mp GF 


And in particular if k is a null set, 


1 L. _ 2 
ep ki; bi; = > CK 

V. Any characteristic number can be considered both as the minimum and the 
maximum for certain problems. If the (m — 1)(u — 1) numbers are arranged 
in a system dy S do S +++ Ne Ss SA S++ SAwm-yu-y,» then X, is 
the minimum* and maximum respectively for those values of u;; which make 
(2 = 1 and satisfy the relation (6) and the relations 


[11(a)] } ki; US) u;; = 0 (« =1,-++,8—1); 
y 

[11(6)] Do kis US wig =O («=841, +++, (m—1)(u—1)) 
y 


respectively. 
Since proofs for all cases are similar, we shall content ourselves with a dis- 
cussion of the last case. From III and relation [11(5)] it follows that 
8 ' 1 
9) — . (x) iS _ («) 
(12) ui; = p> c, l ij (c. ~ mu > hes Uy U5, ) 
on 


and hence that Q. = ci + ¢3 + --- +c? =1. From (12) and (A) we have 


L® (ui;) = ode ¢ hig UY) 


«=1 


and hence from (12) and (9) 
= — > L(uj)u;=mMetrmeeg+--- +c. 


Then Q: — As = (Ar — Ay) C7 + +++ + (Ag — Ay) 2-1 which is zero or nega- 


* As an extension of this result it seems probable that any characteristic function u (2, y) 
of the hyperbolic differential equation (5) (p >0,q > 0) may be regarded as the solution 


. ° ee ° ° : Y 2 
of the problem of finding a minimum and maximum respectively of the integral (pe 
. 
- qui, — ru’) dR for those values of u which vanish on the boundary, satisfy f kwdR =1 
. Re 


and an infinite number of orthogonal sub-conditions of the form J, kU; udR = 0, where z 
R 


ranges over the whole set of solutions corresponding to values of \ less and greater respectively 
than the one considered. The minimum or maximum value would be the corresponding 
characteristic number \}. The minimum property can be easily established for the elliptic 
case (p > 0,q <0), since here the number of sub-conditions is finite; it is probable that 
the maximum property holds also. 
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tive, and zero when ¢; = ¢: = +--+ =¢,-; = 0. That is, Q; is a maximum 
for UY) and takes on value X, in this case. 

VI. The solution U%) of V furnishes the maximum and minimum 1/A, 
respectively for Q2: when those sets of constants u,;; are considered which 
satisfy the relation Q; = 1, the relations (6) and the relations [11(a)], [11(b)] 
respectively. 


VARIATION PROPERTIES OF THE CHARACTERISTIC NUMBERS OF THE 
ALGEBRAIC EQUATIONS 


The functions p(x, y), q(t, y), r(@,y), k(x, y) may be regarded as 
analytic functions of a parameter ¢ in an interval -1=t=1. If in the 
discussion of I, § 1 we allow p’, q’, r’, k’ to approach p,q, r,k, the param- 
eter \’ approaches \ and the formula (10) becomes 


dr 1 dr;; . rN dk;; , 

dt a hij uiy = --SF ’ ha ee ot daa 

(13) : ‘. , 
- oP ii ae ee. pe aq; re . 
—e dt al * (Ui Ui—1;) =e dt we > (Ui; Uij-1) e 


From this formula one can determine the method of variation of the param- 
eter \ with the variation of one or more of the functions p, q, r, and k. 

I. In the set of equations (A) all of the characteristic values \ are increased 
aif one or more of the p’s are increased or one or more of the q’s decreased and 
all are decreased if one or more of the p’s are decreased or one or more of 
the q’s increased. 

II. In the set of equations (A) all of the values || are decreased if one or 
more of the k’s are increased. If \} = 0, the new characteristic value is also 
zero. 

III. In the set of equations (A) all of the values \ are decreased if one or more 
of the r’s are increased. 

It may be observed that since 


aay bi ui; = = (2 = 


my 
the change of \ is an infinitesimal of the same order as the mean change of 
the r’s and k’s. And since neither 


] 
2m? Us; — Uj-1;)” or - 2 (us; — Uj)" 
my ( ij ~13) n mu LH ( ij ij—1) 


can be zero when Q2 = 1, the a of X is of at least as high order as the 
mean changes of the p’s or q’s 
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3. FINITENESS OF CHARACTERISTIC NUMBERS 


The comparison theorems of the preceding section may be used to determine 
intervals within which the values of \ for all m, uw must lie. In the unit 
square R with opposite vertices (0,0), (1,1), constants Pi, Po, Qi, Qe, 
R,, Re, may be chosen so that 


P,=2=Pr, USE, Ry sph. 


aeedt Thee 


From I-III of § 2 it follows that the characteristic values of (A) are inter- 
mediate between those of 


m? Py (visas — Qui; + Uaay) — w? Qi (Mijn — Qui; + Uij-1) 
+ (Ri +d)uiz = 0, 
m® Po (Uina; — 2g; + Uiaasy) — we? Qe (tiggs — Qui; + Uij-1) 


+ (Re +) uy = 0. 


(14) 


(15) 


That the characteristic numbers of (14) and (15) are the differences between 
the characteristic numbers for the sets of equations 


m® P; (vig, — 20; + 1-1) + (2 + 1); = 0, 


(16) ' 
we? Qs (wy — 2w; + wri) +m w; = 0; 
m Po (41 — 20; “fe Vi-1 ) oo (Rs -t- V2) Vj a 0, 
(17) 


we Qs (wyi1 — 2wj + wii) +; = 0, 


respectively, where v;, w; are each sets of variables in one dimension and 
V1, Vi, ¥2, ¥2 are parameters, may be seen on multiplying the first of each 
set by w; and the other by — »;, adding and setting u;; = o;w;. The Xd of 
(14) is vy; — v; and of (15) is v2 — vp. 

It is easily shown that characteristic numbers of (16) are approximations 
for the corresponding differential equations 


(18) Pitzz(t@) + (i +n)o(x)=0, (0) = (1) = 0, 


(19) Qi Wy (y) + w(y) = 0, w(0) = w(1) =0, 
each in one variable and with constant coefficients. These have solutions of 
the form v = sin n, Tx, w = sin ne Ty (Mm, NM. = 1, 2, --+), which oscillate 


nm, — 1 and m2 — 1 times respectively; the corresponding parameter values are 
vy, = (m7)? Py — Ry and », = (m7)? Q:. If we multiply (18) and (19) 
by w(y) and — v(2) respectively, add and set u(x, y) = v(x) w(y) we 
obtain the differential equation 


(20) (Pi uz)e — (Qiuy)y + [Ri + (1 — 1) Ju =0 
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corresponding to (14) and for which explicit solutions sin n; rr sin nz ry 
may be written down. The parameter value \ of (20) is then the difference 
vi — %) = (m7)? Py — Ri — (ne 7)? Q: between the parameter values for 
the equations (18), (19), and is characterized by the number of oscillations 
my — 1, m2 —1 of their solutions. Using the notation n,n, they may be 


ordered as follows: 11, X21, Ar2, Asi, Azz, Auz, ***. Corresponding to the 


“e) 
first we have a solution sin rz sin zy of (20) which does not oscillate, to the 
second sin 272z sin ry which oscillates once as a function of 2, ete. 


The solutions of the equation 
(21) k[ Po tre — Qe Uy + (Re +r)ul = 0 


corresponding to (15) are again of the form sin 7x sin m2 ry but now the 
characteristic number Ayn, is (m1 7)? P2 — Re — (m2 7)? Qe. 

For any given 71, m2 and an arbitrarily chosen e > 0, we may choose m, u 
such that the characteristic numbers of (16) and (17) differ from those of 
(18) and (19) by less than €/2 when m = m,u =u. Hence the characteristics 
Ann, Of (14) and (15) differ from those of (20) and (21) respectively by less 
than e. There is then a parameter value of (A) which lies between 


w(ni Pi —73Q:)-Rit+e and m (nj Po — ni Qo) — Re — €. 


The parameter value Xn n, of the general equation (A) shall be considered as that 
which is obtained by variation when the coefficients of (14) or (15) are varied 
continuously in order to assume those of (A). 

The uniqueness of this determination of the \’s may be shown by the 
following argument. Let p, q, r be chosen as monotone functions of ¢ in the 
interval — 1 =¢ = 1 so that with increasing t the constant coefficients P (t) , 
Q(t), R(t) vary monotonely from P, = P(1), Qe = Q(1), Re = R(1) 
to Pi; = P(—1), Qi. =Q(-—1), Ri = R(—1): the parameter \,,,, (t) 
will increase from the value Ann, (te) = 7 (ni P2 — nz Qo) — Ry to the value 
Anne (ti) = m2 (ni Pi — n3 Qi) — Ri. If the coefficients remained constants 
throughout, for each value of ¢ the solution would be the product of a function 
of x alone oscillating n,; — 1 times by a function of y alone oscillating n. — 1 
times. 

The value of the parameter will decrease if the coefficients vary monotonely 
in any fashion from kP2, kQ2, kR, tokPi, kQ:, kR, even if p(t), q(t), 
r(t) do not remain constant throughout the region for intermediate values 
of the parameter ¢. If the values of the coefficients differ at most only in- 
finitesimally from constants for any value of ¢, the value of the corresponding 
Xn n, Will differ by at most an infinitesimal from the value X,,n,(¢). By the 
introduction of another parameter s to take care of the variation of the coef- 
ficients from the original values P(t), Q(t), R(t), we can have Ann, (t, s) 
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a decreasing function of ¢ for all values of s and a continuous function of s for 
all values of t. But the equation (A) may be regarded as obtained in this 
manner by a monotone variation of the coefficients p(a, y, t, 3), g(a, y, 
t, s), r(x, y,t, 8). Since each of the coefficients p(1,s), g(1, 8), 
r(1,s), p(—1, s), gq(—1, s), r(—1, 8) is a constant for all values 
of s and since Ann, (1, 8), Ann, (— 1, 8) are uniquely determined by mn; and nz 
it follows that Ann, (¢, #) is completely determined in this manner. 

TuHEeorEM I. For each value of m, p there is a characteristic number nn 
of the equation (A) determined by the oscillation numbers n,, nz of the solution 
of (14) or (15). For fixed ny, ne the totality of Ann, lies within a finite interval, 
upper and lower bounds for which may be readily determined. 

Each set of An n,’s will have at least one limit number, and we thus can 
pick out a denumerable set of limit values for the characteristic numbers. 
The bounds of each characteristic number so obtained may be readily calcu- 
lated and definite statements made in regard to the limit points of the set 
for the elliptic and hyperbolic cases. For the former (p > 0, q <0), there 
is a limit point at positive infinity only; for the latter (p > 0, q > 0), there are 
limit points at both positive and negative infinity beside possibilities of more in 


the finite region. 
4. SOME AUXILIARY THEOREMS 


For subsequent use we shall now derive some miscellaneous auxiliary 
theorems relative to maxima and minima of quadratic forms and integrals 
and to the limits of linear forms. 

TueoreM II. The value of the integral of x(du/dx)* dx along the straight 
line joining the points P,(0, M), P2(1, N) CN being less than M ) is less than 
that along any other curve lying above the line and joining these points. 


by 





y 





B (0, M) 











x x 


P, (1, N) P, 








Fic. 1 


To prove this we shall first show that the integral taken along two straight 
lines Q; Qe, Qe Qs is greater than along Q, Q; (Fig. 1). Without loss of gen- 
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erality it may be assumed that the ordinate of Q; is 0. We wish to prove 
that 


f :(= z ay ) ae < i (= 3) dx + f 2(% r Qe ) ae 
2 a2 > J Qe + a3 _. 


a B — Br)? + F 


a3— a Ge = Gi a3 — Ae 


or 


2 


, 


which may be written 


Bi (a3 — ae) Bs (a3 — a) 
Br(as — )™ Bi (a3 — a2) 


But if ag >a, >0 and 0 <f <1 we know that fa; + a2/f > a.+ a. 
Hence our inequality is satisfied when 8i/B82 < (a3 — a1)/(a@3 — a2), which 
certainly holds when Q, lies above the line Q: Q3. 

Since the integrand is never negative the integral along the dotted line 
R, R. whose slope is zero is less than that along the curved line R; R:. In 
the discussion it is then necessary to use only monotone comparison curves. 
But such a function can be approximated as close as we wish by a monotone 
function composed of broken lines. By repeated use of the argument made 


a2 > ae + ay. 


above, our theorem is established. 

Analogous to the division of the plane into rectangles as outlined in § 1, 
one can consider an interval 0, 1 in one dimension as being divided into m 
parts and at the points of division 7/m (i = 0,1, ---,m) variables 2; as 
defined. Within the intervals the function » may be considered as made up 
of straight lines joining the points (7/m,v;). For convenience sake the 
intervals will be taken equal in length. 

TueoreM III. Jf within a subinterval ¢/m, (o + 0’)/m the oscillation of 
the v;’s is M, then m 922" (v; — 14-1)? > mM?/o’. 

For, we have m D922" (v; — v;-1)? = m/o’ [ Det" | 

¢+2’|v; — v:-1| is at least equal to the oscillation, the theorem is established. 

CoroLtuary. When m DX; (0; — v1)? remains limited with increasing m, 
it is impossible that the v’s have a limiting function V (x) which is discontinuous. 

For, at a point of discontinuity of V let M be the oscillation. Given e > 0, 
it is possible to find an m such that in a subinterval of length o’/m < € the 
oscillation of the approximating function is greater than M/2. Then 
m >; (v; — v1) > M?/4e which can be made great at pleasure. This 
being a contradiction, V (2) must be continuous. 

Returning now to the discussion of the two dimensional difference equa- 
tions (A), we shall for the sake of convenience set m = yw and denote by A» 
the maximum of the normalized solution for a given An», and a given m. 


119 . 
v; — V;-1| ? and since 
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Let A» be taken at a point (o, tr) [u(o,7) = An], and denote by (o +0’, 
7 + 7’) the point furthest from (¢, 7) such that there is a chain of points 
connecting them for which u;; > An/2. 

TueoreM IV. If A, is the maximum* of those values u;; for which 


1 9 
22) —- 2... = | 
m? »» 9 : 
we have the formula 
. 
1 ; At o” 
(23) <> m? (ui; — wij)? + m? (ui; — Uizr)?= a 
mS 32* m 


In proving this it may without loss of generality be assumed that o’ = 7’ 


= 0. From (22) it follows that the total number of points where u;; > Am/4 
is less than 16m?/A?. Let us select a rectangle with opposite vertices (o, 7) 
and (¢ + 0’, +r + 32m?/A?® o’) having double this number of points in it. 
On at least half the lines parallel to the y-axis in the rectangle, there is then a 
point where u;; = A,/4 and on every such line a point where u;; = A,/2. 
From Theorem III we have 

1 , enit@es{aYRe _ te 

me me (ig — Misa)? ES 2m > ( 4 ) 32m — 32? m?? 


nj 


~ 


from which (23) follows at once. 
Were then 
a 
lim =o, 
m=a =M 
the value of (23) would increase without limit as A, increases. We can 
investigate the conditions under which 


. Ane’ 
liim— +o, 
am a 


2 


Since in this case o’/m is an infinitesimal of order at least as great as 1/A;,, 
and since the same result would be obtained if we used fractional parts of A» 
other than A,,/2, A»/4, it follows that the point where w;; increases indefinitely 
is isolated. 

~ * A similar theorem can be obtained for the minimum by replacing ui; by — uijz. 


t Were there not at least one such point on one of the four directions from (¢, 7) such 
that ui; > Am / 4 we should have two particular terms, viz., 


1 
= [m? (uer — Ue-ir )? + m? (Uer — Uor-1)?] 


of the sum (23) greater than }A2 but since from (22) A2, < m?, it follows that this is greater 
than A‘ / 2m?; this relation satisfies all the purposes of the theorem. 
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TuHEeorEM V. [f at a point P a function u defined over a unit square R takes 
on a maximum M and if 
f wdR =1, 
R 


foe + uz)dR = M‘o4,/32°, 
R 


then 


where oy is the radius of the greatest circle such that there is a point Q on the cir- 
cumference which may be joined to P by a curve along which u = M/2. 

The proof of this theorem is analogous to that of Theorem IV. 

Let us now deduce some properties of solutions of the elliptic partial dif- 


ferential equation which for the sake of simplicity may be taken in the form 
(24) Ure + Uy + (r+rAk)u = 0. 


If this is multiplied by u and integrated over a field S bounded by a curve C 
there results the formula 


f G+ akyatds = f (at +u3)ds + fuseas, 
s s c On 


where 0u/dn is the normal derivative directed inward; and hence 


u 


(25) max (r+ a) fw dS = f(t + uz) dS + fusas, 
8 Ss Cc 


Using the notation of Theorem V, let K denote a circular conical surface 
having for its axis the z-axis and cutting the plane z = M/2 in a circle of 
radius o,. The surface u(x, y) lies outside the cone at the vertex and 
inside or on the surface at the section z = M/2. Denoting by C, the pro- 
jection on the xy-plane of the intersection of K and the surface u(x, y) and by S 
that part of the region R in the xy-plane bounded by C, then 


2 1 4 9 Y 
(26) f (uz + u3)dS > ra weds. 


For, on passing to polar coérdinates we have 


(27) foe + uZ)dS = [(4 + “) pdpdé = fou pdpdée. 
JS JSS Ss 


From Theorem II we see that this last integral is greater than the value 
computed as if w were on the conical surface, in which case u, would be 
— M/2cy. Hence from (27) 


Me 
(28) f cad + gyas >a pi, (0)d8, 
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where pg, is the value of p on the curve C,. But since uw S M it follows that 


MM? F 
(29) fiw dS = at { as = Mt? | pdpds => f po, (0) dé. 
Ss s 8 “ 0 


Formula (26) follows now from (28) and (29). 
TueorEM VI. The value of the integral 


} (uz + u?)dR 
JR 


for a normalized solution u(x, y) of the elliptic differential equation (24) having 


a maximum M is greater than 
M? 
64? max (r + Ak)° 


For, on the curve (; the derivative du/dn is positive and hence the line 
integral in (25) is positive. Using formulas (25) and (26) we get 


f (uz + ui) dS < max (r + dk ) fw dS 


(30) 
< 4o%, max (r + Ak) f (uz + us) dS, 


and hence 


] 


dou > max (r+Ak)° 


From this formula and Theorem V we have 


2 M4 Om M4 
(31) f ce + u,)dR > > 392 > 642 ais he 


Coro.iary I. The value of the integral f (uz + u5)dR is greater than 
R 
M 


64° max (r + ne) | u2dR 
R 


This may be seen by following through the steps of the reasoning when the 
hypothesis that the solutions are normalized is omitted. 


Coro.tiary II. The maximum value of a normalized solution of the equation 
(24) is less than 8 Vmax (r + dk). 

This follows from (31) and (25). 

Let us next consider some properties of the hyperbolic equation and its 
approximating difference equations. For the sake of simplicity of notation 
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the equations will be taken in the canonical form 


Ury t+ (r+rAk)u =0 
and 
mp (Uj; — Uij—r — Wing + Ui-rj-1) + (rij + Maj) ui = 0, 
the solutions being supposed as before to vanish on the boundary of the 
region. The characteristics are now parallel to the axes and by a simple 
summation of the difference equations we get the formula 


(32) [u(A) +u(C) —u(B) —u(D)] = FD (rs + Mis) ui, 
vu 


where A, B, C, D are the vertices of a rectangle with characteristics for 
sides, the summation extending over this rectangle. 
TueorEeM VII. [If for limited u;; the single sum* 


] , 
= > m? (ui; — Ui-r;)? 


is unlimited for some section parallel to the x-axes it is unlimited for all sections 
and the double sum 


1 ; ae 
— > m? (ui; — Ur; )* 


my “FG 
is also infinite. 
For, if we square the special case of the formula (32) 


1 t 
M(Ust — Use) = M(Uye — Use) — > = (165 + Ake; ) Ue;, 
j=t 


where (s, ¢) is any lattice point and (s, 7) a point such that ¢ > 7, and then 
divide by m and sum over s, we obtain the relation 


Zz m (Unt — Uw)? — p m (Uj, — Ui, )” 
lop 
- adi (tr, k) + D0 f(r, kb) (tee — Wie), 


where the totality of the functions 
; l< 
fir, k) == ff > (te; + Mis; ) Us; 
j=r 


is limited, each being the average of a set of limited functions. The first 
sum on the right being the average of the fi is then limited. If the second 
sum is limited, the sums }; m(uie — wat)? and Di; m(u, — Uia,)® in- 
crease indefinitely together, their difference being finite. If the second sum 
is unlimited, its order of infinity is not more than half that of the sums on the 


ees es re 
* Corresponding to the single integral f u2 dz. 











506 R. G. D. RICHARDSON [October 


left. For, we have 


Do fi(r, hk) (ie — Ui-te) | < max |f;| Ou — vial, 


Fle = tel Sym Ch = te? 
1 t 


Hence the argument holds also for this case. And further if we sum over 
all values of ¢ and multiply by 1/u we are again taking an average and the 
double sum 
=r m? (uUi;— Ui-1;)* 
approaches infinity. 
Similar theorems may be readily established for the sums 


‘ l P 
> (ui; — Uj)? and — > we? (uj; — Uij-1)*. 
j 


mp “F 


In discussing the effect on the solutions of the equations (A) of increasing 
m, pw, let us consider the hyperbolic case (which presents the greater difficulty) 
and for simplicity of notation set p;; = 4;; = 1 and m =y. Adding the 
equations for wj41;, Ui-1j, Uij41, Uij-1 We have 


mM? (Uit2; + Ui-ej — Uij4e — Uij-2) +4 (ry + Maj) Ui; 
(a) = 4( ri; + hij) wig — (ring + Meigs) Wig — (ring + Mea; ) win; 
— (rijga + Mijas) Wiig — (ria + Meij-1) wij, 
while the set of equations for a value of m half the original is 
(8) mM? (Uso; + Ui-2j — Uij4e2 — Uij-2) +4 (ri; + Miz )uy = 0. 


For each value of m the equations (A) will have a set of characteristic num- 
bers X. Let us choose a X for some m and for larger values of m pick out 
(by the method given in § 3) \’s corresponding to this first one, the totality 
of such )’s being limited. On the assumption that r and k and their first 
difference quotients m(r;; — ri;), m(rij — Tij-1), ete., are limited and 
that the values of u;; and the first difference quotients m(u;; — ui-1;), 
m(Uu;; — Uj.) are limited, there will then be an upper bound M for the 
absolute value of these functions and their difference quotients. 

The left-hand members of equations (a), (8) are identical in form and we 
shall now show that the right-hand side of (a) is in all cases an infinitesimal 
depending on 1/m. For, we note that the expression rj; uj; — Ti¢1; Uizaz is in 
absolute value not greater than 

2M? 


! | 
| rigs, = ri3| | wij41| + |Tij| | Uij-1 — U;| < nm’ 
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We have also 


| | 2 2M? 
Vij Uig — Ti-1j Ui-1j < “m ? Viz Ugg — Vijea Uij+i < on ? 
2M? 
55 Ugg — Piz, Usij—s <-—;; 
i J J J 1) m 


and similar relations for the four corresponding terms ini. Hence the right- 
hand side of (a) is less in absolute value than 8M?/m [1 + upper bound of \| | 
and is then an infinitesimal of order at least as high as 1/m. 

As is immediately evident, the same discussion would show that if m were 
multiplied by any power of 2 and the corresponding equations were added 
together in the proper manner, there would result a set of equations each 
differing from (8) only by terms on the right which taken together form an 
infinitesimal depending on 1/m. Such a state of affairs would lead us to 
suspect that we can for the various values of m pick out corresponding solu- 
tions u;; which are very close to one another and have one and only one limit 
function as m increases indefinitely. 


5. THe ELLIPTIC EQUATIONS 


In this section we shall consider the elliptic case of the difference equations 
(A) (p>0,q <0). It will be found possible to trace for various values of 
m, » the relations between the corresponding values of \ and w,; for these 
equations. When m or yp is doubled, the number of equations and the number 
of characteristic values and solutions is approximately doubled. To be more 
exact, it is changed from (m—1)(uh—1) to (2m—1)(u—1) or 
(m—1)(2u—1). Some of the new characteristic numbers are greater 
and some are less than any of the original set, as one may see by studying 
the minimum problem. There are, however, (m — 1)(u— 1) of the new 
set to each one of which there is related one of the old set in that the latter 
is a more or less rough approximation to it and the corresponding members 
of the set of u;;’s are similarly related. If m and yp are further increased, one 
anay trace the new X’s and u;;’s which are related to those already selected. 
And these approach as limit the parameter values and solutions of the elliptic 
partial differential equation 


(33) (puz)2 — (quy)y t+ (r+rAk)u =0, p>o, q <0. 


The details of this discussion we now proceed to give. 

We have already shown in § 3 that for fixed n;, nz the set of parameter 
values is limited for all m,u. From this it is possible to show that the various 
sets of u;;’s corresponding to the d,,n,'s are limited. With this end in view 
let us prove the following 
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TueoreM VIII. If the maximum of the normalized solutions of (A) increases 
without limit with m, pw, the sum 


— Fm? (uj; - 1)” + pw (U5; ‘a Uj; 1)? 


increases without limit. 

To prove this let us make use of the results of Theorem VII and its corollaries 
by passing a surface consisting of portions of planes through the points (7/m, 
i/u, uj). Let us agree to pass planes through the sets of three points cor- 
responding to (i,j), (@+1,j), (4,79 +1), and others through those 
corresponding to (i+ 1,7), (4,j7 +1), (@+1,j7+ 1) and compare the 


integrals 
J wdhk, | (ui + up) dR 
a Ja 


for the total surface found in this way with the sums 


1 9 | 9 9 € 9 
—> Ui; , —> m* (uj; — Uj-1j)? + we? (aay — Uj YP 
mu 


mp FG 
respectively. 


As a first result let us show that the integral 


| (v2 + ui) dR 
~“/R 


is equivalent to 


I ; ' 

(34) —> m* (uj; — Uj-1;)* + we? (Ui; — Uyj-1)*. 

mu S 
If we consider a triangle on the surface defined by the points (7, 7), (¢ — 1,7), 
(i,j — 1) the slope in the x-direction is m(u;; — uj1;) and in the y-direc- 
tion is w(u;; — uj-1). The integrand being a constant and the area being 
1/2mp, the integral has a value 
l 


(35) i m* (Ui; — Uje1j)? + we? (ij — Ujj-1)*).- 


For the triangle determined by (7,7), (-— 1,7), (@-—1,j+1), the 
integral has a value 
l 


(36) — 
; 2mp 


[ me? (1;; — U1; )? + Me (Ui aj4a — Ui1;)*], 


and for the triangle determined by (7,7), (4,7 —1), (@+1,j7—1), the 
integral has a value 


ia l . ‘ . . 
(37) an m* (Uti; — Uiz)” + Mr (Uy; — Ujj-1)*]. 
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By adding (35), the first part of (36) and the last part of (37) we get a term 
of (34) and similar considerations for the other triangles establish this first 
result. 

As a second result let us show that the sum 


is greater than the integral 


where u is the surface consisting of planes. The integral of u? over one of the 
triangles is not greater than 


a max (uj; )° 
nu 


which is less than 
. Per FS | 2 
= ( | Wag] | tiga | H | ise ) 


2mp 


Summing over all the triangles we see that this is less than 


9 
— Fut. 


Mp 
Irom these two results and Theorem VI, Corollary I, we have the formula 
I 9 \e 9 9 M 

— D> m? (uj; — Ui—-1;)® + pw? (ui; — Ui;-1)? > sad? 
, 192? max (r + A‘) a } # u, 
‘mG * 

where M is the maximum of u,;;.. Theorem VIII is then established. 

In formula (8) which when written in full is 


— So [ m? pi; (aig — Uy1;)® — we? qi; (Mj — Uij-1)?] 
a mp G 
(35S) 


| . 
=—) (r;; + ki; ) ui, 
my »» J J ) tj) 
let us consider normalized values of u;;, (20:; ki; ui; = mu), and note that 
l 


| nyne 


which attention is fixed the value of \ is limited, the right-hand member is 


since max |r;;|/k;; is finite and, since for any particular suite of \,,,, on 
limited for all m, w. However, since p;; and — q;; are greater than some 
positive constant, it follows from Theorem VIII that if max w;; increases 
without limit the left-hand member of (38) increases without limit. This 
would establish a contradiction and hence the 

TurorEeM IX. Those solutions of the elliptic case of equations (A) which 


correspond to a particular set of Ann, are limited for all m, pw. 


nye 
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In proceeding to prove that the first difference quotient m(u;; — wj-1;) 
is limited, we shall confine* ourselves to the case where p and q are constants 
and denote them by z and p respectively. It will naturally be necessary to 
put some restrictions on the first difference quotients of the coefficients of the 
equations. Let us assume that the first derivatives of r and k exist and are 
limited. On differentiation of the equation (24) with regard to 2 one obtains 


(39) Wurrr — PUryy + (ru + Aku), = 0. 


We note that u = 0 on the whole boundary of the square, uz = 0 and u,, = 0 
on the sides y = 0, y = 1, and u, = Oand u,, = 0 on the other sides. Hence 
for the special case p = 7, q = p it follows from (24) that uz, = Oand u,, = 0 
on the whole contour. If we multiply (39) by uw, and integrate over the 


square we note from the resulting relation 


(40) f (me, — pu*,)dR = fi (r+rk)u2 + (r+dk), u,uldk, 
R R 





that if wu, is infinite the order of infinity of the left-hand integrand is higher 
by two than the right, which would lead us to suspect a contradiction. 

With this process for differential equations as a guide let us set up by 
subtraction of the equations for i, 7 and i — 1, 7 in (A) a new set of difference 
equations corresponding to (39). For the case p(z,y) = 7, q(2,y) =p, 
after multiplication by m these become 

m3 wl (Wiga; — 2uij + Uns) — (Ui; — 2upaj + U2; )] 
— mp? p[ (Uijga — Qui; + Wij) — (Wiaj — 2uraj + Uinp-1)] 
+ m[ (rig — Taj) +AChig — binj) ) win; 
+ m(ris + Ai;) (uiz — Uns) = 0. 


On multiplication by u;; — ui+;, summation over 7, 7, and division by yu 





(41) 


we get 
m m 
d= [mm (ig; — Wz) — mui; — Urs) P - ™mPlm (Uijz1 — Uiz) 
u 
9 _- @ ‘ 
(42) — m(U-1j5 — Ui-1y-1) P = r De (rig + Meij) (wig — Ui-15) 
Gj 


m 


+7 [ (ri = T;~1;) + d (ki; _ Kea; ) ] (uy a Uj—1; ) U1; 
tj 


Having already shown that the definite quadratic form on the left of (38) is 
limited for all m, yw, it follows that 


m 


2 (uy = wes) 


* This change is made here solely to avoid complicated notation. 
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is limited. From this we may deduce that since r;; + Ak;; is limited, the 
first sum on the right of (42) is limited and further from the formula 


my } (uj; — U1;)? = [Do lu; — u-a;| FP, 
G 


9] 


or 
m ae 2 
m= d (33 — U1; )° = |= Us — Us—1; | 
that 
] ; oe 
M y 


is limited. The second sum on the right of (42) being less in absolute value 


than 


1 
max m| (rij — rij) + ACh; — kien; ) | max) uy; aN Uij — Uj-1;|, 
j a 


ij 
must then remain limited for all m, uw. 
The equation (39) may be written 


(mUz)ex rani ( puz) yy + (r + Ak) uz == (%, + rk.) u 


and regarded as a non-homogeneous equation of the second order in w;z. 
All that is essential in the reasoning of Theorems VI, VII, VIII may be applied 
to this equation and the corresponding difference equations when wu and 1;; 
are limited. Regarding the difference quotient as unit instead of w;; itself, 
it follows at once from a theorem analogous to VIJI that the left-hand side 
of (42) increases without limit with m(u;; — u;4;). This would lead to 
a contradiction. We have then 

THEOREM X. [If the coefficients of the elliptic case of the difference equa- 
tions (A) are values taken from functions which have continuous first derivatives, 
the first differences m(uj;; — Wins), w(Uiz — Uij-1) which correspond to a 
particular set of Ann, are limited for all m, p. 

By setting up difference equations by subtraction in (41) or (A) it may be 
proved in the same way that the three second difference quotients corresponding 
10 Urz, Ury, Uy, are limited. To indicate briefly the method used in one of 
these cases let us select that corresponding to u,, and confine ourselves to 
derivative notation, remembering that every step may be duplicated by 
difference equations. On differentiating (39) with regard to y we get 


(43) TUryrr — PUryyy + (ru + Aku)z, = 0 


and on multiplication by uz, and integration, the formula* 


f ime, — pu? ]dR = fo + rk) ui, dR + fic +rk)z Uy 
(44) R R R 


+ (r+k)yuct (r+ Ak)zy uj) uz dR. 


* The vanishing of u,z and uy, on the contour entails the vanishing of uzzy and Uzyy also. 
Trans. Am Math. Soc. 33 
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Since the right-hand side of (40) is limited it follows that the left-hand side 
and hence that the integral of u?, is limited. The first integral on the right 
of (44) must then be limited and on the assumption that r,,, kz, are limited 
the second is a fortiori limited.* But were u,, itself unlimited the integral 
on the left could be proved unlimited by an extension of Theorems VI, VII, 
VIII, where u,, (and not w) is regarded as the unit. This would establish a 
contradiction. 

And by an extension of the method, higher difference quotients can be 
shown limited if enough derivatives for r and i are assumed to exist. 

From well-known theorems it follows that since w and its first difference 
quotients remain limited for all m, w, there is at least one limit function u(x, y), 
which is approached uniformly by the solutions of the difference equations. On 
the assumption that the functions p, q, r, k have limited third derivatives, the 
function u also possesses continuous first and second derivatives. On account 
of this continuity the function u(x, y) satisfying (33) for an everywhere. dense 
set of points must throughout the region be a solution of this differential equation. 

From relation (8) it follows that each set of parameter values of the dif- 
ference equations approaches the corresponding parameter value of the 
differential equation uniformly. The X’s are roots of an algebraic equation 
of order (m — 1)(u — 1) and being continuous functions of its coefficients 
are continuous functions of the parameter ¢. Since the u;;’s being roots of 
the linear algebraic equations (A) are continuous functions of the parameter 
t on which the coefficients p;;, 9:;, rij, kiz; and X depend and, since they 
approach the function u(2, y) uniformly, it follows that the limit function 
u(x,y) is a continuous function of t. From this we can deduce that there 
is only one limit for the set of u;;'s. For, in the case of equations (16) in one 
dimension with constant coefficients it may be assumed from the theory of 
ordinary differential equations that there is only one limit for each of the sets 
of v’s and w’s. For the case of the equation (14) in two dimensions with con- 
stant coefficients it follows from this immediately that the u;;’s have one limit 
function only. In considering as in $3 the coefficients of (A) as functions 
of ¢ if there is a point fo for which there is a single solution U,,,, (2, y) of (33) 
while for t) + € there are two solutions Unn,(2,y), Ung (2, y) for every 
positive e€, we have as the limit of (9) the orthogonal relation 


Jf Unig t +€,2,y)Unn,(to+¢,2,y)dR =0, 
R 


which on passing to the limit e = 0 gives 
* The function r + Ak with its derivatives and the functions u, uz, u, being limited, the 


2 x . 2 
only inequality needed is f uy dk = [S. Ury | ar | , 
R 








eo ——— 


\ 
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f BU (to x,y)dR=0. 
R 


This is, however, impossible since the solutions are normalized and the result 
is established. We may now formulate 

TuEoREM XI. There exists in the orthogonal elliptic case (q <0, k > 0) 
a denumerably infinite set of parameter values \, with limiting point at infinity 
only, for each of which there is a solution of (33) vanishing on the boundary. 

A similar fundamental theorem may be.established for the polar case of 
(33), viz.; when r(2,y) = 0 and the region R may be divided into a finite 
number of simple regions in the interior of some of which & is positive and in 
the interior of the remainder of which & is negative. The essential point of 
difference in the treatment of the polar case is in showing that the corresponding 
suite of Xn,n,’s for the various values of m, uw is limited. This may be estab- 
of the 
equations (33) are increased by decreasing the size of the region and retaining 


lished by showing that: (a) the absolute value of all parameters Aj,», 
the condition that « vanish on the boundary; (b) on taking for the smaller 
region a sub-rectangle in which k > 0, the value of each of the positive \j,n,'s 
is still further increased by replacing k by its minimum (> 0) in the sub- 
rectangle and r by its maximum. This gives an upper bound (which may 
be calculated explicitly) for each positive An,n,; (¢) on taking for the smaller 
region a sub-rectangle in which / < 0 the value of each of the negative A,,n,’s 
is still further decreased by replacing & by its maximum ( < 0) in the sub- 
rectangle and r by its maximum. ‘This gives a lower bound (which may be 
calculated explicitly) for each negative X,,n,- 


6. THE HYPERBOLIC EQUATIONS 


Preliminary to an investigation of some properties of the hyperbolic differ- 
ence equations let us study some properties of the hyperbolic differential equa- 
tion 
(45) (pure — (quy)y + (r +Ak)u = 0 (p>0,q>0). 


While an analogous discussion can be given for the more general equation, 
we shall, for the sake of simplicity, confine ourselves to the case where p and q 
are constants. The characteristics will then be sets of straight lines making 
equal angles (y = arctan + Vq/p) with the z-axis. It will be assumed 
that r and k are continuous with as many derivatives as desired and that the 
solutions of (45) are continuous and possess first and second derivatives. 

TueorEeM XII. If the sides of a parallelogram ABCD are characteristics of 
the differential equation then 


(46) 2vVpg[u(A)+u(C) —u(B) —u(D)] = f (r+ ne)uaP, 
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the double integral being taken over the interior of the parallelogram. This formula 
holds identically for all such parallelograms and for all solutions of (45). 











Fic. 2 


For, on integrating the equation (45) over P we get 


[ pu. dy — qu,dx = — fo + rk)udP, 
—~7T e/p 


where I is composed of the sides 4B, BC, CD, DA. To calculate the line 
integral along AB, we resolve it into two parts, one along AB and the other 
normal to it. We have the formulas 

= \ P ds , dy = — |_4 ds , 

pt+q , \ p+q 


| x _ en oe 


——— © Seer en 2 
Uz = “No +4q No +q’ t= “No +q) "No +9 


dx 


and this part of the line integral becomes 


f ue Vpq ds = Vpq{u(B) — u(A))]. 
AB 


We have further 


f pu, dy — qu,dx = Vpq[u(B) —u(C)], 
BC 


f pu,dy — qu,dx = Vpq [u(D)—u(C)], 
“cD 


{ pu, dy — qu, dx 
DA 


Formula (46) is then obtained by addition. 
Let us make one or more of the following hypotheses. 


Vpq{u(D) —u(A)]. 
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Hyporuesis I. The solutions of (45) are normalized; that is, ff bu dR 
R 


= constant, where R is a unit square containing the parallelogram. 


If r and ) are limited the integral (r+Ak)u'dR will be limited. It may 
YR 


be assumed without loss of generality that this last integral is 1. 
Hyportuesis II. The integral { (uz + uZ) dR is limited. 
R 


Solutions which satisfy I cannot be infinite at a finite number of points 
in the region, as is readily seen from formula (46). 

THeorEM XIII. Solutions of the differential equation (45) which satisfy 
Hypotheses I and II are limited. 

For, let u take on a large value M at the point A. It follows from I that 
there will be points in the neighborhood of A where u < M/2. We may 
distinguish (Fig. 2) two cases: (a) u(B) = u(D) = M/2, while within the 
intervals AB and AD, u > M/2; (b) along at least one of the characteristics, 
u > M/2. Incase (a) we have from (46) and Hypothesis I that u(C) <1. 
Let 71 and m2 be lengths of AB, AD respectively, and 4:1 = 72. For two 
corresponding points EK and E’ we have on applying formula (46) to the 
parallelogram ABE’ E, the relation u( EL) —u(k’) > M/2—1. The 
double integral of u? + uj is independent of a rotation of axes and if the 
x-axis is taken in the direction AB we have from a theorem of integrals analo- 
gous to Theorem III the formula 


M 3 
E aE? _ ss% 
f (w+ uz)\dxr= i} uz dx = -— . 
E ~w/E 1 
Hence 
M 3 
a 


° 9 2 ’ M 
f (uz + uy) dydxr = = n2 cos 2y = ( > 1 ) cos 2y, 
p 


and since this integral increases without limit with J/, a contradiction estab- 
lishes the theorem in this case. 

Case (b) may be readily handled by the methods used in proving The- 
orem IV. 

Instead of allowing u to increase indefinitely we can fix its maximum equal 1 
and have a suite of functions u for which Hypotheses I and II become 


Hyportuesis I’ lim f wdx = 0, 
R 
Hypotuesis II’ lim | (uz + u;)dx = 0. 
R 
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From Theorem XII it follows at once that there can be no series of solutions 
of the differential equations having these properties. 

Every solution of the differential equation (45) can be approximated as 
close as we wish in the parallelogram P by a solution of the difference equa- 
tions; by defining a surface made up of portions of planes through the points 
u;; one obtains from equation (A) the formula 


— 1 
(47) 2Vpqtu( A) +u(C) — u(B) —u(D)] = 22 De (nist Mis) Wis bemn » 
y 


where the sum is taken over the lattice points within or on the parallelogram 
and €, is an infinitesimal depending on m, p. 

One can, moreover, get a precise formula analogous to (46) for the difference 
equations when p, q are so chosen that the points of the lattice lie on the 
characteristics, that is, if Vp/q is a rational number. Let us derive such a 
formula for the case p = q = 1. The results in the more general case are 
similar. The difference equations for u;; for the case p = 1, q = 1 are of 
the form 

Uieag Hb UiRrz — Nig — Uj = — 1 /mp (rij + Nhe jj) Ui;. 
Let us sum these over a rectangle whose vertices A (01,71), B(o, 72), 
C(o3, 73), D(o4, 74) are at lattice points and whose sides make angles of 
7/4 with the axes. We note that for every point within the rectangle the 
term u;; occurs four times in the sum: once each in a positive sense for those 
equations which enter for points next to the left and next to the right; once 
each in a negative sense for those next above and next below. This is also 
true for those points on the sides of the rectangles which are not vertices. 
For points lying in the next row below AB and to the left, with the exception 
of the end points (o; — 1, 71), (2, 72 — 1), the term occurs twice in the 
sum, once as positive and once as negative. A careful analysis of the remain- 


ing terms gives the formula 
_— + ay + ( Moss, + Ussti15) _ (= + a 


(48) 
_ (og, + Uo grt )= 1/mp Zz ( rij + Ni; ) Uij- 
ij 


Since under Hypothesis II the difference between the values of u at two 
neighboring points is less than a constant G, we have 

2[u(A) +u(C) —u(B) —u(D)] = 1/mp SS (ris + Mj) Wij + mg » 
U 

lMmy | < 4G 


(49) 





In general if Vp/q is rational we get a formula similar to (48). 
Turorem XIV. When vVp/q is a rational number, while 


] 2 
—> kj; Uz = | F- 


Mp i 
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and 


1 
peri >» m?* (Uj; oo Ui—1;)? + we (U5; — Ujj-1 )? 


my GF 
is limited, it is impossible that max u;; increase indefinitely with m, p. 

This is proved in a manner analogous to Theorem XII. While in the 
case where Vp/q is irrational precise formulas corresponding to (48) are not 
available, there seems to the author no reason to think that the theorem does 
not hold in that case. This theorem is based on formula (48) which depends 
in no way on boundary conditions but is identically true (when Vp/q is rational) 
for all solutions and all points of the lattice. Further, the solutions of the 
difference equations are continuous functions of the parameter ¢ on which 
the coefficients p, qg depend. From the principle of continuity we should 
expect a formula similar to (49) to obtain for all p and q. 

As was seen in § 3, for each m, uw the limited function X,,», (¢) is monotone 
in the interval — 1,1. Each function possesses a derivative and the points 
where — dAnn,/dt = M > 0 form an aggregate whose measure is less than 
c/M, ¢ being a constant independent of M, m, uw. For the infinite set of 
functions Ay,n, corresponding to various values of m, wu there must exist at 
least one limit function and this function will be limited and monotone. It is 
well known that such a function possesses a limited derivative except possibly 
at a set of points of measure zero, 

In formula (13) 


dhe ddl er ae I dri; , AX —dki; , 
dt dt mG mag Sag ™ mp “TZ dt mes my GZ dt Mis 
1 dpi; , ms , d4;; ; 
+S - m= dp (ui — Uj;-1;)° — pe dt (uj; — Ujj-1)° 


we may take (as in § 3) the limited derivatives — dp/dt, dq/dt, dr/dt, dk/dt 


positive and argue that the sum 
eae . m” ( Uiz — Us-1;)* + we ( Usjj — Usj-1 y 


is limited for those values of ¢ for which d\/dt is limited. It is then possible 
to enunciate the following 

TuroreM XV. For all values of t except possibly for a set of measure zero, 
the sum 


1 , 
—— Dm? (us — Urs )® + we? (uis — Wij)? 
my G 

as limited. 


If for a value of t for which d)/dt is finite a formula analogous to (48) is 
valid, we can deduce at once that u;; remains limited for all m, u and that 
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there is a limiting function u(2, y) satisfying the differential equation. If, 
for example, for t = to we have both dd/dt limited and Vp/q rational this would 
be the case. 

From another point of view, any solution of the difference equations is an 
approximation to a solution of the differential equation. With increase of 
m,m@, a sequence of limited solutions of the difference equations which at a 
given lattice point are made to take on a value 1 will have at least one limiting 
function which is a solution of the differential equation at that point. But, 
as was seen above (cf. Hypotheses I’, II’), for any sequence of solutions of the 
differential equation it is impossible that 


lim f n° dx = 0, lim fc +uzy)dx = 0. 
Rk YR 


But this seems to imply that for limited d\/dt, normalized solutions of the 
difference equations like normalized solutions of (45) remain limited. 

It may be proved by a method similar to that used in Theorem VII that 
when the double sum 


ym m2 Pid (uj; — Uji1;)? — yo (uj; — Uj)” 

mp G dt 7” =e “as ” 
remains limited as m, u increase, the single sums 0; m(u;; — u;—1;)* for all j 
and Dj m(u;; — u;;-1)* for all 7 remain limited. It follows from the corollary 
to Theorem III that the limiting function u(2, y) of the sets of solutions of 
the difference equations (A) is a continuous function for all lines parallel to 
the axes which belong to the lattice for any m, w and hence is continuous 
throughout the square. 

If we make the hypothesis (which appears a most reasonable one) that a 
formula analogous to (48) can be obtained for any set of hyperbolic difference 
equations (A), there would be for all values of t (except perhaps for a set of 
measure zero) a bounded limiting continuous function u(2, y) which would 
be a solution of (45). If there were a value of ¢ for which such a function u 
did not exist, there would be coefficients as near to the original as desired 
such that for these the differential equation would have a solution. 


Brown UNIVERSITY 











NOTE SUPPLEMENTARY TO THE PAPER “ON CERTAIN PAIRS 
OF TRANSCENDENTAL FUNCTIONS WHOSE ROOTS 
SEPARATE EACH OTHER” 

MAXIME BOCHER 


In the paper with the above title, which appeared in these Transac- 
tions, vol. 2 (1901), p. 428, I obtained a number of general theorems 
concerning the zeros of functions of the form ¢2 y’ — ¢1 y, where y is a solu- 
tion of a homogeneous linear differential equation of the second order. There 
s a further general theorem of the same character which I should have included 
in that paper if I had discovered it at the time, and which I now give as VII’ 
below, it being a counterpart to VII in the earlier paper. I also give three 
special applications of this general result, numbered VIII’, VIII’, VIII”, 
since they are counterparts of VIII. A part of VIII’ was obtained by Fite 
inthe Annals of Mathematics for June, 1917, (see the closing 
lines of his article) by another method; and it was this special result of Fite’s 
that suggested to me the much more general Theorem VII’. 

What follows should be regarded as an addendum to the former paper, 
to be inserted on page 434 at the end of § 2. All references are to that paper, 
to which the reader should turn for an explanation of the notation and a 
statement of the restrictions placed on the functions. 

VII’. If none of the six functions 


oi v2 — doi, Wix2 — W2X1, x1b2 — x2d1, (hi, bo}, {Wi, We}, {x15 xe} 


vanishes in (1), and if the last of them has the opposite sign from the two pre- 
ceding ones, then no one of the functions ®, VY, X can vanish more than once in 
(I): and no two of them can vanish except perhaps ® and VW. If these last two 
functions do both vanish in (1), the root of ® will be greater or less than the root 
of WY according as the product of the six functions written above is negative or 
positive. 

The first part of this theorem is an immediate consequence of IV, p. 431. 
Let us then assume that ® vanishes at 21, V at x2. These two points can- 
not coincide, since then ¢: 2 — ¢2 ¥1 would vanish at this point. 
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Suppose, first, that 71 < #2. Then 
sgn ®(a.) = sgn ®’(2,), 
and, consequently, by (3), p. 430, 
sgn (x2) = — sgn (pi Yo — $2 Yr) sgn {d1, G2} sgn V(x). 


On the other hand, if we substitute in the identity (6), p. 432, first x = 21, 


then x = x2, we find 


— sgn (x1 $2 — x2 $1) sgn V(x), 


Il 


sgn (di 2 — $2 ¥1) sgn X (21) 
sgn ($1 ¥2 — o2¥1) sgn X (a2) = — sgn (Yi x2 — Yo x1) sgn P( 22). 


Multiplying together the last three equations, and remembering that X 


does not vanish, so that 


sgn X (2,1) = sgn X(a) = +1, 
we find 


sgn ®(22) = — sgn [(d1 v2 — bo 1) (x1 2 — x2 1) (Yi x2 — Ye x1) ] 
sgn {d1, do} sgn ®(2z2). 


If here we cancel the factor sgn ®( 22) = + 1 from both sides, and replace 
— 1 by its value sgn {¥1, Yo} sgn {x1, x2}, we infer from this equation that 
the product of the six functions of VII’ is positive. 

In a precisely similar manner, we show that if 2; > 22, the product of these 
six factors is negative, and this completes the proof of VII’. 

The theorem just proved admits of extensions in two directions. A brief 
mention of these will be sufficient. 

First, the assumption that {¢1, 2}, {¥1, Ye}, {x1, x2} do not vanish may 
be replaced (at least if we assume that ¢1, ¢2, ¥1, ¥2, X1, X2 have continuous 
first derivatives) by the less restrictive conditions that these three functions 
do not change sign and that none of them vanishes at all points of a connected 
subinterval of (I). This extension is effected as in similar cases on p. 436. 

Secondly, we may suppose that ®, instead of being of the form @2 y’ — diy, 
is of the form ¢1 y — $22, where (y, z) is a solution of a pair of homogeneous 
linear differential equations of the first order. Making corresponding changes 
in the definitions of YW and X, we obtain in this way a counterpart to XVI 
in these Transactions, vol. 3 (1902), p. 207, the proof being closely 
parallel to that we have here given. 

We make a first application of VII’ by letting 


7 


@=y, WV=y, X=py’+quy=-y’, 
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so that 


dit2 —dofi = —1, Yixxe—¥ex1i = X1%2 — x2¢1 =P, 


w| ad fi 
{oi, do} = 1, {vi, Yo} = 4, {X1» X2} = e| =( e) + |. 


Theorem VII’ now becomes* 
VIII’. If 
d (p F 
q>0, —E> Pe ~ ft, p nor wre zero, 
dx\ q 
then no one of the functions y, y’, y’”’ can vanish more than once; and no two of 
them can vanish except y and y’. If these two functions do both vanish, the root 
of y is greater or less than the root of y’ according as p is negative or positive. 
As a second application of VII’, we let 
se97, Tey, Go<s 
and thus find 
VIII”. If 
d ({p 
q <0, —i=-i<-—l, Pp nowhere zero, 
dx\ q 
then no one of the functions y, y’, y’’ can vanish more than once; and no two of 
them can vanish except y’ and y"’. If these two functions do both vanish, the 
root of y’ is greater or less than the root of y’’ according as p 1s negative or positive. 
Finally, letting 
G=7, Lvs, Co-s", 
we obtain 
VIIl’”’. If 
d{p 
q <0, —{=-)>-—l1, p nowhere zero, 
dx\ q 
then no one of the functions y, y’, y”’ can vanish more than once; and no two of 
them can vanish except y and y”. If these two functions do both vanish, the 
root of y is greater or less than the root of y'’ according as p is negative or positive. 
HARVARD UNIVERSITY 
July, 1917 
* The latter part of this theorem is the result of Fite above referred to, except that Fite 
requires merely 
d (p 
bs ( e) ~ 
dx \q 
This generalization, however, is not a substantial one, since the equality sign may be removed 


by an arbitrarily small variation of p, and such a variation, if the root of y’ is held fixed, 
produces only a slight change in the root of y. 











A THEOREM FOR SPACE ANALOGOUS TO CESARO’S THEOREM 


FOR PLANE ISOGONAL SYSTEMS* 
PERCEY F. SMITH 


Investigationst in the geometry of ordinary differential equations of the 
first order have yielded many interesting results, especially in the properties 
of the osculating circles of the integral curves of such equations. In particular, 
Scheffers (1. ¢.) has introduced two simple lineal element transformations, 
which have been well described by Kasner{ as “turns”’ and “ slides,” re- 
spectively. A turn converts a lineal element into one having the same point 
and a direction making a given angle with the original direction. Under a 
slide an element is displaced a given distance along the line of the element. 
Applied to the integral curves of 


(1) F(x,y,p) =9, 


a turn of each lineal element satisfying the equation leads to a second equation 
whose integral curves form a system of isogonal trajectories of the integral 
curves of (1). A slide leads to equitangential trajectories. Cesdro§ has 
shown that the osculating circles of the curves of an isogonal system at any 
given point will pass through a second common point. Scheffers (I. c.) has 
demonstrated a companion theorem for an equitangential system of curves, 
namely, the osculating circles of the curves of such a system which touch a 
given line will also touch a second line. 

In approaching questions of this nature for space, it is clear that an under- 

. . . . . . i 
standing must be reached as to the direction the investigation should take.|| 

* Presented to the Society, December 2, 1916. . 

t The following papers may be referred to in this connection: Scheffers, Jsogonalkurven, 
Aquitangentialkurven und complexe Zahlen, Mathematische Annalen, vol. 60 
(1905), p. 491; Gundelfinger, On the geometry of line elements in the plane with reference to 
osculating circles, American Journal of Mathematics, vol. 33 (1910), p. 153. 

t Kasner, The group of turns and slides, etc., American Journal of Mathe- 
maties, vol. 33 (1910), p. 193. 

§ Geometria intrinseca, 1896. 

|| Kasner, for example, has proposed and solved a problem for space, different from that 
discussed here, in his paper Equitangential congruences of curves in space, Rendicontidel 
Circolo Matematicodi Palermo, vol. 35 (1913), p. 283. 
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To consider partial differential equations of the first order and their integral 
surfaces is a first thought, but to consider osculating spheres of such surfaces 
as analogous to the osculating circles of plane curves is wide of the mark, 
since at a general point on a surface there is no osculating sphere. Further- 
more, a turn, or a slide, of a lineal element is definite in so far as the point, 
or the line, of the element is unique. It is natural to define a turn of a surface 
element as a rotation of the element through a given angle about any axis through 
the point of the element and lying on its plane. Likewise by a slide we may 
understand a translation of the element in its plane a given distance in any 
direction. But in the first instance, the axis of the turn, and in the second, 
the direction of the slide, are not definite. A study of the geometry of partial 
differential equations of the first order, to which this paper is devoted, clears 
up the vagueness referred to above, and leads to theorems which appear to 
be as fundamental in connection with such differential equations as those 
cited for ordinary differential equations. 


1. PRELIMINARY FORMULAS 


Consider a surface element (x, y, z, p, q), and a consecutive element 
(2a+dx, y+dy, z+dz, p+dp, ¢+ dq) united with the first, that is, 
such that 
q (1) dz = pdx + qdy, 


the ratios of the infinitesimals* satisfying 


> dx dy dz _dp_ d& 
(2) a Bp pat+qB 6 «€° 


The pair of united elements will establish ©! curvature elements (2x, y, z, 


P,4q,v,8,t), where the last three codrdinates satisfy 
(3) ; 6=ra+s8, €=sa+ ip. 


For any surface to which the pair of united elements belong, the curvature 
of normal sections at (2, y, z) containing the direction 6x : dy : 62 is 


| 1 r 6x? + 2s bx by + t by’ 
(4) R on 9 9 22 - 
t vi+ p+ ¢ 6s 
Eliminating r and s by using (3), this becomes 
(5) 1 1 B6 ba” + ae by? — s(Bbx — aby)? 
5 ra SS cre ———— a ‘ 
R V1 -+ p" + g° ap 6s 


* Cf. a paper by the author On osculating element-bands associated with loci of surface elements, 
these Transactions, vol. 11 (1910), p. 302. 
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The pair of surface elements under discussion present two directions (usually 
distinct) which possess obvious geometrical interest and simplicity. These 
directions establish a line-pair lying in the plane FE of the first element 
(x, ¥, 2, p,q) and passing through its point P(x, y,2z). The first direc- 
tion is that of the line joining the points of the united elements. Here 


(6) bx: by: 62::a:B: pat QB. 


The second direction is that of the line of intersection of the planes of the 
elements. For this line 


(7) bx : dy: 62 =e: —6: pe — gd. 
For the first direction, equation (5) becomes 


(8) l 1 ab+ Be 


Ro V+ p+ ee +h + (pa + 9B) 

The disappearance of s indicates that the normal curvature of every surface 
containing the pair of elements is the same in the direction in which these 
elements are united. The normal sections of these surfaces in this direction 
have therefore a common osculating circle. Hence we may associate with a 
pair of united surface elements a first osculating circle whose radius is given by 
(8) and which has the property described in the preceding sentences, namely, 
it is the osculating circle of the normal section of any surface containing the element- 
pair taken in the direction in which the elements are united. 

teserving for a later section (§ 3) a similar discussion for the direction (7), 
pass now to the application of the preceding exposition to the partial differ- 
ential equation (assumed not linear) 


(9) F(r,y,2,p,q) =9, 


in which, as usual, p = 02/dx, g = 0z/dy. Let S be an integral surface of 
this differential equation, P(a2, y,z) a point on S, and C the characteristic 
on S through P. Infinitely many integral surfaces of (9) contain the char- 
acteristic C and all touch S along C. On all of these surfaces the direction 
conjugate to that of the characteristic at P is common.* That is, in each 
surface element satisfying (9) there are two significant directions, namely, 
(1) the direction of the characteristic determined by this element, and (2), the 
common conjugate direction above described. If we consider a pair of united 
elements satisfying (9) whose points are on the characteristic determined 
by the first element, these two directions are obviously precisely those dis- 


* Carathéodory, Zur geometrischen Deutung der Charakteristiken einer partiellen Differential- 
gleichung erster Ordnung, Mathematische Annalen, vol. 95 (1904), p. 377. An 
excellent exposition of the geometry involved here is given by Goursat in his Legons sur l’in- 
tégration des équations aux derivées partielles du premier ordre (Paris, 1891), pp. 181-88. 
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cussed above (equations (6) and (7)). For this pair of elements, the infini- 
tesimals dx, dy, etc., satisfy the well-known differential equations of the 
characteristics (or characteristic bands) 

dx dy dz _ dp - dq 

F, F, pFp+qF, —¥F:—pFk: —F,—-qF:’ 


Pp 





(10) 


where subscripts indicate partial differentiation with respect to the variables 
expressed. 

To obtain a first osculating circle for any surface element satisfying (9), 
we may use formula (8) in which a, 8, etc., of equations (2) are to be replaced 
by the corresponding denominators in (10). The property possessed by this 
circle is clear. The normal section of any integral surface containing this surface 
element taken in the direction of the characteristic will be osculated by this circle. 
In the following discussion the phrase “ first osculating circle for the surface 
element ” will be employed to designate this circle. 

To derive a convenient expression for the radius of the first osculating 
circle, it is desirable to introduce ‘“ homogeneous element coédrdinates ” by 
setting 

Pi P2 


11) = = mw EE, 
( ) P P3 : Ps 


The given differential equation (9) now becomes 
(12) F(2,y,2,p,q) = G(x, y,2, pi, po, ps) = 9, 


where G is homogeneous and of degree zero in pi, p2, ps. Hence, if we set 
G; = 0G /dp; (a = P 2; 3) ’ 


(13) pi Gi + po Go + p3 G3 = 0. 


Furthermore, F, = — p3Gi, Fy = — psGe, pF» + qFq = — psGs. Then 
we easily obtain from (8), (2), and (10) the symmetrical formula* for the 
curvature of the first osculating circle 


dia 1 GiG.+GoGy + 63G, 
RR, o(@G7+634+ G4) ’ 
where w? = pi + p; + pi. 
This convenient formula being secured we proceed in the following section 


*It is not without interest to derive an expression for the curvature of the osculating 
circle at (x, y) of the integral curve of the ordinary differential equation F(z, y, p) = 
determined by a lineal element (x, y, p) satisfying this equation. In fact, if we place 
p = — pi/ pz, and proceed from the usual expression for curvature, it is easy to show that 
the resulting formula is precisely (14) written for two dimensions instead of three. In other 
words, the first osculating circle here introduced is a literal generalization to three dimensions of 
the osculating circle of the plane. 
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to apply to each of the surface elements of equation (12) the transformation 
designated a turn. 

In the analytical developments which follow it will be necessary to refer 
to a result which may be assumed from differential geometry. 

If (a1, a2, a3), (li, le, ls), amd (Ax, Ae, As) are the direction cosines of 
three lines mutually at right angles, and if such directions are chosen on these 
lines that | 


a, A: 3) 


(14’) l; ly l; = aa l 
Ar Ae Az 
(rather than the alternative value — 1), then, indicating by primes differ- 


entiation with respect to a common parameter, relations will hold of which 
the following is the type: 


jf ly l. l; ls ls | 
(15) a; =I;\h Ig I] —dsiAr Ye Asl- 
ids de Ke Ai Ae ASI 


2. THEOREMS FOR AN ISOGONAL SYSTEM 


Turn now to the differential equation ; 
(16) G(x, Y,2, Pi, Po, Ps) = 0, 


and consider one of its surface elements (2, y,2, 1, P2, P3). We give to 
the normal () of this element a definite direction by choosing its direction 
cosines \;, Ae, Az such that 


(17) Ay = pi/o, Ae = pPo/w, As = P;/w, wo=+ Vpi + pi t+ pi. 


To the tangent (/) of the characteristic through (2, y, z) in this element 
we assign as positive direction that whose direction cosines are (see (10)) 


(8) L=G6,/6, h=G6./6, 4=6,/6, G=4+V@4+@6+4. 


For F,: Faq: pFp + qFq = Gi: Ge: Gs. 

We next establish an axis (a) through (2, y, z) in the plane of the element 
at right angles to (/) with direction cosines (a ;, a2, a3) and such that (14’) 
holds. A rotation of the element about the axis (a) is now entirely definite, 
and we may set down the 

DEFINITION. A turn of a surface element satisfying a partial differential 
equation of the first order is a rotation through a given angle about the axis (a). 

The equations for a turn 7’, through an angle 6 are obviously, 


(19) XN; = A; cos 6 + 1; sin 6 (i =1,2,3), 


in which the direction cosines of the normal to the element in its new position 








Min 


CPS PRS wick. 











as 





PRET ci. See 
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are (A;, 43, A,). The turned elements will also satisfy a partial differential 
equation of the first order. To obtain this equation we replace \; and /; in 
(19) by their values as defined in (17) and (18), and from the three equations 
thus obtained and (16) eliminate p,, pe, and p3;. Then placing 


(20) N=p/o', wo =+VR +r +r (i= 1, 2,3), 
the result of the elimination may be indicated by 
(21) G' (2,9,25 Dis Per Ps) = 0, 


in which the parameter @ is implicitly involved. ‘Thus equation (21) repre- 
sents a system of differential equations, which will be called an zsogonal 
system. Justification for this terminology will appear presently. 

The capital fact in connection with the turn 7’, given by (19) is this. For 
a fixed point (x, y, z) the turns form a group. In other words, the direction of 
the characteristic in a turned element is perpendicular to the axis (a). To 
establish the truth of this statement the following considerations may suffice. 
The planes of the surface elements of the differential equation (16) with a 
common point P envelop a cone,*—the elementary cone of Monge. The 
elements of this cone are the tangents of the characteristics through P. A 
turn converts the tangent planes of this cone into those of the elementary 
cone of the transformed equation. The elements of the original cone also 
“turn” into those of the transformed cone. This simple geometric fact is 
not common knowledge but may readily be established. 

An important consequence of the preceding result is the following. If 
(1;, U3, I) are the direction cosines of the characteristic of G’ = 0 in the 
turned element (2, y, 2, Pi, Ps, Py), then 
(22) I; = — X; sin 6 + 1; cos 0 (i =1,2,3), 


since this characteristic is at right angles to the axis. Here /; = G;/G’, and 
G” = DiiG?. 

It will now be possible to prove for the isogonal system (21) a theorem 
analogous to that of Cesdro for plane isogonal systems. The preceding 
presentation should have made clear the following situation. 

(1) Turning an element of the original partial differential equation G = 0 
gives rise in the isogonal system G’ = 0 to a pencil of elements whose axis is 
the axis of the turn. 

(II) The characteristics in this pencil of elements are at right angles to the 
axis. 

(III) The first osculating circles for the elements of this pencil are coplanar, 

* For simplicity of statement, the equation is assumed not linear. The conclusion, how- 


ever, is true also for this case. 
Trans. Am, Math, Soc, 34 
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their common plane being at right angles to the axis. Moreover, these circles 
have by definition one common point,—the common point of the pencil of 
elements. 

To show that these osculating circles have a second common point is the 
next step. For this purpose, apply formula (14) to the curvature of the first 
osculating circle in the isogonal system (21). This leads to 
(23) ca i, G, + Gz G, + Gs G! 
RK (G6 +G6;, + Gye) 


Making the convenient and permissible assumption, 
(24) w=’, 

the transformation 7’, (19) becomes 

(25) pi = p,cos8+1;w sin 6. 


The partial derivatives of G = 0 and the transformed equation G’ = 0 
are related as follows. First, 


= 
(26) G, = G. + sin dus Go 
Now, by (22), 
_ G; Pi. 
(27) — = ——sin@+/]; cos @. 
G @® 


In consequence, the sum in (26) is 


l 
(28) v(-2 sino Dv.F + cos 0 le 5 x): 


From Dp; /; = 0, and SJ? = 1, this expression vanishes. Hence 


y? 


(29) G, = G., Gy = * G, aes 
Again, 


Gy = Gj cos 0 + sin 0D > pi ole)» 


dl; 


= Gj cos 6 + © sind EG; l; + sin erase, 


, Gr’ . * 9 
= G, cos 6+ mi sin 6 cos 6 + G’ Ll, sin? 0 


the last term arising from the fact that, since © p;l; = 0, then 


Eas —h. 
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Hence, referring to (25), the final result is 


(29’) G, = G, cos 0 + g’ sin 6, 


with similar expressions for Gz and G3. 
From these equations we find at once 


(30) Gt = > GF. 
a i 


Thus the relations established in (27), (29), and (30), make it possible to write 
formula (23) in the form 
_ (GG) | 


n G, 
(31) —, ~ cos 0 — (Pi Gs) 


; = ———- gin 6 
R wl? w? G ‘ 


where 
(G; G,) = G, G, + Gs G, + G3 G,, (p1 G,) = Di G, + P2 Gy + P3 G,. 


From (31) follows the fact that the centers of the first osculating circles 
lie on a line. To see this it suffices to take as axes of rectangular coérdinates 
the normals to the elements for 6 = 0, and 6 = 90°, and to remember that 
R; cos 6, Rj sin @ are the rectangular codrdinates of the center of the oscu- 
lating circle with respect to these axes. Hence is established 

THEOREM 1. When a surface element of a partial differential equation of the 
first order is turned about an axis in its plane drawn at right angles to the char- 
acteristic, the first osculating circles for the transformed elements form a coazal 
system. 

The system of transformed equations 


(21) G'(2,Y,2, Di» Po» Ps) = 0 


may properly be called as zsogonal system from the following considerations. 
Assume S to be an integral surface of the original equation G = 0 (i.e.,6 = 0). 
Let us direct our attention upon an orthogonal trajectory D of the charac- 
teristics upon the surface S. Under 7’, each surface element of S along the 
curve D is turned through the angle @ about an axis tangent to D. Assumea 
fixed value for 6. The turned elements satisfy an equation of the isogonal 
system and form a surface band intersecting S along D under the angle 6. 
By the Cauchy existence theorem this band will lie on a unique integral 
surface of the isogonal system. Then it is clear that any given integral 
surface of G = 0 is cut under a given arbitrary constant angle @ by a one- 
parameter family of integral surfaces of the isogonal system (21), the curves 
of intersection on the given surface being the orthogonal trajectories of the 
characteristics on that surface. This result may be expressed in 

THEOREM 2. Any integral surface of an isogonal system will be cut by a 
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one-parameter family of integral surfaces of that system under a given arbitrary 
constant angle along orthogonal trajectories of the characteristics. 

Consider the special case when the characteristics of the original equation 
G = 0 are lines of curvature on all integral surfaces. The condition therefor 
is found by expressing that the direction on an integral surface conjugate to 
that of the characteristic is at right angles to the characteristic. This conju- 
gate direction (common to all surfaces containing the characteristic*) is 
given by (7) when 6 and ¢ are the denominators of dp and dq, respectively, 
in (10). In homogeneous coérdinates it is readily found that, for this common 


conjugate direction, 

(32) ba: by : bz = po G, — ps Gy: ps Gz — 1 Gz: pr Gy — po Ge. 
The condition sought, namely, G, 6a + G. dy + G; 6z = 0, reduces to 
(33) a, Gz + a2 Gy +0a;G, = 0. 


But each term in (33) is unchanged by 7,. Hence follows 

THEOREM 3. In the isogonal system G’ = 0 if the characteristics are lines of 
curvature on all integral surfaces for any value of 0, then this is true for all values. 

The class of partial differential equations of the first order for which the 
characteristics are lines of curvature on all integral surfaces forms therefore 
an invariant system under 7’',, and has been studied from this point of view 
by Liebmann.t Theorem 1 takes a special form in this case (as proved by 
Liebmann), since the first osculating circles are now great circles on the 
principal spheres of one system. ‘The paper of Liebmann deals with these 
spheres only. 

In a plane isogonal system the integral curves may be arranged in isogonal 
nets, or two one-parameter systems such that each curve of one system will 
intersect all those of the other under the same angle 6. We inquire if cor- 
responding double families of surfaces exist among the integral surfaces of 
the isogonal system G’ = 0. Let S be an integral surface of G = 0. For 
convenience, denote the family (Theorem 2) cutting S by S,. Consider those 
characteristics on the surfaces S, which intersect one of the characteristics 
on S. The locus of these «! characteristics is a surface orthogonal to S, 
and hence arises a system of surfaces T orthogonal to S intersecting it along 
the characteristics, and also containing the characteristics on the S,. If 
now S belongs to a one-parameter family So of integral surfaces of G = 0 
each of which stands in the same relation to the system S, as S, then the T 
and S, surfaces must intersect everywhere orthogonally. Furthermore, this 

* See p. 524. 

+ Aequitangential- und Isogonaltransformationen der partiellen Differentialgleichungen D,2, 
Rendiconti del Circolo Matematico di Palermo, vol. 29 (1910), 


p. 139. 
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fact is to be independent of @. Consider, in particular, 6 = 90°. The T, 
the S,, and the So surfaces now form a triple orthogonal system and hence 
intersect along lines of curvature. That is, a necessary condition is that the 
characteristics shall be lines of curvature on the Sp» surfaces. But in this 
case the characteristics are lines of curvature on the S, for all @ (Theorem 3). 
Since the 7 and S, surfaces for every @ are orthogonal, it appears that the 
characteristics on the S, must be lines of curvature on the 7’ surfaces for 
every 0, that is, the 7’ surfaces must be plane or spherical. A necessary con- 
dition for the aforesaid arrangement is thus readily seen to be that all integral 
surfaces shall be developable or canal surfaces. The isogonal net of plane 
curves is, therefore, not duplicated for the case under discussion in any general 


sense. 
Certain points in the preceding argument have excluded the linear equation 
(34) Xp+Y¥q-Z=0, 


in which X, Y, Z are functions of z, y, z, only. The isogonal system arising 
here is 


(35) Xp+ Yq -—Z+sin 0 V1 + 7° +f VV24 24 27 = 0, 





in which the primes on p and g have been dropped. The elementary cones in 
(35) for a given value of 6 are congruent cones of revolution whose axes are 
the tangents of the characteristics of (34). Equation (35) may be written 
in the alternative form 

! 4 bd if 9 
(36) I; : > i cos? 6, 
the left-hand member being the sum of the squares of the three determinants 
of the matrix. For @ = 90°, real integral surfaces exist only if the charac- 
teristics of (34) form a normal congruence, and in that case, there are a simple 
infinity of such surfaces only. 


3. THEOREMS FOR AN EQUITANGENTIAL SYSTEM 
In the differential equation 
(9) F(x,y,2,p,q) =0 


we may inquire as to the locus of the points of the surfaces elements having a 
common plane 
37 z=pr+quyte. 

This locus is plainly the curve* arising by regarding (9) and (37) simul- 
taneous equations with p and q constant. We may call this curve the ele- 


* Caratheodory (I. c.) calls this curve “die Elementplankurve.” 
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mentary curve of the plane. What is the direction of the tangent of this curve? 
The infinitesimals dz, dy, and dz at a point on the curve satisfy 


(38) F,dx+ F,dy+ F,dz =0, dz = pdx + qdy, 
from which the ratios dx : dy : dz are found to be 
dx dy dz 
(39) a 
—F,-—qF. F.+ pF. qF:—-— pF, 


This direction in the element (2, y, z, p, q) is precisely the second of the 
two directions referred to in § 1, namely, that direction on all integral surfaces 


containing the element which is conjugate to the direction of the characteristic,— 





the common conjugate direction, as we shall call it. 

It is well now to place in juxtaposition the following statements which 
appear dual. 

(a) The planes of the elements satisfying F = 0 and having a common 
point envelop a cone, the elementary cone of the point. The elements of the 
cone are tangent to the characteristics through the point. 

(b) The points of the elements satisfying / = 0 and having a common 
plane lie on a curve, the elementary curve of the plane. The tangents of this 
curve are the common conjugate directions in the plane. 

These facts justify us in assigning to the tangents of the elementary curve 





of a plane a réle in the geometry of partial differential equations equal in 
importance to that of the tangents of the characteristics. This we do in the 
following developments. 

Let (¢1, ¢2, ¢3) be the direction cosines of the common conjugate direction. 
Then by (32), . 

Pi Pe Ps 

(40) C1 3C23C3 = Ge Gy Ge’ 
the notation meaning that the c’s are proportional to the corresponding 
second order determinants of the matrix of the right-hand member. 

In § 2 we defined a turn of a surface element of a partial differential equa- 
tion of the first order. We now define a slide. 

Derinition. A slide of a surface element satisfying a partial differential 
equation of the first order is a translation through a given distance in its plane 
in a direction perpendicular to the common conjugate direction. 





Applied to a given equation 
G(x, Y,2, Pr, pe, ps) = O, 


a slide S, through the distance a transforms the elements into those satis- 
fying a new equation 


(41) G’ (x’, y’, 2’, Diy Po, Ps) = O. 
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The transformed elementary curve in any plane is obviously a parallel 
curve of the original elementary curve. In fact, it is precisely this phase of 
the matter which suggested the direction of the slide. An immediate result 
is expressed in 

THeorEeM 4. Under a slide the common conjugate direction in an element 
remains parallel to itself. 

If (01, o2, 73) are the direction cosines of the line along which the point 
of the element moves, then the equations of the slide S, are 


(42) zr’ =x+ a0, y’ =y+ao2, z’ =2z2+ 03. 


For a given equation G = 0, the slides S, form a group. 

It is our purpose to prove a theorem for space analogous to Scheffers’ result 
already quoted for plane equitangential systems. To do this a second oscu- 
lating circle must be associated with a given pair of united surface elements. 
Consider a developable surface upon which these elements lie. The direction 
of the generator of this developable is the common conjugate direction (c). 
From equations (3), § 1, we find the relation* 


(43) aB (rt — s°) = de — (ad + Be)s. 
For a developable surface, rt — s? = 0, ort 
be 
(44) = OB + be 


Putting this value of s in equation (5), we find for the curvature of normal 
sections of the developable, 

(45) | 1 (6 6a + € dy)? 

5 — - — oe ae Se ee 
j R, V1 fe P+” (ad + Be) ds? 

Apply this to the direction (¢). For this purpose we merely have to set 
a; = 61/58, o2 = bdy/és. Replacing in (45) the values of a, B, 4, e€, ete., 
in homogeneous coordinates as in § 1, simple reductions lead to the result 
(46) 1 (a; G, + 02 Gy + 03 G, y 

) Dp = een? ie . - Y i y ce Y / y ad 

Rz w (Gy; G, a Gs G, + Gs G.) . 
giving the radius R, of the second osculating circle. 

It will be well to describe this second osculating circle. The integral surfaces 
along a characteristic have in common a characteristic surface band (in Lie’s 
terminology). This characteristic band will lie on a developable, the char- 
acteristic developable. The second osculating circle for an element (2, y, z, 


* See the paper by the author already cited. 
{1 We assume a8 + de + 0, the meaning of which is that the characteristic and common 
conjugate directions in the element (xz, y, z, p, q) are distinct. 
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p,q) of the characteristic band is the circle osculating the normal section of 
the characteristic developable which is at right angles to the generator. The 
first osculating circle is, moreover, the osculating circle of this developable 
in a normal section along the characteristic. In fact (45) reduces to (8) 
when 6x :dy = a:8. The two osculating circles are therefore associated 
with this developable. Moreover, when the characteristics are lines of curva- 
ture, the two osculating circles are identical. For the generator of the develop- 
able is now at right angles to the characteristic. If ¢ is the angle between 
the characteristic and common conjugate directions, then, by Euler’s Theorem, 
R, = R, sin’ ¢. 

Return now to an element (2, y, 2, 1, P2, Ps) of the differential equation 
G = 0. Consider the second osculating circles of the «! transformed ele- 
ments under S,. These circles lie in a plane (perpendicular to the conjugate 
direction (c)), and touch the line which is the path of the point of the trans- 
formed element. But these ~! circles also touch a second line, as we now prove. 

Under the transformation 


, 


(42) z’ =2+a0\, y =yt+aaz, 2 =z+403, 


with p; : po: ps remaining unchanged, the given equation G = 0 is trans- 
formed into the system 


(41) r(x’, y’, 2, Pi, Po, Ps) =. 
For this equation the curvature of the second osculating circle is, by (46), 


1 CG +n & Gi)? 
46’) 1 _ (61 Ge + 02 Gy + 93 Gy) 





The partial derivatives of G = 0 and G’ = 0 under (42) are related as follows: 


= , _ 7” Oo; ” 00% 7” 003 
(47) G, = Gn + a( a. + Gy =~ 18! On ). 
The directions (c), (a7), and (A) (equation (17)), are mutually at right 
angles. Assume 
0, Gq G3 


(48) Ai Ae A3;= +1. 


;}C, C2 C3 | 


Then the application of (15) gives 





Mt Ae A3 

00; Ci Co C3 | (¢ =1,2,3). 
(49) yy =— “| Ae, Oc2 des | 
|\Ox Ox Ox) 


Denoting the determinant in the right-hand member by A,, equation (47) 
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becomes 
(50) G, = G, — aA, (Gre, + Gye, + Gres). 


Similar equations hold for G, and G,. 
Now ¢: Gz + ¢2 Gy + ¢3 G, = 0 by (40). Hence from (50) will follow the 
equation 


(51) (GL (+ Gy’ C2 + G,’¢3) (1 - a(cy An + Co A, + C3 A,)) = 0, 


which must hold for all values of a. In consequence, the first factor must 
be zero, and hence G,, G,, and G, remain unaltered by a slide. The same fact 
was observed for the turn S, (equations (29)). Next, 


’ yr 0o71 ’ 002 ’ 003 
52 = @G@ hs tes —+64.— }. 
(6 ) Gr +a ( G5 + 6,5 + * Opi 
Referring again to (15), we may write 
~ 0c; 
(53) Set = AN + Bes, 
in which A and B are determinants. Moreover, A = — o;/w, when equations 


(17) are used. We thus find the final form 
(54) Gi = 14+ 3 (rGz)a, 


and similar equations for G; and G3, the same abbreviation (p; G,) being 
used as in (31). 
Substitution in (46’) gives 


ay - —— (01 Gz)” _— 

aie Ry ~ @(Gi Gz) + a( pr Gz) (01 Gz) Je” 
the parentheses standing for sums of three terms as before. But this equa- 
tion shows that the centers of the second osculating circles lie on a line. For 
the codrdinates of their centers referred to the line (a) and the diameter of 
the circle for a = 0 as axes, are R} and a, and, in these codrdinates, (55) is the 
equation of a straight line. Hence 

THEOREM 5. When a surface element of a partial differential equation of the 
first order is slid in a direction at right angles to the common conjugate direction 
in the element the center of the second osculating circle describes a line. 

The system of transformed equations 


(41) G’ (2’, y’, 2’, Pi, Po, ps) = O 


may aptly be called an equitangential system. For let S be an integral surface 
fora = 0. Under S, each surface element of S is displaced a distance a in a 
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direction at right angles to the common conjugate direction. Consider there- 
fore on S the system of curves (2) which are the orthogonal trajectories of 
the system conjugate to the characteristics.* Then each element of S slides 
along a tangent of a (>) curve. Draw the system on S conjugate to the 
(=) curves, and consider one of these curves,—call it y. Circumscribe to S 
along y the tangent developable ['. The surface elements of S along y 
become under S, for any arbitrary given a a surface band along a curve Ya 
upon the tangent developable T such that the corresponding surface elements 
of S and of this band have a common distance of length a. This band along 
Ya determines a unique integral surface S, of the system G’=0. This 
surface and S have a common tangent developable and the distance on any 
generator of this developable between the points of contact with S and S, 
is constant and equal to a. That is, S and S, are equitangential surfaces. 
Accordingly we have 

THEOREM 6. Associated with any integral surface in an equitangential system 
is a one-parameter family of integral surfaces each of which is an equitangential 
surface of the original surface. The curves of contact of the common tangent 
developables are conjugate to the orthogonal trajectories of the sys:em conjugate to 
the characteristics. 

In the preceding section it appeared that the partial differential equations 
of the first order for which the characteristics are lines of curvature on all 
integral surfaces form an invariant class under 7,. A corresponding theorem 
is true when turns are replaced by slides. For the condition (33) applied to 
(41) is 

a, Gz + a2 G, + a3G, = 0. 
Now a = 133 — I$ \2 = 123 — 12, ete., and the /; are proportional to 
the G; of equation (54), say 1; = p G;. Then 


, , , a ' : 
a, = p (Gers — Gs de) + <3 (pi Ge) (o2%s — 3X2). 


Hence 
(56) a; Gz + az Gy + a3 G. = p’ (a1 G, + a2 Gy + a3G.), 


yrs 


lp yo v2 / yw2 2 2 ’ 
where p’ = VG, + G3 4+ G; VG, + Gy + G3. We may therefore state a 
. rT’ . 

companion theorem to Theorem 3 as 

TueorEeM 6. In the equitangential system G’ = 0 if the characteristics are 
lines of curvature on all integral surfaces for any value of a, then this is true for 
all values of a. 

This result was established by Liebmann in the paper cited above. 

From formula (55) it appears that Rj vanishes for one value of a, that is, 

* The point is here that the system conjugate to the characteristics takes the réle assumed 
by the characteristics in § 2. 
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(G, G,) = 0 for the corresponding surface element. In deriving (45) excep- 
tion was made of the case ad + Be = 0, namely, when the common conjugate 
direction and the characteristic direction coincide. This exceptional case 
occurs, however, normally for an element of the transformed system, the 
vanishing of (G; G,) being precisely this case. To include all cases, there- 
fore, we may assume in (46) that (G,G,) may vanish, and no flaw in the 
argument developed will be occasioned thereby. 


4, GENERAL THEOREMS 

The result established in § 2 as stated in Theorem 1 suggests an investigation 
along the following lines. 

In the differential equation 
(12) G(x, Y,2, Pry Po, ps) = O, 
the first osculating circles for an element (2, y, 2, p1, Po, p3) and the ele- 
ments derived by turning it form a coaxal system passing through a common 
second point P(#, 7, Z). Are these circles also first osculating circles at P 
for a second isogonal system of differential equations? ‘This query is the more 
natural since Scheffers has shown, in the memoir cited above, that the cor- 
responding theorem holds for a plane isogonal system. 

Begin by determining the direction of the axis of the coaxal system. Let 
C be the angle between this axis and the normal to the surface element 
E:(2, y¥, 2, Pi, P2, pz) at P. Obviously C is the value of 6 in (31) when 
R; becomes infinite, and hence 


ia w (1; Gz) 
§ C.= ; 
(57) tan (p: Ge) 


For convenience, denote the first osculating circle for the element E by ceo. 
The radius J, of this circle is given by 
1 _ (hG;) 
R, w(G)* 
Since the circle ¢p is to be the first osculating circle at P for an element 
E :(#,9, 2, pi, pe, ps) of an equation 


(58) G(#,9, 2, pi, Be, Bs) = 0, 


simple considerations show that the relations between the codrdinates of E 
and EF are 
pi = p; cos 2C + wl; sin 2C (¢ =1,2,3), 


R : Rh. 
(59) Z=2+—(fi-m), G=y¥+—(h— Pe), 


R 
2=2+— (ps — ps). 
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In other words, equation (57) results from the given equation G = 0 by 
the transformation (59). 

Express now the partial derivatives of the original equation in terms of 
those of the transformed equation (58). Omitting details, the relations are as 
follows. 


(60) G.=G6,+- “Ris H + 2C,(11 91) — sin 2C (pi a2 a32) K, 


and similar equations for G, and G, resulting when the subscript x has been 
replaced by y, etc. Moreover, 


H = (71 Gz) — (p1 G;) = 2 sin C(wcos C(h G,) — sin C (pi G;)), 
(61) = 26. 4G, I= 6, + Gs, l= 26, + Gi, 
gi = — p; sin 2C + wl; cos 2C (i =1,2,3), 
K = (ahh) = Cay G;) + (a Gh), 


and (p @ @3,) is a determinant with the usual notation. The partial deriva- 
tives with respect to p,, p2, and ps3 are expressed by equations of which the 
following is typical. 


)» 


R, 
G; = a+(4 ), H + 2C;(Iig,) — 2 sin? CI, 


(62) 


+sin 20( 2 (41) - 2 (pits) ~ (pravanr) K ). 


Consider now the necessary conditions attached to the assumption that co 
is the first osculating circle of the element E of the differential equation 
G = 0. Using, as above, a superior bar on symbols associated with G=0, 
then, since the characteristic in E must be tangent to ¢ at P, we must have 


(63) 8) = |; = sin 2C — 1; cos 2C =— g;/w (4 =1,2,3), 


in which (G) = + YG? + G2 + G3. Then from (61), 


I =» 
(hg) === »( 1,G,) = Fe [ 2G ey + ee | 


(G) (G 
~ - (GG) (Bs _ B) 
(G) w w }* 


Hence the necessary condition (J; gi) = 0, since R; must equal R,. 


(64) 
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Recognizing next that the angle C must be the same whether associated 
with G = 0, or G = 0, we have, comparing with (57), 


J= 
r 


<= tanC. 


This condition leads easily to H = 0. To complete the discussion, multiply 
equation (62) by a;, and form the sum (a;G,). Using the conditions 
H = (Iig:) = Oalready found, and remembering that (a; G,) = (a1 G:) =0, 
we find that K = 0. Then from equations (61) we see that (a; G-) = @, 
and since, fr.m (59), G, = G., also (a,;G,) = 0. But this condition holds 
when and only when the characteristics are lines of curvature on all integral 
surfaces, as shown in (33), Conversely, assuming (a;G,) = 0, we may 
readily show that (a1 G) = H = (I,g,) = K =0, and hence this is a 
sufficient condition. The details in this demonstration may be omitted since 
the sufficiency condition is implicity demonstrated in the paper of Liebmann 
referred to above. Thus we have 

THEOREM 7. The first osculating circles for a pencil of corresponding surface 
elements of an isogonal system G’ = 0 have the same relation to a second isogonal 
system at their second common point of intersection when and only when the 
characteristics are lines of curvature on all integral surfaces in both systems. 

For an equitangential system the discussion proceeds along similar lines 
with a similar result. Recalling that the configuration arising by continuously 
sliding an element is a flat straight surface band, we may state 

THEOREM 8. The second osculating circles for a flat straight band of cor- 
responding elements of an equitangential system have the same relation to a second 
equitangential system along their second common tangent when and only when 
the characteristics are lines of curvature on all integral surfaces in both systems. 

When an element of G = 0 is turned, a relation between the radii of the 
first and second osculating circles may be derived in the form 

1 _ (a1 Gz)? 


5 ou 2)" 
(65) RR t+ l= w* (G?) ’ 


and / is an invariant of the turn 7,. Similarly, when the element is slid, 


(66) R, = R2+ - . = i oF ; 
te (o1 G, ) 
and L is invariant under S,. The conclusion may be drawn at once that the 
center of the second osculating circle describes a conic when the element is 
turned, a similar property holding for the center of the first osculating circle 
when the element is slid. 
In the preceding discussion the derivatives G,, G,, and G, remained un- 
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changed under both transformations 7, and S,. The interpretation of this 
fact possesses geometric interest. When an element is turned, the common 
conjugate direction moves in a plane whose normal has the direction estab- 
lished by the ratios G,: G,:G,. When an element is transformed by S,, 
the line of centers of the second osculating circles has this direction. 

One final detail may be mentioned. When an element is slid, the charac- 
teristic direction turns about a fixed point in the plane of the element. That 
is, the tangents to the characteristics form a pencil. Denote the vertex of 
this pencil by V. Direct the attention upon an element at the point P. 
When this element is turned, the line of centers of the first osculating circles 
will pierce the plane of the element in a point W. It may be shown without 
difficulty that P is the mid-point of V and W. 

SHEFFIELD ScientiFic ScHOOL 
May, 1917 





ERRATA, VOLUME 18 
Page 73. T.H. Hitpepranpt. On a theory of linear differential equations in 


general analysis. 
Page 79, line 25, the expression 


> Cik Yokj (x) 
&=) 


should be replaced by 


b York (a) Ckj - 
k=1 
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